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Preface 


In  this  thesis  I  have  sought  to  create  a  computer  model  of  the 
NATO  standard  continuously  variable  slope  delta  voice  encoding  system 
with  sufficient  flexability  tq,  permit  continued  study  of  the  standard's 
specifications  and  tolerances.  This  investigation  has  started  the  pro¬ 
cess  of  evaluating  the  proposed  NATO  standard,  however,  additional  study 
is  necessary  to  determine  the  standard's  adequacy  to  assure  system 
interoperability . 

I  wish  to  ‘thank  my  thesis  advisor,  Capt.  Kizer,  and  the  members  of 
the  thesis  committee,  It.  Col.  Carpinella  and  Capt  Seward,  for  their 
assistance,  guidance,  and  tolerance  during  the  course  of  this  project. 

Jeffrey  A.  Lersch 
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Abstract 

A  computer  model  of  the  continuously  variable  slope  delta  voice 
encoding  system  specified  in  the  draft  STANAG  on  "Analogue/Digital  Con¬ 
version  of  Speech  Signals  for  Tactical,  Digital,  Area  Communications 
Systems",  dated  June  1978,  is  developed  and  implemented  in  FOkTRAN  IV, 
The  model's  performance  is  then  characterized  in  terms  of  idle  channel 
noise,  total  harmonic  distortion,  intermo dulation  distortion,  signal- 
to-noise  ratio,  and  frequency  response.  For  each  of  these  attributes, 
the  system's  performance  is  presented  graphically  and  compared  to  the 
criteria  established  in  the  draft  standard.  The  model  is  then  exercised 
by  varying  the  system  parameters  to  the  limits  imposed  by  the  standard 
and  the  resulting  performance  compared  to  the  previously  determined 
ideal  system  performance.  The  results  show  that  the  performance 
characteristics  measured  are  most  sensitive  to  the  primary  integrator 
response  and  output  filter  response  when  the  system  parameters  are 
restricted  to  the  range  allowed  by  the  draft  MATO  standard. 


INVESTIGATION  OF  CONTINUOUSLY  VARIABLE  SLOPE 
DELTA  (CVSD)  i-'O  D  U  L  A  TO  R  /  D  E  I-'.O  D  U  L  A  TO  R  COEPATAblL  ITY 


I.  Introduction 

A  draft  NATO  standard  on  the  analog  to  digital  conversion  of 
speech  signals  using  continuously  variable  slope  delta  (CVSD)  modula¬ 
tion  is  presently  being  circulated  among  the  military  services  for  com¬ 
ments,  The  proposed  standard,  "The  Analogue/Digital  Conversion  of 
Speech  Signals  for  Tactical,  Digital,  Area  Communications  Systems," 

June  1978,  seeks  to  assure  transmission  systems  interoperability  by- 
standardizing  the  system  architecture  and  setting  tolerances  on  key 
system  parameters.  The  draft  standard  (see  Appendix  A)  oonsists  mainly 
of  end-to-ena  system  performance  criteria,  primarily  signal-to-noise 
ratios  and  amplitude  response  characteristics.  No  standards  are  specific¬ 
ally  established  for  transnission-end/reception-end  mismatch  performance. 

The  Air  Force  Communications  Command,  AFCC/OA,  has  questioned 
whether  the  limited  number  of  specifications  given  are  adequate  to 
assure  system  performance  when  the  CVSD  encoding  equipment  is  not  per¬ 
fectly  matched  to  the  decoding  equipment.  Are  the  tolerances  specified 
sufficiently  narrow  to  assure  no  serious  signal  degradation  when  the 
modulator  and  demodulator  parameters  differ  by  the  maximum  amount 
allowed  by  the  draft  standard?  This  is  the  question  that  this  inves¬ 
tigation  seeKS  to  answer. 

Problem  Statement  Determine  the  adverse  effects  on  the  transmitted 
signal  and  their  severity  when  the  CVSD  encoder  and  decoder  parameters 
differ  within  the  limits  allowed  by  the  draft  STANAG  on  "The  Analogue/ 
Digital  Conversion  of  Speech  Signals  for  Tactical,  Digital,  Area  Com¬ 
munications  Systems,"  June  1978, 

Approach  The  approach  of  this  investigation  is  to  perform  a  computer 
simulation  of  the  CVSD  analog  to  digital  conversion  system  then  evaluate 
the  system's  performance  under  varying  external  and  internal  conditions. 
Initially,  a  basic  mathematical  analysis  of  the  system  components  is 
performed  and  mathematical  models  of  the  CVSD  encoder,  decoder  and  the 
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input  and  output  filters  are  developed.  These  models  are  then  translated 
into  computer  subroutines  and  coded  in  FORTRAN,  In  the  next  section,  the 
tests  used  to  characterize  the  model  are  described.  These  tests  consist 
of  the  standard  voice  frequency  measurements  as,  idle  channel  noise,  total 
harmonic  distortion,  intermodulation  distortion,  signal -to-noise  ratio, 
and  frequency  response.  The  system  is  first  characterized  with  the  encoder 
and  decoder  parameters  matched.  Each  test  is  performed  at  frequencies  and 
amplitudes  across  the  normal  active  range  of  the  system.  After  system 
performance  under  ideal  conditions  is  established,  the  system  parameters 
are  allowed  to  vary  across  the  ranges  allowed  by  the  draft  standard 
and  the  degradation  of  the  transmitted  signal  by  encoder  and  decoder 
parameter  mismatching  evaluated.  The  results  of  the  testing  are  then 
analyzed  to  determine  which  parameter  mismatches  most  seriously  degrade 
system  performance  and  to  determine  if  the  degradation  is  serious  enough 
to  prevent  signal  transmission . 
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1 1 .  Analog/Digital  Conversion  System  Model 

The  basic  analog/ djp.it.al  CVSD  system  defined  in  the  draft  stan¬ 
dard  is  shown  in  the  block  diagram  in  Figure  1.  It  consists  of  four 
major  components,  the  input  and  output  filters,  the  CVSD  encoder,  and 
the  CVSD  decoder.  Each  of  these  components  is  discussed  in  the  following 
sections . 

CVSD  Encoder  Operation  The  CVSD  encoder  structure  is  shown  in  Figure 
2,  The  bandlimited  signal  from  the  input  low-pass  filter  is  applied  to 
one  input  of  the  comparator  and  sampled  at  the  clock  rate,  either  16  or 
32  kb/s.  Each  input  sample  is  compared  to  an  estimate  of  the  signal 
generated  by  the  encoder  feedback  network  from  previous  input  samples. 

In  this  model,  the  comparator  output  is  +1  if  the  input  sample  is  greater 
than  the  signal  estimate  and  -1  if  the  input  sample  is  less  than  the 
estimate.  This  polar  signal  is  converted  to  binary  (+1  =1,  -1  =0)  and 
forms  the  transmitted  data  signal.  To  generate  the  next  signal  estimate, 
the  polar  signal  from  the  comparator  is  routed  to  the  input  of  the  slope 
overload  detector. 

Slope  overload,  as  defined  for  this  system,  is  the  condition  when 
the  last  three  comparator  outputs  are  identical,  either  all  +l's  or 
all  -l's.  This  indicates  that  the  input  signal  amplitude  is  rising  or 
falling,  respectively,  faster  than  the  encoder  can  track  using  the  pre¬ 
sent  step  size.  Other  systems  define  slope  overload  by  different  length 
strings  of  identical  comparator  outputs.  Strings  of  two  or  four  identical 
bits  are  also  commonly  used  to  indicate  slope  overload.  The  last  three 
comparator  outputs  are  stored  in  a  shift  register  within  the  slope  over¬ 
load  detector  and  combinational  logic  circuits  used  to  determine  if  a 
slope  overload  condition  exists.  The  slope  overload  detector  output  controls 
the  position  of  the  switch  shown  in  the  block  diagram.  Under  normal 
conditions,  when  slope  overload  does  not  exist,  the  switch  is  in  the 

V  position.  Upon  occurrence  of  slope  overload,  the  switch  position 
min 

is  changed  and  V  is  applied  to  the  input  of  the  syllabic  filter, 
m&x 

The  syllabic  filter  is  generally  a  simple  single  pole  RC  filter 
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Figure  1.  The  CVSD  Signal  Processing 


whose  output  is  defined  as  the  step  size  of  the  CV3D  encoder.  The  syllabic 

filter  controls  the  response  of  the  system  to  the  envelope  of  the  speech 

sipe-1  being  processed.  Prolonged  anpli  ration  of  V  .  to  the  input 

min 

of  the  syllabic  filter  causes  the  output  to  decrease  to  a  minimum  non¬ 
zero  value  that  approaches  V  .  .  Under  continuous  slope  overload  con- 

min 

ditions,  V  is  continuouslv  aprlied  to  the  input  of  the  syllabic  filter 
max  ” 

causing  the  filter  output  to  increase  to  an  average  value  approaching 

V  .  The  magnitude  of  the  syllabic  filter  output  is  used  to  control  the 
max 

amplitude  of  the  output  pulse  of  the  pulse  moduli'. tor .  Polarity  of  the 
pulse  is  controlled  by  the  last  output  of  the  comparator. 

The  primary  integrator  responds  to  the  square  w-ve  signal  from  the 
pulse  modulator  and  its  output  forms  the  signal  estimate  used  by  the 
comparator .  At  the  end  of  each  clock  period  a  new  estimate  is  avail¬ 
able  to  be  used  by  the  comparator  in  generating  the  next  binary  output  ’ 
and  the  next  signal  estimate.  The  primary  integrator's  response  con¬ 
trols  the  maximurr  analog  signal  frequency  that  can  be  processed  Trough  the 
the  CVSD  analog,  to  digital  conversion  system .  In  figure  3,  are  shown 
sample  waveforms  at  each  stage  of  the  analog  to  digital  conversion 
process . 

F.ncodor  Algor'  ■  hm  The  mathematical  description  of  the  CVSD  encoder  oper¬ 
ation  is  largely  a  description  of  its  component  filters,  the  primary 
integrator  and  the  syllabic  filter.  One  of  the  system  characteristics 
specified  by  the  draft  standard  is  the  primary  integrator  response.  The 
impulse  response ,  in  its  simplist  form,  is  given  as, 

a(t)  =  e“2irfclt  (1) 

where 

f  ,  =  the  ool<  frequency  of  the  filter  in  hertz 

cl 

The  primary  inte~.ro tor  innut  signal  is  the  square  vave  output  of  the 
pulse  modulator,  which  for  a  single  pulse  can  be  described  as, 

a  (t )  =  0  t  <  0  and  t  >  T 

=  a  0  £  t  *  T 


vJ 


slope 

pulse  syllaDic  overload  .  primary 

binary  modulator  filter  detector  comparator  integrator  analo; 

output _ output _ output _ output  output  output  mput 


sample  rate 


Figure  3.  Sample  Waveforms  at  Various  Points  in  the  Encoder 


where 

T  = 


sample  rate 


The  primary  integrator  output  is  determined  by  convolving  the  filter 
impulse  response  with  tne  input  signal. 


c 

x(t)  =  a  ( t. )  *  Qf(t)  =  J“ a  (  t)  -  r)  dr 

—  oo 

=  0  ,  for  t  <  0 


r  -2n f  ti 

[l-e  cl  J  ,  for  0  S  t  S  T 


27Tf 


cl 


-27Tf  T,  -27Tf  (t  -  T) 


(3) 


2irf 


a  [l-e  cl  ]  e  ‘"“cl"'  x '  ,  for  t  >  T 
cl 


Since  the  primary  integrator  output  iu  of  interest  only  at  the  end  of 
each  clock  period,  when  it  is  used  for  comparison  with  the  input  analog 
signal,  the  continuous  equations  developed  above  can  be  simplified 
as  follows.  For  t  =  nT, 


r2-  (i  -  a)  an  ”  1  ,  n  =  1,2,  ...  ,11 

TL 


(41 


where 


-2nf  ,  T 
a  =  e  cl 


k.  =  2irf 
1  cl 


Using  suoorpoaiti on ,  t^e  primary  integrator  output  as  the  result  of  a 
series  of  N  Dulses  can  be  described  as, 


><1 


—  (1  -  a)  +  afI_t  [l  -  a) a  +  ...  +  ^  (i  -  a)aU  L] 


IT— 1 


Eri-n 
k1 


^  n  (1  -  a)  a"  ,  for  n  -  1,2, 


•  ,  u  (3) 


n=0 
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This  expression  ran  also  be  defined  recursively ,  depending  only  on  the 
present  input  and  the  last  output.  This  definition  can  be  used  to 
simulate  the  CV3P  encoder  on  a  computer. 


=  Via  +  (1  -  a)  ^ 


(6) 


The  analysis  of  the  syllabic  filter  output  follows  identically 
that  of  the  primary  integrator.  The  impulse  response  of  the  syllabic 
filter  is, 


0(t) 


rPrr  t 

U 


] 


(7) 


where 

t  =  the  time  constant  of  the  syllabic  filter 
c 

The  recursive  expression  for  the  syllabic  filter  output  is, 

b,  -  b.-E’  ♦  (1  -  «  ^  (8> 

2 \* 
t 

c 

exp  [-  —  ] 
c 

either  V  or  V  .  ,  the  syllabic  filter  input 

max  mm 

Encoder  Computer  Subroutine  Equations  (6)  and  (8)  are  implemented  in 
the  subroutine  used  to  perform  the  CVSD  encoding  for  this  investigation. 
Figure  3  is  a  flowchart  of  the  subroutine  used  for  encoding  and  Appendix 
B  is  the  FORTRAN  code  used.  All  of  the  system  defining  parameters  are 
transmitted  to  the  subroutine  through  the  calling  statement.  Encoding 
is  performed  on  an  array  basis.  The  analog  signal  to  be  analog  to  digital 
converted  is  first  sampled  at  the  clock  rate  and  the  samples  placed  in 
the  input  array,  which  is  of  size  1  x  N,  where  N  is  the  number  of  samples. 
All  of  the  samples  are  encoded  by  the  subroutine  and  the  binary  data 
stream  placed  in  the  output  array  before  the  subroutine  returns  control 


where 

k. 


Figure  A.  CVSD  Encoder  Subroutine  Flowchart 


LO 


Figure  A  (continued ) .  CVSD  Encoder  Subroutine 

Flowchart 
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to  the  calling  program. 

Two  of  the  primary  system  defining  parameters ,  VT'AX  and  VAIN 
must  be  generated  by  subroutine  VMAXOPT  (appendix  I)  before  the  encoder 
subroutine  is  cailed.  It  should  be  noted  that  trie  constants  and 
derived  in  equations  (5)  and  (8)  are  not  specifically  included  in 
the  program  statements  defining  the  primary  integrator  and  syllabic 
filter  responses  but  are  expected  to  be  included  in  the  values  calculated 
for  WAX  and  WIN. 

The  variable  DC  in  the  subroutine  is  the  duty  cycle  of  the 
slope  overload  detector.  This,  variable  is  not  used  during  the  encoding 
process.  Instead,  it  is  used  only  by  VHAXOPT  when  the  values  of  VI 'A X 
and  WIN  are  being  determined. 

CVSD'  Decoder  Oreration  The  CVSD  decoder  circuit  is  identical  to  the 
encoder  feedback  circuit,  A  block  diagram  of  the  decoder  is  shown  in 
figure  5.  The  only  difference  between  the  decoder  and  the  encoder  is 
that  the  decoder  has  no  comparator.  The  binary  signal  from  the  encoder 
is  applied  directly  to  the  slope  overload  detector  and  the  output  signal 
is  taken  from  the  primary  integrator.  The  signal  estimate  generated  in 
the  decoder  is  identical  to  that  generated  in  the  encoder,  if  the  parameters 
of  each  unit  are  i  dr-nv  ical .  However,  at  the  decoder  the  signal  estimate 
is  of  interest  at  all  times  and  not  just  at  the  sample  periods,  as  the 
decoder  signal  estimate  is  the  approximation  of  the  analog  signal  trans¬ 
mitted  by  the  CVTD  encoder.  Figure  6  shows  sample  waveforms  at  various 
points  within  the  decoder.  The  waveforms  are  identical  to  those  shown 
in  figure  3,  except  that  the  decoder  primary  integrator  output  is  shown 
as  a  continuou  ,  signal. 

Decoder  Algorithm  Since  the  decoder  circuit  is  identical  to  the  "Header 
circuit  without  the  comparator,  the  mathematical  analysis  level o: e  i  fur 
the  encoder  is  also  applicable  to  the  decoder.  One  exception,  however, 
is  that  the  simplification  used  to  obtain  equation  (4)  is  not  -/morally 
applicable  to  the  decoder  since  the  primary  integrator  output  '  ihe 
decoder  is  required  to  be  continuous.  The  recursive  expression  for  the 
decoder  output  is, 
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sample  rate 


i 


Figure  6.  Sample  Wave  forms  at  Various  Points  in  the  Decoder 


-  [y(<  [N-:]7)]  e"1!1  + 


(1  -  e  "1  )  (9) 

kl 


for  (N-l )T  <  t  <  NT 


where 


sample  rate 


=  the  primary  integrator  input  for  (M-l)T  <  t  <  NT 
=  the  pole  frequency  of  the  primary  .in  to  ora  tor  it)  hertz 


Analysis  of  the  decoder  syllabic  filter  is  identical  to  that  of  the 
encoder  syllabic  filter  and  equation  (.8)  also  applies  to  the  decoder. 


An  An-1^  +  (1  ^  k2 


where 


i?  =  —  ,  An  =  the  syllabic  filter  output 


p  =  exp 


t 

c 


V.,  =  either  V  or  V  .  ,  the  syllabic  filter  input 

N  max  run 


Decoder  Computer  "ubroutin- '  Using  equations  (6)  v-ii  (8),  the  decoding 
subroutine  is  implemented  as  shown  in  the  flowchart  in  figure  ",  Equation 
(9)  is  not  used  since  the  straight  line  approximation  provided  by  the 
".llcomp  plotter  rrovides  a  sufficiently  accurate  representation  of  the 
decoder  output  for-  this  investigation.  Except  for  the  elimination  of 
the  comparison  step  used  in  the  encoder  subroutine,  the  decoder  sub¬ 
routine  is  nearly  identical  to  that  of  the  encoder.  All  comments  applic¬ 
able  to  the  encoder  subroutine  are  also  applicable  to  the  decoder  sub¬ 
routine.  The  r'ORTRAU  code  for  the  decoder  subroutine  is  attached  in 
Appendix  C. 
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Figure  7,  CVSD  Decoder  Subroutine  Flowchart 
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Pncodor  and  Pncodcr 


From  the  expressions  developed  in  the 


preceding  sections  describing  the  CYSD  encoder  and  decoder,  is  can  he 
seen  that  there  are  four  parameters  that  determine  the  characteristics 
of  the  encoder  and  decoder.  They  are,  f  for  the  primary  integrator, 
t  for  the  syllabic  filter,  and  V  and  V  .  whose  value  determines  the 
magnitude  of  the  step  sizes, 

Tlie  draft  standard  specifies  the  value  of  f  ^  explicitly  in  para¬ 
graph  3,2,  When  the  primary  integrator  consists  of  a  single  pole  filter, 
the  value  of  f  is  required  to  he  between  100  and  300  Hz. ,  Other  poles 
and  zeros  can  be  added  to  the  primary  integrator,  in  accord;  nee  with  the 

draft  standard,  however,  only  f  is  required.  In  this  investi gation , 

cl 

the  single  pole  version  of  the  primary  integrator  is  used  in  the  CVSC 
encoder  and.  decoder  models. 

For  the  syllabic  filter,  the  draft  standard  does  not  specify  the 

value  of  t  directly.  Instead,  t  is  specified  in  terms  of  the  decoder 
c  c 

output  signal  when  a  given  input  is  applied  to  the  encoder.  When 
the  analog  input  signal  at  300  Hz  is  stepped  from  -42  d3m0  to  0  dBmO, 
the  decoder  output  signal  is  required  to  achieve  90%  of  its  final  value 
within  2  to  4  milliseconds  after  the  output  signal  starts  to  rise.  ("0TF : 
For  this  system,  the  standard  specifies  the  reference  test  level  point 
to  be  -4  dEr.i .  So,  a  -42  dlmO  is  actually  -46  d3n .  All  measurements  taken 
in  this  investigation  are  stated  in  dBmO,  unless  explicitly  stated 
otherwise . ) 

The  values  of  V  and  V  ,  are  also  not  specified  directly  by  the 
max  min 

draft  standard,  but  are  specified  in  terms  of  the  syllabic  filter  output. 

The  syllabic  filter  output,  which  has  previously  been  defined  as  the 
step  size  of  the  encoder  and  decoder,  is  required  to  ue  linear  as  a 
function  of  the  slone  overload  detector  duty  cycle.  The  slope  overload 
detector  duty  cycle  is  defined  the  ratio  of  the  number  of  tines 
slope  overload  is  detected  to  the  number  of  samples  in  the  sane  period. 

In  paragraph  3,4  of  the  draft  standard,  the  step  size  is  shown  us  varying 
linearly  as  the  duty  cycle  ranges  from  0  to  .5.  The  step  size  ratio,  the 
ratio  of  the  syllabic  filter  output  when  an  800  Hz,  0  dBmO  sign-  1  is 
applied  to  the  encoder  input,  to  the  syllabic  filter  output  when  the 
encoder  input  is  grounded  is  required  to  be  34  dB  ±  2  dB.  This  specification 
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in  combination  with  the  specifications  for  f  and  t  determine  the  values 

cl  c 

of  V  and  V 

max  min 

Due  to  the  fact  that  the  parameters  interact  with  each  other,  the 

values  of  t  ,  V  ,  and  V  .  need  to  be  determined  recursively.  A  value  of 
c  max  min  J 

f  ,  is  chosen  within  the  range  given  by  the  standard  and  an  estimate  of 
cl  " 

t  chosen  near  its  expected  value.  The  syllabic  filter  determines  the 

system  response  to  the  amplitude  modulation  of  a  voice  signal.  As  the 

highest  frequency  in  the  envelope  of  the  voice  signal  is  generally  about 

100  Hz,  t  is  estimated  to  be  the  reciprocal  of  this  frequency  or  .01. 

A  nominal  step  size  ratio  is  given  by  the  draft  standard  to  be  34  dB . 

These  three  parameters  are  used  to  calculate  the  values  of  V  and  V 

max  min 

Figure  8  is  the  flowchart  of  the  program  that  calculates  these  values 

using,  the  subroutines  shown  in  figures  9  and  10,  then  plots  the  resulting 

syllabic  filter  output  as  a  function  of  slope  overload  detector  duty  cycle. 

Initially,  estimated  values  of  V  and  V  .  are  used  and  the  slope 

max  min 

overload  detector  duty  cycle  and  system  step  size  ratio  calculated  when 

an  800  Hz,  0  dBmO  test  signal  is  input  to  the  CVSD  encoder.  If  the 

calculated  values  are  not  within  the  tolerances  specified,  V  and  V  . 

max  min 

are  adjusted  and  the  calculations  repeated.  This  process  is  continued  un¬ 
til  values  of  V  and  V  .  are  determined  that  produce  a  slope  overload 
max  min 

detector  duty  cycle  of  .5  ±  1%,  and  a  step  size  ratio  within  . 01 %  of  the 
input  value , 

After  determining  the  values  of  V  and  V  .  ,  the  entire  CVSD 

max  min 

system  is  tested  to  determine  if  the  rise  time  requirement  is  met  using, 

the  parameters  that  have  been  calculated.  The  flowchart  of  the  test  program 

is  shoun  in  figure  11.  To  determine  the  system  rise  time,  a  test  signal 

consisting  of  alternate  series  of  500  samples  of  an  800  Hz,  -42  dBmO  sine 

wave  and  500  samples  of  the  800  Hz  sine  wave  at  0  dBmO,  The  initial  series 

at  -42  dBmO  initialize  the  storage  elements  of  the  slope  overload  detectors 

in  both  the  encoder  and  decoder  and  get  the  system  into  an  initial  steady- 

state  condition.  After  processing  the  test  signal  through  the  system, 

the  output  signal  is  plotted  in  the  vicinity  around  one  of  the  steps  in 

input  signal  power.  The  system  rise  time  is  then  determined  graphically. 

Figure  12  shows  a  sample  output  from  this  program  for  both  the  16  and  32 

kb/s  sample  rates.  This  test  was  performed  with  f  ,  =  100  Hz,  step  size 

cl 

ratio  =  34  dB,  t  =  .01  for  the  16  kb/s  sample  rate,  and 
c 
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Figure  13 


Input-,  and  Outgirt  I_ow-PasP  i  1  hr r s  The  last  of  the  najor  continents  making 

up  tlie  CVSD  analog/ digital  conversion  system  are  the  input  and  output  fil¬ 
ters.  These  filters  are  used  to  limit  both  the  input  and  output  sig¬ 
nal  spectrum  to  the  voice  band  frequencies  only.  For  telephonic  communi¬ 
cations,  the  voice  band  is  generally  considered  to  be  those  frequencies  less 
than  3600  Hz,  lor  optimal  system  performance,  these  filters  should  have 
a  very  sharp  cut-off  and  high  loss  characteristics  in  the  stop  band. 

The  purpose  of  the  input  filter  is  to  limit  the  input  signal  spec¬ 
trum  to  prevent  aliasing  due  to  the  sampling  process.  When  the  input 
signal  is  sampled,  in  addition  to  the  input  spectrum,  the  output  spectrum 
also  contains  sum  and  difference  frequency  components  centered  around  the 
sample  frequency.  If  the  input  spectrum  were  to  contain  frequencies  very 
much  larger  than  the  desired  spectrum,  aliasing  or  interference  would 
occur  v/hen  the  difference  frequencies  fell  into  the  baseband  spectrum. 

For  this  model,  the  input  filter  is  considered  to  be  an  ideal  low-pass 
filter.  The  test  signal  generator  output  spectrum  is  limited  to  the  voice 
band  frequencies  only,  with  no  components  falling  outside  that  range.  This 
simulates  a  low-pass  filter  with  zero  insertion  loss  in  the  pass  band  and 
infinite  Iocs  in  the  stop  band. 

The  function  of  the  output  filter  is  also  to  limit  the  signal  to 
the  voice  band,  however,  in  this  case,  the  components  outside  the  original 
input  spectrum  are  produced  by  the  non-linearit ies  of  the  processing 
system.  The  output  filter  smooths  the  signals  and  eliminates  the  harmonic 
components  above  the  voice  band.  This  filter  may  have  the  same  character¬ 
istics  as  the  input  filter  or  may  have  a  narrower  pass  band  to  improve 
performance.  The  filter  chosen  for  this  model  is  a  maximally  flat,  linear 
phase  symmetrical  finite  impulse  response  (FIR)  filter.  The  mocel  of 
this  filter  was  developed  by  J.F.  Kaiser  of  the  Digital  System  Research 
Department  of  the  Bell  laboratories.  Reference  1  provides  more  complete 
documentation  of  the  filter  model.  There  are  two  parameters  that  define 
the  response  of  the  filter,  beta  and  gamma.  Figure  14  shows  the  response 
of  the  filter  generated  by  this  program  and  where  beta  and  gamma 
are  defined.  Beta  is  the  normalized  center  frequency  of  the  transition 
region  and  gamma  is  the  normalized  width  of  the  region.  Normalization 
is  with  respect  to  the  sample  rate.  This  filter  was  chosen  for  its  flat 


response  in  the  p  i:-:.;  bind  in  order  to  minimize  disturbance;  of  the  OVSD 
r  nco dt ■  r  / deco  :o r  response  ,  since  those  arc  the  nrirory  system  cor.ponon  ta 
under  investigation .  The  parameters  chosen  for  the  filter  are  /}  =  .1875 
and  y  -  .1  for  the  10  kb/s  sample  rate  and  (3  -  .1  and  y  =  .1  for  the  32 

kb/s  sample  rate.  Table  ii  lists  the  filter  coefficients  generated  by  the 
program  .and  figure  IS  shows  the  frequency  response  for  both  filters. 


Figure  14  .  Maximally  Flat  FIR  Fi  lter  Response 

Character; Stic  and  Parameter  Definition 

The  center  of  the  transition  region  for  the  16  kb/s  filter  is  3  kHz  and 
3,2  kHz  for  the  3?  kb/s  filter.  As  can  be  seen,  the  maximally  flat  char¬ 
acteristic  is  achieved  at  the  expense  of  stop  band  loss.  However,  it  will 
be  shown  in  the  performance  results  that  the  system  performance  meets 
moot  of  the  criteria  specified  by  the  draft  standard  in  spite  of  the 
poor  filter  performance. 

Filter  Subroutines  The  filter  program  developed  by  J.F.  Kaiser  is 
used  to  generate  the  FIR  filter  coefficients,  however,  it  has  been  mod¬ 
ified  to  be  a  subroutine  that  returns  the  coefficient  values  to  the  calling 
program  instead  of  printing  them  out.  These  coefficients  are  produced  by 
FI  TRGF.H  then  used  by  subroutine  FIITER  to  actually  filter  the  si>  nal  input 
to  the  filter.  The  maximum  number  of  coefficients  that  can  be  r ro d- iced 
by  FI-TRGEH  without  program  modification  Is  200,  Subroutine  FlI'i'FR 
delays  the  output  signal  by  200  sample  period:;  so  that  it  has  at  least 
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TABI.E  II 


FIR  FILTER  COEFFICIENTS 


16  kb/s  snnnle  rate 


P( 

i) 

Z 

.32*')  -  9  2926 

R( 

2) 

Z 

.29:- 2  7.3  32^1 

B  ( 

3) 

z 

.1  3"  T  656  2 

P( 

4) 

z 

-  .  7 S T F  9477  b 

B  ( 

5) 

z 

'.657791 724 

R  ( 

A) 

z 

-.11836  75375 

B  ( 

7) 

= 

.  ,237i  35732 

B  ( 

3) 

z 

. .  2  377  p-iir 

B  ( 

9) 

= 

.  ‘j  )  ’3  3^  '<■  29 

P< 

1C) 

z 

1  27.  7  R  9264 

B  ( 

11) 

z 

-.  V?',  4 1 c  -  2 

P( 

12) 

= 

.  ;  j  21  ?  7-r3'  7 

P( 

13) 

z 

.  j  •"  -*  3  -“3621 

R  ( 

14) 

= 

.r  ' 135 948!  9 

B( 

lO 

r 

-.  '  o  124  22<  2? 

E< 

13) 

z 

0  3'  390  5 

P( 

17) 

z 

-. .  1  ,19770... 

P( 

IP) 

z 

.00  4  32956  1 

P( 

19) 

r 

.  27  519  26 

P( 

20 

z 

•  1.)  5  26  2  4 

P( 

21) 

z 

31339 

5( 

2?) 

z 

-.  .  1  1-9  :0  1 

P( 

2  3) 

= 

,  ^  ’  1  %9il 

P( 

24) 

z 

.7  ?  0175639 

P( 

20 

= 

.  .  n  6?8<  •: 

P( 

23) 

= 

146  44 

V( 

2^) 

r 

-  •  j  i  '»  216  2  6 

B  ( 

2^) 

z 

-.  ■'  )  '7  ,  *16521 

B  ( 

2-') 

- 

.  '  )  .'C. 014  29 

P( 

30 

- 

o  nr  "ie45 

B  ( 

3 11) 

z 

.:-7  IOC  7136!: 

P( 

32) 

r 

''r  -;,u7 

P( 

33) 

z 

- .  1  7  ; 1 '  7  o  C  c  q 

P( 

34) 

= 

:)r  110  36 

E  ( 

35) 

= 

-  .  *  1  "i  1  '  i .  ** 

P< 

36) 

= 

.  );) •  )  jl  2 

B( 

37) 

z 

. litn'i 

P( 

39) 

r 

.c  n  j :  i-iro 

B  ( 

39) 

z 

.  -  I*'  )0L.  j 

P( 

4C) 

z 

,  r  1*;1'  O'.,  7 

e< 

41) 

- 

.  ,)J7  9)10. 

p< 

42) 

z 

.  J  3  )  J  MKU 

32  kfc/  n  rar.plc  rate 


r  / 

1) 

Z 

.194584523? 

B( 

?) 

r 

•  13'' 345 g£  67 

R( 

3) 

z 

.1  427570  h 

9( 

4  ) 

z 

.  .,9  23'  2  21'  26 

O 

3) 

z 

.2423271121 

P( 

5) 

= 

143332322 

P( 

7  ) 

= 

1651  9  16-3 

P( 

9) 

r 

-..22 7 7843?7 

D( 

O 

z 

-.Cl 9759612 3 

P( 

1  O 

- 

-.11  031  021 

P( 

11) 

z 

-.i)3:  ■>  6  77  ’  2 

P  ( 

12) 

= 

.  6  ’  17 3112  1 

n  ( 

I7) 

r 

.19  35  2  3 p4  '9 

B( 

14) 

z 

..  j  3  2  T  3  8  2  T  3 

P( 

1  5 ) 

r 

.  7  217  3:20 

B( 

10 

z 

.-l’-  H33  3 

P( 

17) 

= 

.*  K2"066  35 

D( 

15) 

z 

-.  .  1  0  21  )£2  2 

D( 

19) 

z 

-  .  ■  j .  2  3  1  39  7  , 

?( 

20 

z 

0  2  '4  3956 

°( 

21) 

r 

! 1 O 38 275 3 

9  ( 

22) 

r 

- .  :  : ?  8 « i  -7 

P( 

23) 

= 

-.Oil  j  393114 

B  ( 

24) 

z 

-.t  )  7  31-06  4  3 

B( 

20 

- 

- . O  696  5  3 

E  ( 

26) 

= 

-  .‘  1  'i  I  2  612  8 

e( 

27) 

z 

,'T  )3 z:  2 

B( 

26) 

= 

-.  •*  V')  7  2LZS 

9< 

29) 

z 

-  9  264  7 

P( 

30 

z 

.  )  !1 5  6 

F  ( 

31) 

z 

:  J  6  3  4 

R  ( 

72) 

z 

-.  /I  O  -  9%.  7 

P( 

33) 

z 

- .  n  ’  ■  7 . 6  o :  i 

E( 

34) 

= 

-.  •  i  .  j  .  O'O  Q 

Q< 

35) 

= 

-.  i  .)  ■:  7  m  o 

B  ( 

36) 

z 

-  ,  :  1  '  J  o  C 4  *4 

B  < 

37) 

z 

- . ;  *  0  3  r  P )  i '  c. 

p  ( 

30 

z 

To:  o 

P( 

30 

z 

.  :oo  9  7 . o 

9  ( 

4  0) 

z 

.  r  .1  9  1  '  r* U  u 

?( 

41) 

z 

.  i  TJ  ) ;  vjUrj 

P( 

42) 

z 

.  '1  J  J  -  )  w  ^  ^ 

^  a 


D 


FREQUENCY  (HZ) 


Figure. 15b,  CVSD  System  Model  Output  Filter  Response 
for  32  kb/s  Sample  Rate 
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ure  16.  Signal  Filtering  Subroutine  Flowchart  (FILTER) 
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Figure  16,  (continued) 
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200  input  signal  samples  can  be  used  by  the  filtering  algorithm.  liqua¬ 
tion  (10)  shows  the  filtering  expression  implemented  by  the  FI  IT EH 
subroutine . 


f 


the  nth  output  sample 

the  nth  input  sample 

the  ith  FIR  filter  coefficient 

the  number  of  filter  coefficients 
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III. 


Performance  Tests 


A  model  oT  the  continuously  variable  slo.<e  d.-l!  en.ilog  to  dig¬ 
ital  conversion  system  is  constructed  from  the  component  models  described 
in  the  previous  sections.  Figure  17  shows  tbr  test  configuration  sim¬ 
ulated  by  the  computer  model  used  in  this  investigation.  The  system 
under  test  is  shown  in  figure  1.  In  this  simulation,  the  test  signal 
generator  is  a  subroutine  that  generates  samples  of  a  sinusoidal  signal 
that  can  be  composed  of  up  to  two  frequency  componcnts  at  individually 
specified  amplitudes.  The  standard  test  signal  used  in  the  performance 
tests  is  an  800  Hz  sine  wave  at  -20  dBnO,  unless  otherwise  stated.  As 
previously  indicated,  the  reference  signal  level  is  -4  dBm.  All  power 
measurements  are  made  relative  to  this  level.  The  test  signal  is  generated 
as  an  array  of  5000  samples  for  most  of  the  tests  performed.  This  array 
is  then  processed  through  the  system,  the  output  array  of  each  system 
component  becoming  the  input  array  of  the  next.  The  final  system  output 
signal  is  then  processed  to  determine  the  various  signal  characteristics. 
System  performance  is  measured  in  terns  of  the  commonly  used  voice  frequency 
tests  as,  idle  channel  noise,  total  harmonic  distortion,  intermodulation 
distortion,  signal-to -noise  ratio,  and  frequency  response.  These  tests 
are  first  performed  with  the  CVSD  encoder  and  decoder  parameters  matched, 
then  performed  again  with  various  combinations  of  encoder  and  decoder 
parameters  to  show  ho w  system  performance  degrades  under  mismatched 
conditions . 


Idle  Channel  Noise  Test  Idle  channel  noise  is  measure  of  the  basic  amount 
of  noise  that  the  processing  system  adds  to  the  output  signal.  The  output 
signal  power  is  measured  while  the  system  input  is  grounded.  Any  non¬ 
zero  power  measured  is  the  idle  channel  noise.  Before  measuring  the  idle 
channel  noise,  however,  the  system  insertion  loss  is  first  set  so  that  the 
standard  test  signal  experiences  no  change  in  power  after  being  processed 
through  the  system.  Idle  channel  noise  is  then  measured  as, 

ft 

IC"-  111  > 

n=l 
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Figure  17.  Simulated  System  Test  Configuration 


Figure  18.  Idle  Channel  Noise  Program  Flowchart 
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where 


the  output  signal  amplitude 
the  number  of  samples 


y 

N 


n 


Figure  18  is  the  flowchart  of  the  idle  channel  noise  test  used  to  measure 
the  CVSD  system  performance. 


Total  Harmonic  Distortion  Test  Total  harmonic  distortion  is  one  of  the 
measures  of  system  non-linearity.  The  CCITT  procedure  for  measuring 
total  harmonic  distortion  is  to  input  a  single  frequency  test  signal 
near  the  center  of  the  system's  pass  band  and  measure  the  magnitude  of 
the  harmonic  components  in  the  output  spectrum.  Total  harmonic  dis¬ 
tortion  is  then  calculated  by, 


THD 


V- 


E3  * 


+  E, 


x  100% 


(12) 


where 


E2'  E3*  *’*  ,EN  =  of  the  harmonic  signal 

components  in  the  output  spectrum 

=  the  RMS  voltage  of  the  primary  signal  component  in  the 
output  spectrum 

N  =  the  largest  harmonic  within  the  system  pass  band 


Figure  19  is  the  flowchart  of  the  computer  program  used  to  calculate 
total  harmonic  distortion.  After  processing  the  single  frequency  sine 
wave  test  signal  through  the  CVSD  system,  the  output  signal  spectrum 
is  calculated  using  a  fast  fourier  transform  (FFT)  .  Due  to  the  limitations 
of  the  FFT,  the  standard  test  signal  is  not  used,  instead,  a  1000  Hz 
signal  at  -20  dI3m0  is  used.  The  FFT  procedure  can  only  measure  sig¬ 
nal  components  at  multiples  of  the  minimum  frequency  resolution  which  is 
determined  by  the  number  of  samples  in  the  FFT  window.  In  this  case,  the 
window  length  wa s  specified  to  be  255  samples,  which  allowed  a  frequency 
resolution  of  62,5  Hz  at  the  16  kb/s  sample  rate  and  125  Hz  at  the  32  kb/s 
sample  rate.  A  1000  Hz  test  signal  was  chosen  as  being  both  compatible 
with  the  FFT  and  a  commonly  used  test  signal  in  voice  frequency  measure¬ 
ments  . 


Figure  13.  Tota1.  Harmonic  Distortion  Ten'  ProrrE.:; 
Flowchart  (THD) 
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In‘^r:--iril  ■■  1  iori  Oj_l  *or‘  ^Q''1  The  total  harmonic  di  short  ion  test  often 

does  not  five  a  complete  idea  of  the  system  response  non-1  i near i tier. . 
Intermodulation  distortion  is  another  measure  of  non-linearity  used  in 
voice  frequency  system.  The  CCITT  procedure  of  measuring  interr.odula tion 
distortion  is  to  input  a  composite  test  signal  made  of  of  two  sinusoidal, 
signals  of  equal  amplitude.  The  frequencies  of  the  two  cisnals  are 
separated  by  an  amount  that  the  difference  frequency  is  within  the 
pe.os  band  of  the  system.  In te modulation  distortion  is  t,l,en  calculated 
fcy» 


whe  re 

F.  =  the  RMS  voltap.e  of  she  di  -  I v.r,- ■nc»-  :  r*-o  •  <v  :■ 

dii, 

in  the  output,  spectrum 

E  =  the  RMS  voltap.e  of  the  •  ir-s:  fre  :  y  v  ■  •  <  •  ■  • 

output  spectrum 

Eg  =  the  RMS  vo  1  tare  of  tic  second  ‘me -mien  <>  •.<■■  ■  ■  • 
output  spec+rum 


Fif.ure  2.0  shows  the  flowchart  of  the  pro r ram  d  t 

modulation  distortion  for  the  CTIT  system .  The  :  ;i, 
similar  to  the  total  harroni  c  distortion  ; ax*;  r 
tral  components  used  from,  the  KK7  out;).’  are  •  •> 

and  the  difference  frequency.  The  tort  mi  nal  . 
sists  of  7b 0  Hz  and  1000  Hz  compon^m  * . ■  a*  -  ., 


Signal-to-Ho  i re  Ratio  Measurement  Th<-  i  ,-m,u  2 -to-noi  rn”; 
measure  of  how  accurately  the  system  be  in*-  characteri 
input  signal ,  A  test  signal  is  process  through  tie  sjvt-er  a 
ing  output  sipnal  compared  to  the  input  signal  after  eo:--.;.-n:-n '  .  • 
the  system  insertion  loss  and  ni final  delay.  SHR  is  then  cal  cm’ 
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Iii  i  t  i  a  1  i  ze 


Figure  20.  IntermoOulation  rintortion  Test  IYo;'i\> 
riowchurt  (  IHTI::::-TD ) 
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n-1 

X 

Y) 

■■■  the 

nth  input  r: 

liyna.l  sample 

"  n 

the- 

nth  outnut 

sir.ua1.  sample¬ 

V 

=  tlie 

total  nur.be 

s'  o’"  sampler. 

Fis'.urv  21  r’nows  the  flow  d.-tr  t  for  the  signal -to-noi se  ;  ror.ran  used  to 
characterize  the  f’VSD  uyr performance .  The  standard  BOO  Hz  test 
signal  is  used  to  per  Torn,  the  initial  syste:  characterization . 


Frequency  i  --c.w  'O  -i-.sre -'-’•it  Two  neiiwlr  of  performing  frequency 
response  measure-. on t  are  used  in  this  invest!. e  tion .  "he  first  is  flat 
weighted  measure:,  en !  which  is  used  to  determine  the  frequency  response  of 
the  entire  CVSP  system  sine  this  is  the  me thod  for  which  the  draft  stan¬ 
dard  specifies  per  for-:. nee  criteria.  A  second  method  is  the  frequency 
selective  measurement  of  the  response  character! sti os.  This  method  i: 
used  to  investigate  the  frequency  response  of  the  CVSD  encoder  and  decoder 

only, 

rial  weighted  frequency  response  measurement  is  performed  by  input- 
tinp  a  sin: le  frequency  sine  wave  tost  sifnal  at  a  constant  amplitude,  then 
measuring  the  r-yr.tem  output  signal.  lower.  The  output  signal  includes 
coi  iponontr,  at  frequent  ion  other  than  the  test  signal  frequency,  however, 
the  power  of  ti,  •  entire  composite  signal  i  r,  measured  without  filtering. 

The  measured  pain  variations  are  then  ploi.ed  and  scaled  such  that  the 
BOO  Hz  measurement  is  0  ri3 .  The  flowchart  of  thrt  rrORrain  using  this 
procedure  if;  shown  in  figure  02. 

Frequency  selective  mea statement  of  frequency  resp-nse  uses  the 
s  «me  test  ;  i-ocedure  except  only  the  magnitude  of  the  output  signal  com¬ 
ponent  at  the  test  fre  quency  is  nnasum  .i.  T!  e  o  -.her  comro’w  ri  ts  of  the 
conrof  It”  output  signal  are  riot  included  in  this  measurement ,  The 
measurement;';  are  then  scaled  and  riot4 ted  such  that  the  1000  Hz  measure¬ 
ment  is  0  *53.  Figure  23  is  the  flowchart  of  the  program,  to  perform  the 
free  ••'•ncy  selective  meas.urrrr  rnts .  T!ie  1000  ’.\z  measurement  is  used  as  the 
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Figure  22 


(continued)  Pror.rnm  DGAIN  Flowchart 


Figure  22.  (continued)  Program  DGAITI  Flowchart 
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5 


reference  vain**  since  the  fast  Courier  transform  used  to  calculate  the 
output  signal  spectrum  cannot  measure  the  component  at  800  Hz. 


1 
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IV.  Test  Results 


The  results  of  the  tests  described  in  the  previous  section  are 
presented  here.  Each  test  was  first  performed  with  the  standard  800  Hz 
test  signal,  while  the  CVSD  encoder  and  decoder  parameters  were  matched. 
This  test  characterized  the  ideal  system  performance  with  the  system 
parameters  at  their  nominal  values.  Next,  the  tests  were  performed  allowin 
the  system  parameters  to  vary  across  the  ranges  shown  in  Table  I  and  using 
test  signals  that  varied  in  frequency  and  power  across  their  normal 
dynamic  ranges,  while  still  maintaining  encoder/ decoder  match.  Finally, 
the  test  were  repeated  again  with  the  encoder  parameters  held  constant 
at  one  extreme  of  the  permissible  values  and  the  decoder  parameters  allowed 
to  range  to  the  opposite  extreme.  Each  test  was  performed  changing  one 
variable  at  a  time  while  the  other  were  held  at  their  nominal  values. 

Idle  Channel  r'oise  The  results  of  the  idle  channel  noise  tests  are 
shown  in  figure  24.  For  each  sample  rate,  the  idle  channel  noise  perfo: — 
mance  improves  as  the  step  size  ratio  increases.  This  results  from 
the  decrease  in  minimum  step  size  as  the  step  size  ratio  increases.  The 
output  signal  depends  entirely  on  the  minimum  step  size  when  the  system 
input  is  zero  or  grounded,  fince  the  minimum  step  size  is  defined  to 
be  non-zero,  the  output  signal  will  alternate  positive  and  negative 
around  zero  attempting  to  approximate  the  zero  input  signal.  The  smaller 
the  deviation  from  zero,  the  less  the  power  in  the  output  signal  and  the 
better  the  idle  channel  noise  performance.  System  performance  exceeds 
the  criteria  specified  in  the  draft  standard.  Idle  channel  noise  is 
-88  dBmO  vs.  the  specified  -50  dBmO  at  16  kb/s  sample  speed  and  -97  dBmO 
vs.  -60  dBmO  at  32  kb/s. 

Encoder/decoder  parameter  mismatch  has  no  effect  on  idle  channel 
noise.  This  is  a  result  of  the  fact  that  no  matter  what  the  encoder's 
parameters,  the  output  will  always  be  alternatin'  ones  and  zero;'  when  the 
encoder  input  is  grounded.  Therefore,  the  input  signal  at  the  decoder 
will  always  be  the  same  and  the  output  signal  will  only  be  affected  by 
the  decoder  parameters.  The  idle  channel  noise  performance  under  mis¬ 
matched  conditions  will  be  the  same  as  shown  in  figure  24  where  the 
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Figure  24.  CVSD  Signal  Processing  System  Idle 
Channel  Poise  Performance 


i 
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parameters  are  those  of  the  decoder. 

Total  Harmonic  Pi  ntort.ion  The  draft  standard  specifies  no  maximum 
total  harmonic  distortion  for  the  CVSD  system,  however,  it  is  generally 
accepted  that  distortion  levels  of  less  than  20a  will  not  usually  be 
objectionable  to  the  system  users.  As  can  be  seen  from  the  test  results 
shown  in  f inure  2b,  system  \ier Torrance  at  the  1C  kb/s  sample  rate  exceeds 
this  limit  by  4-? System  performance  when  the  -sample  rate  is  increased 
to  3?  kb/s  improves  substantially.  The  total  harmonic  distortion  level 
droos  to  approxim.atlcy  651,  Figure  26  shows  that  the  distortion  level 
at  both  sample  rates  is  relatively  constant  for  all  input  power  levels 
within  the  normal  operating  ran"?  except  at  the  very  low  power  levels. 

When  the  encoder  and  decoder  parameters  are  mismatched,  the  total  har¬ 
monic  distortion  performance  shows  some  degradation  as  f inures  27,  23, 
and  2r<  indicate.  The  largest  .-.mount  of  deviation  from  the  matched 
system  performance  occurs  at  the  very  low  power  levels  where  the  impact 
will  have  the  least  effect.  As  figure  2o  shows,  total  harmonic  distortion 
is  most  sensitive  to  mismatches  of  the  encoder  and  decoder  primary  inte¬ 
grator  pole  frequencies,  Syllabic  filter  time  constants  and  step  size 
ratios  have  minimal  impact  on  the  system  performance  when  mismatched, 
however,  all  have  the  most  impact  at  the  very  low  input  power  levels. 

Intermodul f ion  Pi n tort  ion  Intermodulation  distortion  performance  for 

the  system  model  with  nominal  parameter  values  is  shown  in  figure  30. 

As  is  the  case  with  the  total  harmonic  distortion  test,  the  draft  standard 
provides  no  perfronance  criteria.  In  general,  intermodulation  levels 
of  more  than  4-55-  will  be  objectionable  to  a  system  user.  At  the  16  kb/s 
sample  rate,  the  intermodulation  distortion  measured  ranges  from 
1  to  5%  depending  on  the  syllabic  filter  time  constant  used.  The  dis¬ 
tortion  falls  to  approximately  3  51  when  the  sample  rate  is  increased  to 
32  kb/s.  Figure  31  shows  the  system  intermodulation  response  as  the 
input  signal  power  is  varied.  System  non-linearities  cause  the  distortion 
levels  to  rise  at  the  very  low  signal  levels  and  at  the  high  input  power 
levels.  Across  the  normal  operating  levels  between  -30  dBmO  and  -30 
dBmO,  the  distortion  is  generally  less  than  55',.  When  the  encoder  and 
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Figure  25.  CVSD  Signal  Processing  System  Total  Harmonic 

Distortion  Performance  with  hncoder  and  Pocouer 
Parameters  Matched  (Test  Signal  =  1000  Hz,  -POdr.mO) 
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gure  26a,  CVSD  System  Total  Harmonic  Distortion  Performance 
vs,  Input  Signal  Pov;er  with  Encoder  and  recoder 
Parameters  Hatched  at  16  kb/s  Sample  Rate 
(1000  Hz  Test  Signal) 
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Figure  26b,  OVSD  System  Total  Harmonic  Distortion  Performance 
vs.  Input  Si,-nal  Power  with  Encoder  and  Decoder 
Parameters  Hatched  at  32  kb/s  Sample  Kate 
(1000  Hz  Test  Signal) 
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Figure  27  a.  CVS!)  System  Total  !!nr: .onic  Cir.tortion  Performance 
vs.  Input  Si  coal  Power  with  Encoder  and  Decoder 
Stop  Size  Ratios  tisi/.y  tched  at  16  kb/s  Sample  Rate 
(1000  Hz  Test  Signal) 
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.  CVSD  System  Total  Harmonic  Distortion  Performance 
vs.  Input  Signal  Power -with  Kneoder  and  Decoder 
Step  Size  Ratios  :i snatched  at  32  kb/s  Sample  Rate 
(1000  Hz  Test  Signal  ) 


F ip, lire  27  b 


Figure  28  a.  OVSD  System  Total  Harmonic  Pistortion  Perforrnanc 
vs,  Input  Signal  Power  with  Encoder  and  Decoder 
Syllabic  Filter  Time  Constants  Si snatched  at 
16  kb/s  Sample  Rate  (1000  \h:  Test  Signal) 
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Figure  20  b,  CVSD  System  Total  !:  ir:;o’ii  -  -tor*,  ion  Pep  -n.'.- 
vs  .  Input  Signal  rov.*<-r  v •  *  h  :  ■)  o  :..-r  .-nr.  '  .  r 
Syllabic  Fil  ter  Time  Foil,  i  ;  i.c-'.c  .1  at 

32  kb/s  Sample  Rate  ( ’■  F:’.  ) 
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and  deco  Ip r  are  r.ii  smutched,  the  i n te rnodu la t i  on  distortion  perfor¬ 
mance  shows  the  came  characteristics  as  the  total  harmonic  distortion. 

The  step  sice  ratio  and  syllabic  filter  time  constant  have  a  minimal 
impact  on  system  performance  as  snow  by  the  results  in  figures  32  and 
33,  Then  the  primary  integrators  have  differing  pole  frequencies,  the 
internodulation  distortion  measurements  show  more  deviation  from  the 
matched  system  performance  as  shown  in  figure  34 . 

mi  cnal-to-f:oise  Rati  o  Figures  35  to  43  show  how  the  syst^ia  ‘<.-1 
signal-to-roise  performance  changes  with  variations  is  system  parameters 
and  input  test  signals.  The  SNR  performance  under  matched  conditions 
with  the  standard  test  signal  shows  very  little-  variation  with  differing 
values  of  step  size  ratio,  syllabic  filter  time  constant,  and  primary 
integrator  pole  frequency.  Fignal-to -noise  ratio  vs.  input  frequency 
performance  meets  the  criteria  set  by  the  draft  stnadard  across  most  of 
the  voice  band,  encoder/ decoder  mismatches  of  stop  •; •  ■'  ..tio  and  syllabic 
filter  tine  constant  have  very  little  impact  on  system  performance,  A 
nisr.atch  of  the  primary  integrator  pole  frequencies,  however,  have  a 
much  larger  effect  on  system  performance.  The  SHF.  is  degraded  below 
the  criteria  set  by  the  draft  standard,  with  the  largest  deviation  from 
the  ideal  performance  occuring  at  the  lower  frequencies .  At  the  16  kb/s 
sample  rate  the  SNR  in  lowered  by  about  5  dB  and  at  the  32  kb/ s  speed, 
about  4  dB. 

kirn  a  1-to-no  ise  ratio  perforr.ionrr  vs.  input  signal  power  fails 
to  meet  the  criteria  established  in  the  draft  standard.  At  input  levels 
below  approximately  — 10  d3m0  the  model's  performance  falls  below  the 
desired  level.  Trends  in  system  performance  as  the  result  of  variations 
in  the  system  para  ao t«-rs  can  be  observed  i...  spite  of  this  poor  .-orfoi-iance , 
a-.  :  how.  3-.  fiyur-  43  to  "3,  fetched  system  performance  minimal  change 
as  the  system  parameter-:  -cried  across  the  ranges  specified  in  Table 

I,  When  the  encoder  and  decoder  are  not  matched,  SNR,  like  tO-:l  harmonic 
distortion  and  inter,  odulat  ion  distortion  ,  si  lows  litll-  via  l  ■  >  •  .  f;x>  > 

1- matched  siv.  performance  except  at  the  very  low  signal  lev  Is, 

SNR  shows  the  most  c’acs-r-  from  ideal  system  performance  when  the  primary 
integrator  nole  frequencies  are  mismatched. 
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Figure  32b.  C.VSD  System  Intermodulation  Distorti  on  fcribrsance 
vs.  Input  Signal  Power  with  Kn coder  and  r .coder 
Step  Size  Ratios  Mismatched  at  3?  kb/s  S-ivele  Rate 
(750  and  1  (>00  flz  Test  Signal) 
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Figure  33d  .  CVSD  System  Intern©  dulati on  Distortion  Performance 
vs.  Input  Signal  Power  with  encoder  and  Decoder 
Syllabic  Kilter  Time  Constants  Mismatched  at  16  kb/s 
Sample  Rate  (750  and  1000  rlz  Test  Signal) 
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Figure  33b,  CVSD  System  Intermodulation  Distortion  Performance 
vs.  Input  Signal  Power  with  Encoder  and  Decoder 
Syllabic  Filter  Time  Constants  Mismatched  at  32  hb/s 
Sample  Rate  (750  and  1000  Ha  Test  Signal) 
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Figure  34a.  CVSD  System  Intermodulation  Distortion  Performance 
vs.  Input  Signal.  Power  with  Encoder  and  Decoder 
Primary  Integrator  Pole  Frequencies  .Mismatched 
at  16  kb/s  Sample  Rate  (750  and  1000  llz  Test  Signal) 


Figure  34b.  CVSD  System  Inte modulation  Distortion  Performance 
vs.  Input  Signal  Power  with  Encoder  and  Decoder 
Primary  Integrator  Pole  Frequencies  Mismatched  at 
32  kb/s  Sample  Rate*  (7bC  and  1000  1 1/.  Teat  Signal) 
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Figure  35  ,  CVSP  Signal  Processing  System  Signal-to-Poiue 

Performance  with  Encoder  and  Decoder  Para: >oters 
Matched  (Test  Signal  -  800  Hr. ,  -SOdBmO) 
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CVRD  System  Ripnal-to-iloise  Performance  vs. 
Frequency  with  Encoder  and  Decoder  Parameters 
Matched  (-20  dPir.Q  test  signal)  at  16  kb/s 
Sample  Rate 


CVRD  System  Signal-  to  -  •  *6  i  ;3<?  Per  f'orimn'-t-  v.  , 
Frequency  with  Encoder  and  Decoder  Parameters 
fetched  (-20  dHmu  test  siena J  )  at  32  kb/s; 
Sample  Rate 
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Figure  37a.  CVSD  System  Si ~nal-to-Foise  Performance  vu, 
"requenry  with  Encoder  and  Decoder  Step  Size 
Ratio  "isnatehed  at  16  kb/s  Sample  Rate 
(-20  dBmO  test  signal) 
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■'igure  37b.  CVSD  System  Signal-to-floise  Perfonuncc  vs. 

Frequency  wi^h  Fncod*»r  and  Decoder  Step  Size 
Ratio  Mismatched  at  '2  kb/s  Sample  Fiatc 
(-20  dBmO  test  sign  ;1  ) 
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Figure  3g  a.  CVSD  System  Si .qna  1  -  to -’foi se  Performance  vs. 

Frequency  with  Encoder  ana  Decoder  Primary 
Integrator  Pole  Frequencies  *  Unmatched  at  16  kb/s 
Sample  Rate  (-20  dBr.O  Test  Signal) 


Figure  30  b.  CVSD  System  Sirnal-to-'ioise  Performance  vs. 

Frequency  with  Encoder  and  Decoder  Primary 
Integrator  Role  Frequencies  Mismatched  nt,  3;’  kV/s 
Sample  Rate  (-20  dlknO  Test  Sipnal ) 
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Figure  40a.  CVSD  System  Si  ma  l  -  to  -?b  i  se  Performance  vs. 

Input  Signal  Am;  lltudo  ri th  Encoder  and  Decoder 
Parameters  Matched  at  1C  kb/ s  Sample  Rate 
(800  Hz  test  ;;  i  rnai  ) 
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Figure  40Vi.  CVSD  System  Signal -to -Hoi so  Performance  vs. 

Input  Signal  Amplitude  with  Encoder  and  Decoder 
Parameters  "a tehee  at  32  kb/s  Sample  Kate 
(800  II 7.  test  signal) 
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;  rec  ;■  n.~  v  :  r-rou.  m  The  system  freq ucncy  resiO i..;e  characterist.;  an  ere 
speci  fie  ’•  in  the  <:ra  ft  s  tan  d:  ire  when  n.c-arurt-rt  by  the  flat  wci.fhtod  method. 
The  results  of  the  tests  performed  usin/r  Ur*  s  technique  ere  shown  in  fie  — 
uros  44  to  47.  Tilt  computer  model's  performance  complies  with  the  specific' 
criteria  r-xet  rt  at  the  top  end  of  the  voice  band.  The  model  ’  s  response 
roll:)  off  sharply  at  approximately  3  kHz.  V.Tiile  the  draft  standard  allows 
some  ro!  1-off,  the  response  is  not  allowed  to  brook  sharply  un’il  0 
Variations  in  the  system  parameters  have  very  little  effect  on  the  response 
characteristics  when  the  encoder  and  decoder  are  matched.  Under  mis¬ 
matched  conditions,  frequency  response  performance  follows  the  same 
pattern  as  that  established  in  the  previously  described  tests.  The  step 
size  ratio  and  syllabic  filter  tine  constant  have  minimal  i mpac  1  ,  while 
the  primary  integrator  pole  frequency  causes  increased  deviation  from 
the  matched  ays  ter.  performance . 

Usinq  tie  frequency  selective  measurement  technique,  the  response 
characteristics  of  the  CVS!;  encoder  and  decoder  connected  back-uo-Dack 
without  the  input  or  output  filters  were  measured.  The  rc  sul  ts  are  shov.Ti 
in  f if. lire s  48  to  tl  .  These  tests  show  that  the  response  rolls  off 
at  %  the  sample  frequency  at  both  the  lb  and  3k  kb/s  sample  rates.  At 
4  kHz  for  the  it  Kb/s  sample  rate  and  8  kHz  for  the  3?  kb/ s  sample  rate 
the  frequency  response  cui  ves  break  sharpls.  Incoder  sr.d  decoder  mis¬ 
match  has  practically  no  effect  on  the  frequency  response  of  these 
system,  components.  The  pr  unary  integrator  pole  frequency  shows  sliphtly 
ore  effect  on  the  response  than  the  other  parameters.  Its  effect  is 
mainlv  at  the  very  Jow  frequencies  in  the  voice  band. 
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Figure  45  a.  CVSD  Systen  Sain  Variation  vs.  Frequency  vnth 

Encoder  and  Decoder  Step  Size  Ratios  Mismatched 
at  IS  kb/s  Sample  Rate  (-20  dBnO  Test  Signal ) 


t'.VSD  System.  Sain  Variation  vs  .  Frequency  v/i  th 
Encoder  and  Decoder  Step  Si :;e  Ratios  Mi  nmatchod 
at  32  kb/s  Sample  Ha  to  (-20  dltaO  Test  Si.-nc>  1 ) 


Figure  40  b  . 
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Figure  47a.  OVSD  System  Gain  Variation  vs.  Frequency  with 
Encoder  and  Decoder  Primary  Integrator  Pole 
Frequencies  Mismatched  at  16  kb/s  Sample  Rate 
(-20  dBmO  Test  Signal) 
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Figure  47b.  CVSD  System  Gain  Variation  vs.  Frequency  with 
Encoder  and  Decoder  Primary  Integrator  Pole 
Frequencies  Mismatched  at  32  kb/s  Sample  Rate 
(-20  dDmO  Test  Signal) 
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Figure  46  a . 

CVSD  Encoder /Decoder  Hack— to— Pack  Frequency 
Response  Performance  v:i  t.h  "ntched  Parameters 
at  16  kb/s  Sample  Rate  (-20  aPmO  Tost  Signal) 

Step  Size  iiitio 

n  -  32  d:-' 

O  -  34  dB 

A  ~  30  dB 


FREQUENCY  l HZ) 

r'irure  46  b.  C VS D  Encoder/Pecoder  Bnck-to-Back  Frequency 
Response  Performnce  v/ith  (latched  Parameters 
at  3?  kb/s  Sample  Rar.e  (-2.0  dBmO  Test  Signal) 


FREQUENCY  t HZ  J 

Figure  50  b.  CVSD  Encoder/Decoder  Back-to-Back  Frequency 

Response  Perfornar.ee  v/ith  SyllaDic  Filter  T:  ne 
Constants  Mi  snatched  at  3?  kb/s  Sample  !\a  te 
(-20  dBmO  Test  Sipnal ) 
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Figure  51a.  CVSD  Encoder/Decoder  Back-to-3ack  Frequency 
Response  Performance  with  Primary  Integrator 
Pole  Frequencies  mismatched  at  16  kb/s  Sample 
Rate  (-20  dBmO  Test  Signal) 
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Figure  51b.  CVSD  Encoder/Decoder  Bad; -to— Back  Frequency 
Response  Performance  with  Primary  Integra  to  r 
Pole  Frequencies  mismatched  at  32  kb/s  Sample 
Rate  (-20  dBmO  Test  Signal) 
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V,  Cone ]  m;  i  ons  and  Recommendations 

The  test  results  show  that  the  computer  model  meets  most  of  the 
performance  criteria  set  by  the  draft  standard  when  the  system  parameter  i 
are  matched  an  a  at  their  nets  in  >1  values.  The  output  filter  docs  not 
have  the  stop  band  loss  characteristics  specified  in  the  standard  and 
as  a  result,  the  system  performance  is  marginal.  Si rnal-to-no i an  r  it: os 
fall  below  the  established  criteria  when  the  input  power  is  less  than 
-10  dS:nO.  In  spite  of  this,  the  general  effects  of  variations  in  the 
system  parameter  valuer;  and  encoder/decoder  mismatches  can  be  observed 
in  the  test  results. 

1,  When  the  encoder  and  decoder  are  matched,  chan  rc::  in  the  icy 
size  ratio,  .-yliotic  filter  time  constant,  and  the  primary  integrator 
pole  frequency  values  within  the  tolerances  allowed  by  the  draft  standard 
have  a  tie.  ■,!  i  gable  effect  on  the  transmitted  signal . 

2,  If  the  step  size  ratio  or  the  syllabic  filter  time  constant 
are  not  the  same  in  both  the  encoder  and  the  decoder,  the  effect  of  the 
mismatch  on  Vac  transmitted  signal  is  neplir.able  except  at  input  power 
levels  loss  than  -32  dBmO.  At  these  levels,  the  effects  would  not  be 
noticeable  to  the  system  users, 

3,  System  performance  is  most  sensitive  to  encoder  and  decoder 
primary  integrator  pole  frequency  mismatches.  All  the  performance  tests 
show  a  larger  deviation  from  the  matched  system  performance  when  the 
primary  integrators  are  mismatched.  This  tyro  of  parameter  mismatch 
dominates  missis tch.es  of  the  other  parameters. 

4,  The  frequency  response  of  I  .hr  r.ysV'.  is  dterm.inei  largely 
by  the  output  filter  when  the  pass  band  is  restricted  to  less  than 

’(  the  sample  rate .  Above  'A.  the  sample  rate,  the  response  is  determined 
by  the  CVS!)  encoder  and  decoder  pass  bind.  The  frequency  selective 
measurement  of  the  encoder /-decoder  r  e-ponse  siv i  that  the  syce.-:.;  is 
incapable  of  meeting  the  draft  standard  gain  variation  vs,  frequency 
criteria  rivcn  in  figure  7a  of  Appendix  A.  The  encoder  and  decoder  alone 
have  a  re.  onse  the!,  falls  off  sharply  above  A  kHz  for  the  16  kb/s  sample 
rate,  while  the  standard  requires  that  the  response  not  fall  off  more 
than  5  until  C>  kHz  is  reached. 

b.  The  specification:  in  i  tolerances  <’iven  in  the  dr.xft  standard 

i“>(  j 


.  r.'  ;  'ecua  ( ■'  to  ::'r  rare  reasonable*  system  nerforync'-  g.-ir  vc  i  c<*  -  ’  , 

Ti.p  t.f.  it'.  sni  1  1 .<• -.1  scnals  wi  1  i  s  .;«•  ;s  degradation  due  1.0  ; -a  )'a:v 

!■■  i  . !  brl  !•'  r.  between  the  encoder  ami  u>»cmler  ■  Tn  most  on  set',  l  he  &o'!"'ndi*  i 
is  minimal,  however,  .  Vi  '.ion-r  '  .tin;  would  be  necessary  to  determine 
if  the  system  resixinso  would  continue  to  meet  the  draft  standard  criteria 
un  dor  ni s mu to he d  condi tions . 

Recomr.en da tj  on s 

1.  Since  this  model's  performance  is  marginal  under  ideal 
conditions,  para  meter  mi  :.;:uteh  causes  the  performance  to  fall  below 

the  criteria  set  in  the  draft  standard.  Testin',  should  be  repeated  using 
a  filter  that  lias  higher  stop  band  loss.  The  additional  testing  should 
concentrate  on  primary  integrator  response  mismatches  between  the 
encoder  and  decoder,  since  the  other  parameters  have  little  effect  on 
the  syten  response. 

2.  System  tolerance  to  bit  errors  in  the  transmission  system 
was  not  tested,  A  mismatch  between  t'ro  encoder  and  decoder  may  cause 
increased  sensitivity  to  transmission  errors.  Testing  to  establish  the 
system  response  to  transmission  bit  error  rate  my  be  desireable. 

3,  The  draft  standard  specifications  for  gain  variation  in  the 
region  between  4  kHz  and  C  kHz  for  the  10  r.b/ s  sample  rate  should  be 
modified.  Performance  should  be  allov.’‘"-u  to  roll  off  as  sharply  as 
possible  a  love  3,5  k’V. . 

4,  Testing  was  p<-r formed  using  continuous  sinusoidal  test  signals 

only.  If  quasi-ana lo  signals  are  expected  to  be  used  on  the  CVP"  svs  ten, 
teat  in.  s.h  on’  1  be  repeated  using  this  typ  •  of  signal  a  veriou.,  input 
Veyin*  s  ,  Piece  the  .H.'"H  1  -orit'.i.  ‘  pends  on  the  high  sannle  to 

sun ole  ■  orr"  Is  sion  c:  ..nr  a'-  te"  i  sties  of  tic  ./vice  signal,  quasi- analog 
si.'n'il:.;  can  he  expected  i.o  suffer  noro  rleg.rad.it ion  as  the  re  •  <11  of 
encoder/  decode.-  s. i.seatohes ,  Ho  re  restrictive  tolerance.?  i  -.ay  be  nr-o- 

es.s  ry  to  ns.  urn  a  le...u.i  i  •»  rori’orisance  with  these’  signals .  In  a  :.<li  lion  , 
more  standard;  g  y  be  necessary ,  such  as  delay  distortion  spec i fj cations . 
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1.1  Analogue/digital  conversion  of  telephone  signals  (speech  or 

other  voice-band  signals  shall  be  performed  by  a  delta 
coder/ decoder  using  syllabic  companding  controlled  by  a 
three  bit  logic  . 
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2,  Four-wire  to  Four-wire  Audio  Frequency  Characteristics 

2 .1  Relative  Level  at  Points  A  and  3 

The  relative  levies  at  points  A  and  B  shall  be  -4  dBr. 

2.2  The  absolute  level  is  calculated  by  the  equation  dBm  =  dBr  +  dBmO. 

2.3  Impedance  at  Points  A  and  B 

The  nominal  value  of. the  impedance  at  points  A  and  B  shall  be 
600  ohms , 

2.4  Return  loss  at  Points  A  and  B  against  600  ohms 

The  return  loss  at  points  A  and  B  shall  be  a  16  dB  in  the 
frequency  range  from  300  Hz  to  3400  Hz  against  a  load  resistor 
of  600  ohms  with  an  input  level  of  -20  dBmO. 

2.5  Symmetry  at  Points  A  and  B 

Points  A  and  B  shall  be  balanced  and  not  referred  to  ground, 
i.e.  shall  be  floating. 

3,  Details  of  the  Coder  and  Decoder  Circuits 

3.1  Input  and  Output  Audio  Filters 

For  frequencies  above  6  kHz,  each  filter  shall  have  an  atten¬ 
uation  of  2:  25  dB. 

3.2  Frequency  Response  of  the  Principal  Integrator 

The  ideal  amplitude  frequency  characteristic  between  points  F 
and  G  is  shown  in  Figure  3, 


f,  =  100  Hz  to  300  Hz 

0 

N. 

1 

\si 6  d3/octave  f2  (optional )  &  2700  Hz 

dB 

N 

f g  (optional )  X  6400  Hz 

A 

\ 

12  dB/octave 

6  dB/octave 

f 

f 

f 
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Fig.  3  -  Ideal  Amplitude  Frequency  Characteristic  of 
the  Principal  Integrator 
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3.3 


3.4 


3.5 
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Modulation  Level 


A  signal  of  800  Hz  and  0  dBmO,  applied  to  point  A  of  the  coder 
shall  give  a  duty  cycle  (mean  proportion  of  binary  *1 •  digits 
at  point  D  each  one  indicating  a  run  of  3  equal  bits  at  point 
C)  of  c ,  =  0,5  at  point  D  of  the  modulation  level  analyzer 
MLA). 


Compression  and  Expansion 


In  the  coder  and  decoder  the  quantizing  step  size  q  which 
drives  the  principle  integrator  at  Point  F,  shall  have  an 
essentially  linear  relationship  to  the  duty  cycle  at  point  D 
of  the  MLA  integrator  (see  Figure  4). 


Nominal  Values 


-  Limiting  Values 

(curve  as  nearly 
linear  as  possible 

q  =  Quantizing  Step  Size 
°  at  0  dBmO 


1  +  98  c^ 
50 


Fig.  4  -  Relation  between  MIA  Output  Duty  Cycle  and 
Size  of  Quantizing  Steps 


It  follows  that  the  ratio  of  the  quantizing  step  size  at  point 
F  corresponding  to  a  aucy  cycle  of  c^  =  0.5  at  point  D  of 
the  MLA  integrator  at  the  minimum  step  size  q  shall  be 
34  dB  (provisional  tolerance:  +  2dB), 

Companding  Speed 

The  following  is  valid  for  the  condition  that  C  is  connected  to 
C' .  When  an  800  Hz  sine  ave  signal  at  point  A  is  suddenly 
changed  from  -42  dBmO  to  0  dBmO  the  output  signal  at  point  B 
shall  reach  90%  of  its  final  value  within  2  mS  to  4  r.;S, 


NOTE:- 

The  MLA  integrator  circuits  of  the  coder  and  decoder  shall 
have  the  same  characteristics  and  hence  the  same  companding 
speed. 
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Procedure  for  Testing  the  Delta  Decoder 

The  test  bit  sequence  generator  is  connected  to  the  decoder 
input  point  C'  (see  Figure  2). 

Testing  is  performed  by  means  of  periodical  test  bit  sequences 
(listed  in  Table  1)  which  result  in  audio  signals  at  800  Hz  at 
the  decoder  output  point  3.  The  800  Hz  levels  at  point  B 
shall  conform  to  the  values  given  in  Table  1 . 

When  the  signal  at  point  C'  is  switched  from  the  periodical 
test  bit  sequence  to  the  periodical  test  bit  sequence  g,  then 
the  output  signal  at  point  B  shall  reach  90 %  of  its  final 
value  within  5.5  mS  to  11.5  mS.  When  the  signal  at  point  O' 
is  switched  from  the  periodical  test  bit  sequence  g  to  the 
periodical  test  bit  sequence  a,  then  the  output  signal  at 
point  B  shall  reach  10%  of  the  value  of  the  periodical  test 
bit  sequence  g  within  4  mS  to  8mS. 

NOTE:-  for  clarification 

For  an  RC  circuit  in  the  MIA  integrator  with  time  constants  of 
4  mS  fbr  both  charging  and  discharging,  the  envelope  charac¬ 
teristic  of  the  output  signal  at  point  B  is  shown  in  Figure  5. 

For  the  case  of  switching  the  signal  at  point  O'  from  the  seq¬ 
uence  g  to  sequence  a,  the  amplitude  at  the  beginning  of  dischar¬ 
ging  is  at  the  first  moment  after  switching  higher  -  by  a  factor 
of  50  -  than  the  final  value  which  is  reached  asymptotically. 

The  final  value  equals  -42  dBmO,  i.e.  0.00794,  the  amplitude 
at  the  beginning  of  the  discharging  is  hence  0.397  (c^  =  0). 

The  value  of  10%  is  then  reached  at  5.76  raS. 


Fig.  5  -  Envelope  Characteristic  of  the  Output  Signal 
at  Point  B  (Half  the  Envelope) 
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4 .  Electrical  Performance  at  Points  A  and  B 

4.1  General 

The  required  values  under  4.2  to  4,8  are  valid  for  the  condi¬ 
tion  that  C  is  connected  to  C’  . 

For  measurement,  the  input  (point  A)  and  the  output  (point  B) 
are  to  be  terminated  with  600  ohms,  and  signals  whose  fre¬ 
quencies  are  sub-multiples  of  the  sampling  rate  shall  be 
avoided.  Accordingly,  where  a  nominfil  test  signal  frequency 
of  800  Hz  is  indicated,  the  actual  frequency  shall  be  slightly 
different;  a  preferred  value  is  820  Hz,  but  frequencies  from 
804  to  860  Hz. 

The  measurements  according  to  Sections  4.2  to  4.5  shall  be 
performed  selectively. 

4.2  Insertion  Loss  between  Points  A  and  3 

The  insertion  loss  between  points  A  and  3  at  800  Hz  with  an 
input  level  of  0  dBmO  shall  be  0  dBnO  -  2  dB.  The  insertion  loss 
contributed  by  the  transmit  and  receive  sides  shall  not  ex¬ 
ceed  one-half  of  the  value. 

4.3  Attenuation  Distortion  with  Frequency 

The  attenuation  distortion  relative  to  800  Hz  measured-’With  an 
input  level  of  -20  dEmO  applied  to  point  A  shall  be  within  the 
limits  of  Figure  G.  The  distortion  contributed  by  the  trans¬ 
mit  side  alone,  measured  at  point  G  of  the  codex’,  shall  not 
exceed  the  limits  indicated  by  the  broken  lines  in  Figure  6, 

4 .4  Variation  of  Gain  with  Input  Level 

The  deviation  of  the  output  level  compared  with  the  value  at 
-20  dBmO  shall  not  exceed  the  limits  given  in  Figure  7  for  a 
frequency  of  800  Hz, 

4.5  Idle  Channel  Noise 


Idle  channel  noise  at  16  kbits/s: 

The  idle  channel  noise  at  point  B  shall  not  exceed  -45  dBmOp. 
The  level  of  any  single  frequency,  measured  selectively,  shall 
not  exceed  -50  d3m0  in  the  frequency  range  from  0.3  kHz  to 
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Idle  channel  noise  at  32  kbit/s :- 

The  idle  channel  noise  at  point  B  shall  not  exceed  -60  dBmOp. 
The  level  of  any  single  frequency,  meaured  selectively,  shall 
not  exceed  -65  dBmO  in  the  frequency  range  from  0.3  kHz  to 
16  kHz. 

4 .6  Variation  of  Quantization  and  Harmonic  Distortion  with  Innut 
Level 

The  distortion  shall  be  mesured  unweighted  with  a  sinewave 
test  signal  at  800  Hz.  With  such  a  signal  applied  to  point 
A,  the  ratio  of  signal  to  distortion  power  at  the  output 
point  B  shall  be  above  the  limits  of  Figure  8. 

4 .7  Variation  of  Quantizing  and  Harmonic  Distortion  with  Frequency 

The  distortion  shall  be  measured  unweighted  with  a  sinewave 
test  signal  of  -20  dBmO.  With  such  a  test  signal  applied  to 
point  A,  the  ratio  of  signal  to  distortion  power  at  the  output 
point  B  shall  be  above  the  limits  of  Figure  9, 
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Test  Frequency  800  Hz 


f 


Fig.  9b  -  Quantizing  and  Harmonic  Distortion  with  Frequency 
at  a  Digit  Rate  of  32  kbit/s 


RATO  UNCLASSIFIED 

I  13- 


ov  oooooooo  oo  o  o  o  o  u  o  oo  o  ou 


APPErrorx  b 

CVSP  rncodj_n£  9ubrou 

su*aoutinc  ENCoociiiTruT.ouTPirr.N.Fs.FCt.rci.Fca.Tc.uMAx.uNiH.DC) 

• . CVSO  CNCCCER  SUSROUTIHE - - 

THE  S08SOUTINC  CONVERTS  A*  INPUT  TIME  FUNCTION  TO  AM  OUTPUT  BINARY 
DATA  STREAM.  JOT*  INPUT  AND  OUTPUT  is  DONE  TNSC'JCM  MJvRAYS. 

ll(lttt«l«IMIIIil  VARIABLES  ttXXTtXtXXtXttXTtttttTXXtXTttttTXtXtXX 

INPUT  •  AN  ARRAY  CONTAINING  THE  INPUT  TINE  FUNCTION  8 APPLES. 

OUTPUT  •  AN  ARRAY  CONTAINING  THE  OUTPUT  BINARY  DATA  STREAM. 

N  •  TVC  NUMBER  OF  SAMPLES. 

FCl,  fca.  FC3  *  ROU.-OFT  FREQUENCIES  OF  THE  PRIMARY  INTEGRATOR. 

TC  •  THE  TIME  CONSTANT  OF  THE  SYLLABIC  FILTER. 

ALPHA  •  T>€  DECAY  RATE  OF  THE  PRIMARY  INTEGRATOR. 

•ETA  •  THE  DECAY  RATE  OF  THE  SYLLABIC  FILTER. 

EH  *  THE  SIAM  Cr  THE  DIFFERENCE  BETUEOI  THE  CURRENT  IfFVT  AND 
THE  CURRENT  ESTIMATE. 

EM  •  THE  SIGN  OF  THE  DIFFERENCE  ONE  THE  PERIOD  AGO. 

EN3  •  THE  SIGN  OF  THE  DIFFERENCE  TWO  TIME  PERIODS  AGO. 

UIIN  •  THE  MINIMUM  INPUT  TO  THE  SYLLABIC  FILTER. 

UNAX  .  THE  MAXIMUM  INPUT  TO  THE  SYLLAtIC  FILTER. 

FS  •  THE  SAMPLE  RATE. 

DCLTAN  •  nc  CURRENT  STEP  SIZE. 

DIF  •  THE  DIFFERENCE  ISTUEEN  TW  CURRENT  INPUT  AND  THE  CURRENT 
ESTIMATE. 

m  •  THE  CURRENT  ESTIMATE. 

DC  •  THE  DUTY  CYCLE  OF  TX  SLOPE  OVERLOAD  DETECTOR  FOR  THE 
CURRENT  INPUT  STRING. 

ctntmmmimminanttsutunumimmnnnmmmtnnt 

C - SUBROUT  IX  START - 

C -  INITIALIZE  VARIABLES  AND  ARRAYS 

REAL  1  (PUTIN ) 

INTEGER  OUTPUT! N) 

DATA  XHS*/,  EUl/fy,  &8s*s.  Pl'3. 141SKCS3S/ 

SUM  •  $. 

DCLTAN  •  UMIN 

C -  CALCULATE  DECAY  RATO  OF  ENCODER  FILTERS 

ALPHA  *  EXP  <-U.  *  PI  *  FCJ  /  FJ>> 

•ETA  •  EXP  (-(3.  «  PI  /  TC  *  FSJ> 

0 -  START  ENCODING 

OO  54  I  •  t.H 

C - CALCULATE  THE  OUTPUT  OF  T*  COMPARATOR 

dif»  iNPurrd*  -  xn 
CN  «  SIGN! 1 . .DIF  J 

C - GCfCRATE  (CW  ESTIMATE 

M4  •  ALPHA  *  XN  ♦  (1  -  ALPHA)  <  DCLTAN  BEN 
U  *  UMIN 
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C -  GENERATE  THE  (CXT  OUTPUT  OF  THE  HCK  OVERLOAD  DETECTOR 

If  ( { ( (EN  .AND.  EN1 )  .AND.  EN2>  .EO.  1.)  .OR. 

1C ( (EM  .AND.  EN1>  .AND.  EN3 )  .EO.  -1.))  V  •  UHAX 

IF  cv  .eo.  vaax)  sun  ■  sun  ♦  i 

C -  GENERATE  NEXT  STEF  SI2E 

DELTAH  •  SETA  <  DELTAN  ♦  <1  -  I ETA )  t  u 

C -  SHIFT  THE  r.OFC  OVERLOAD  DETECTOR  SHIFT  REOISTEX 

EN3  •  EN1 
EN1  •  EH 
OUTPUT (I )  •  EH 

C -  PCU*  TO  SINARY  CONVERT 

IF  (EH  .EO.  -1)  OUT PUTC I)  •  • 

SR  CONTINUE 

C -  CALCULATE  SLOPE  OVERLOAD  DETECTOR  DUTY  CYCLE 

DC  ■  SUN  /  N 


RETURN 

END 
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APPENDIX  C 

cvsn  Decodin'-  Subroutine 


subroutine  Dcco«ci<iWT.eirmrr.N.r*,rci.rci.rc3.TC,un(M<.uw»,oc> 

- - - -CVSD  DECODING  SUBROUTINE - 

this  subroutine  decodes  the  binary  data  stream  contained  in  the 

INPUT  ARRAY  AND  PUTS  THE  OUTPUT  TINE  FUNCTION  SAMPLES  IN  THE  OUT¬ 
PUT  ARRAY. 

cuuimuiuuitiiuini  variables  iiuuiinininouuiuiiiiiinn 
C  IWN/r  •  AM  ARRAY  CONTAINING  THE  INPUT  BINARY  DATA  STREAM. 

0  OUTPUT  •  AN  ARRAY  CONTAINING  THE  OUTPUT  TIRE  FUNCTION 
C  N  •  THE  NUMBER  OF  SAMPLES 
C  FS  *  THE  SAMPLE  RATE 

C  FCi*  FC3.  FC3  •  ROLL-OFF  FREQUENCIES  OF  THE  PRIMARY  INTEGRATOR 
C  TC  •  THE  TIME  CONSTANT  OF  THE  SYLLABIC  FILTER. 

C  XN  •  THE  CURRENT  OUTPUT  TIME  SAMPLE 
€  ENl  ■  THE  SIGN  OF  THE  DIFFERENCE  ONE  TIME  PERIOD  AGO 

C  DC  -  THE  SIGN  OF  T>C  DIFFERENCE  TVO  TIME  PERIODS  AGO. 

C  UMAX  •  THE  MAXIMUM  INPUT  TO  THE  SYLLABIC  FILTER 

G  UNIM  •  THE  MINIMUM  IWUT  TO  THE  SYLLABIC  FILTER 

C  DELTAM  •  THE  CURRENT  STEP  SI2E 
C  ALPHA  •  THE  DECAY  RATE  OF  THE  PRIMARY  INTEGRATOR 
C  BETA  •  THE  DECAY  RATE  OF  THE  SYLLABIC  FILTER 
C  DC  •  THE  SLOPE  OVERLOAD  DETECTOR  DUTY  CYCLE. 


; - GET  NEXT  tWJT  BIT  AND  BII*AY  TO  POLAR  COTAJCRT 

CN  •  INPUT! Il 

IF  ( INPUT! I )  .CO.  •)  CN  •  -1 

C -  GENERATE  NEXT  OUTPUT  TINE  SAMPLE 

XN  •  ALPHA  *  XN  ♦  <1  -  ALPHA)  t  DEL  TAN  SEN 
V  •  WHIN 
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C - GENERATE  THE  NEXT  OUTPUT  OT  THE  SLOPE  OH*RLOAD  DETECTOR 

IF  {(«<  CN  .A*D.  EH1 )  .A>0.  EH2 )  .CO.  1.)  .OR. 

1(((EN  .AND.  EN1)  .»«D.  E'3>  .EG.  -1.))  U  •  VMX 
If  (U  .CO.  VflAX)  SUN  •  SUN  ♦  1 

C -  GENERATE  NEXT  STEP  SIZE 

DCLTAH  •  SETA  <  DCLTAH  ♦  (1  -  SETA)  S  U 

C - SHIFT  TW  SLOPE  ONLOAD  DETECTOR  SHIFT  REGISTER 

EH3  •  CN1 
CMt  >  EN 
OUTPUT (I)  •  XN 
SO  CONTINUE 

C -  CALCULATE  T>«  SLOPE  OUERLOAD  DETECTOR  DUTV  CYCLE 

DC  •  sun  /  N 

RETURN 

END 
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APPENDIX  D 
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FIR  FilterinR  Subroutine  (FIITl'R) 


SUBROUTI*  FILTERfXT,  H.  HP.  I) 

•' . . . -fin  FILTER  SUBROUTINE - 

THIS  SUBROUTINE  FILTERS  AN  INPUT  TINE  FUNCTION  SP"PLC  STRING  USINQ 
0  FILTER  COEFFICIENTS  CtNEFATED  JV  FN  EXTT'-.nrL  FILTER  CTMSATOR 
C  ROUTINE.  THE  F  ILT.T“*  £D  TII2  FUNCTION  f'PLES  FEE  FtC-CTO  IN  THE 

C  $M(  ASRAV  PS  T>-:  IN-UT  F.O  fi.;TU'.N  TO  TrZ  C.RLtI  >3  F..CC..PN.  CUE 

0  TO  TXE  FILTER  TECNNKI.E,  E.3  SABLES  AT  THE  lEOINNINO  PND  END 

0  OF  THE  SAHPLE  STSINO  A*£  LOST. 

ciitiimiuiiimuunm  variables  msimnmnmxtmnt txmxn 

C  XT  •  Tvc  PRRPV  CONTAINING  TVE  INPUT  S*»5»IE  STRINQ  AND  AFTER  PRO- 
C  CCSSINO,  THE  FILTERED  SAMPLE  STRING. 

C  N  •  THE  NUMBER  OF  SAIVLES  IN  THE  INPUT  PRRPV 

C  NR  ■  THE  NUMBER  OF  FILTER  COEFFICIENTS 

C  I  •  PM  PRRPV  CONTAINING  T>«  FILTER  COEFFICIENTS 

conttmotumtnmmQutntmnmmntmnnmtisuMtnm 

C - SUBROUTINE  START - 

C -  INITIALIZE  VARIABLES  AND  ARRAYS 

DIMENSION  XT(N).  B(W>) 


FILTER  THE  IWVT  SAMPLE  STRING 

DO  IN  I  •  1.46M 
K  •  2M  ♦  I 
m  SR  J  ■  l.NP 

IF  (J  .EO.  IJ  SUM  •  *<J)  *  XT«J 

SUM  •  SUM  ♦  B(J)  t  (XT(K  ♦  J  -  1)  ♦  XTCK  -  4  ♦  1)1 

CONTINUE 

XT<I )  •  SUM 

CONTINUE 

RETURN 

END 


V 
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APPENDIX  E 


FIR  Filter  Coefficient  Generating!  Subroutine 


subroutine  fltrocnikta.oaima.np.i) 


— MAXIMALLY  FLAT  FILTER  PROQRAM- 


TMlS  PROCRAN  OUTPUTS  THE  FIR  FILTER  COEFFICIENTS  CALCULATES  »V 
SUBROUTINE  MXFLAT.  THIS  ROUTINE  EIT-:R  riTURNS  THE  calculated 
COEFFICIENTS  TO  TVE  CALLING  fWIW  OR  P-JIMTS  CUT  THE  ERROR  MES¬ 
SAGES  UH£N  THE  COEFFICIENTS  CAPVVOT  BE  DETERMINED  DUE  TO  THE  CHOICE 
OF  IWVT  PARAMETERS. 

THIS  SUBROUTINE  AND  THE  MXFLAT  AND  PATPRX  ft!»R0UTI«S  USED  TO 
GENERATE  THE  MAXIfJLtV  FLAT  FIR  FIL~T  COTFFICJE.yTS  A'S  ADAPTED 
FROM  A  PROGRAM  MU2LCPED  IV  J.  F.  ICRIC-R  OF  £ ILL  LABORATORIES. 

THIS  PROCRAM  UAS  PL  Llf-CD  IN  *PROCRAMS  FOR  DIOITAL  SIC/iAL  PRO¬ 
CESSING*  IV  THE  IEEE  P*£SS. 


cttnmmmunnni  uariajles 

C  KTA  •  THE  NORMALIZES  CENTER  FREQUENCY  OF  THE  TRANSITION  IANS 
C  RAIMA  -  THE  NORMALIZES  UIDTH  OF  THE  TRANSITION  IANS 
C  HP  •  THE  NUMBER  OF  FILTER  COEFFICIENTS 
e  I  •  AN  ARRAY  CONTAINING  THE  FILTER  COEFFICIENTS 
C  LIMIT  -  THE  LARGEST  NUMBER  OF  FILTER  COEFFICIENTS  ALLOWED 
C  I ERR  •  THE  NUMIER  OF  THE  ERROR  MESSAGE 

C  A  I  C  •  WORKING  ARRAYS 


— SUBROUTINE  START- 


11 
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INITIALIZE  UARIAILES  AND  ARRAYS 

DIMENSION  A!299),B(299),C!299> 

LIMIT  -  20* 

CALL  NXFLAT(B£TA.  WIM,  NP.  A,  I.  C,  LIMIT,  I  ERR) 
PRINT  RESULTS 

IF  (IERR  .OT .  1)  UR  ITS  (6,  OSS! )  BETA,  OAIMA 
FORMAT! *  FOR  ICTA  •  '.FS-D.*  AND  GAMMA  •  *,FS.3) 

CO  TO  U9,  E«,  39,  49),  1EJS 
RETURN 

L»ITE(S,SSS7) 

FORMAT! •  BETA  NOT  IN  RANGE  9.  -  ,S*> 

STOP 

URITEIS.ODPS) 

FORMAT! •  CAPVM  NOT  IN  RANGE*  > 

STOP 

VRITE!S.SCSS) 

FORMAT! •  GANNA  TOO  SHALL,  MIN  IS  .94T*) 

STOP 

END 
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APPENDIX  F 

Subroutine  NXFIAT  -  Part  of  FIR  Filler  Generator 


subroutine  mxflakie.  qa.  np,  a.  i.  c.  unit,  icon 
o — . . sumoutine  pwflat - 

THIS  SUIROUTINE  COMPUTES  THE  COEFFICIENTS  OF  A  MAXIMALLY  FLAT  FIR 
LINEAR  PHASE  FILTER. 


imtmtiinmtiuniu  variables 

K  •  CENTER  OF  THE  TRANSITION  REOION,  RANGE  •  #.  TO  .S 
FREOUENCY  IS  NORMALIZED  TO  THZ  SAMPLE  RATE. 

OA  •  WIDTH  OF  THE  TRANSITION  REOION.  WHERE  THE  OUTPUT  AMPLITUDE 
DECREASES  FROM  SSX  TO  6X. 

LIMIT  -  THE  MAXIMUM  NUHIER  OF  COEFFICIENTS  IN  TIC  FILTER 
»  •  THE  ARRAY  CONTAINING  THE  FILTER  COEFFICIENTS 
IESR  •  ERROR  MESSAGES 

l,  normal  fetus* 

8.  SETA  NOT  IN  RANGE 

3,  GAMMA  HOT  IN  RANGE 

4.  GAMMA  TOO  SMALL,  LESS  THAN  .M 

C  A  >  WORKING  ARRAY 

C  C  •  WORKING  ARRAY 

C  K  •  NURtER  OF  ZEROS  AT  NVOUIST  FREQUENCY 

C  L  •  NUMSER  OF  ZERO  DERIVATIVES  AT  ZERO  FREQ 

C  NT  -  FILTER  HALF  ORDER  ■  NP  -  1 


csmsstnnswmmtmoummummtomnmainsmnnun 

C - SUIROUTINE  START - 

C -  INITIALIZE  VARIABLES  AND  ARRAYS 

DIMENSION  A(LIMIT),  MUMIT),  CUINIT) 

I ERR  •  1 

MP  ■  • 

TWOPI  •  ».  *  ATAN(l.G) 

IF  <<IE  .IE.  «.)  .OR.  (IE  .CE.  .SI)  QO  TO  SS 
■N  -  AM INI (2.  I  BE.  1.  -  3.  *  IE) 

IF  ( (CA  .LE.  «.)  .Clt.  CCA  .CE.  SMI)  90  TO  M 
KT  •  INK  1.  /  (4.  1  CA  I  CA)) 

IF  (NT  .OT.  IE8)  CO  TO  Ul 

ac  •  (i.  ♦  cesrruopi  *  IE))  /  i. 

GLIM  •  LIMIT 

CALL  RATPRXIAC,  NT,  K.  NP.  QUHI 

n  «  a  *  hp  -  i 

IF  (K  .EG.  •)  K  •  1 
©—  COMPUTE  MAGNITUDE  AT  NP  POINTS 
CU)  >  1. 
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Subroutine  RATPRX  -  Part  of  FIR  Filter  Generator 
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SUM  •  SUH  ♦  DELTA 

H  COMTIME 

C -  CALCULATE  DUTY  CYCLE  AMS  AUERACE  STEP  SIZE 

C 

CO<I-l>  •  1.  /  I 

STEPU-1)  •  SUM  ✓  SAA. 

1M  COMTIME 

C -  CALCULATE  T>C  AUERACE  STEP  WAUE  FOR  A  DUTY  CYCLE  Of  2ERO. 

C 

DELTA  >  A. 

SUM  >  A. 
u  •  umn 
DO  20A  J  •  l.SAA 

DELTA  •  EETA  *  DELTA  ♦  (1.  -  SETA)  I  U 
SUM  •  SUM  ♦  DELTA 

MA  COMTIME 

CD(SA)  •  A. 

STEP(SA)  •  SUM  /  SAA. 

0 -  PRINT  AMS  PLOT  T>«  RESULTS. 


C 


l  AAA 


PRINT  *.  •  FOR  THE  FOUOUINO  SYSTEM  PARAMETERS! • 

PRINT  t.*  SA WU  RATE  •  \FS.‘  IPS* 

PRINT  t„*  TC  •  *.TC 
PRINT  *,  •  FC1  ■  *,FC1 

PRINT  «,  ■  RATIO  •  ‘.^TIO 

PRINT  *.*  THE  FOLLCUINO  SYSTEM  PARAMETERS  ARE I* 

PRINT  «,*  UMAX  •  ■.L'.VX,*  UftlN  •  *,WUM 

PRINT  *,*  MAXIMUM  STEP  SIZE  •  *,ST£P(1).*  MIMIMUM  STEP  SIZE  •  *, 
lSTEP(SA) 

IF  (IT  .CT.  1)  CO  TO  SAA 
CALL  F ACTOR (.S) 

CALL  PLOT ( 2 . ,  2..  -3) 

CALL  SCALEICO.lO.S.it.l) 

CALL  SCAL£(STEP.S.d,r^.n 

CALL  AXIS(R.e,e.A.i:-i.’FP  SIZE  (U>.i3,6.A.9*.A.STEP(Cl>,STEP(53>) 
CALL  AXIStO.A.A.A.ltnDVTV  CYCLE.-li.iA.A.A.A.CDlSl >,00(63  J > 

CAU.  F2CT(A.,  A.,  It.,  A.,  3) 

CONTINUE 
I CHAR  -  IT  -  1 

CALL  LIN£(CD,STEP»SA,1,1A,ICHAR1 

CONTINUE 

CALL  PLOTE(N) 

end 
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VMAX  Calculation  Subroutine 


SUIROUTIX  MAXOPTIMAX.MIH.FS.FCi.TC.RATIO) 


0  THU  SUBROUTINE  CALCULATES  the  CALL'S  Of  UNAX  END  MIN  BA-5ED  ON 
*'  THE  IN°UT  S*.*5IE  RAT",  SVLLAOIC  SILTS'!  TC.  FM’MV  1HTTC  ATOR 

ROLL-OFF  FSSC'.LhCV  (SCI).  I  J  TX  .nTIO  BiTL'TV.H  THE  P.'Xir.Ufl  STEP 
SIZE  AND  him  .1  S*sa  513?  C<J7?VT  C-  T:  E  tv.'-  3IC  SIL*"A.  ■’HIS 
CALCULATION  !S  FZ  (•*C'.-:8  AT  A  !.  ..-T-'tNOi  FRSO'  ..(CY  CF  '  »  HU  AMO 
SIGNAL  smiTctS  CS  0  n:.T8.  TX  U.ILU-S  CS  L.  J<  AM)  U.'IN  ARE 
CALCULATED  SUCH  THAT  THE  DUTY  CYCLE  OS  THE  SLCHE  OVERLOAD  DETECTOR 
OUTPUT  IS  .5. 

iittinittisisiiimnmi  uariailes  mmmmmuonaiaimnii 
MAX  «  the  AAXIRUA  IXUT  TO  THE  SYLLAIIC  FILTER 
MIN  •  TVS  AIHIRUR  It*»UT  TO  THE  SYLLAIIC  FILTER 
FS  •  The  SAMPLE  RATE 

FC1  •  THE  ROLL-OFF  FREGUDCY  OF  THE  PSIRARY  INTEGRATOR 

TC  •  THE  TIRE  CONSTANT  OF  THE  SYLLAIIC  FILTER 

RATIO  •  THE  RATIO  rETVEEN  THE  RAXJRUR  STEP  SIZE  AND  THE  RINIRUR 
STEP  SIZE  OUTPUT  Of  THE  SYLLAIIC  FILTER 

TS  •  AN  ARRAY  CONTAINING  THE  TEST  SIGNAL  SARPLES 

•IHOUT  •  AN  ARRAY  CONTAININQ  THE  BINARY  OUTPUT  OF  TX  ENCODER 

PEAK!  -  THE  PEAK  UALUE  OF  THE  TEST  SICNAL  ARPLITUDE 

DIF  •  THE  DIFFERENCE  SSTUEEN  THE  SLOPE  OUERLOAD  DETECTOR  DUTY 
CYCLE  AND  THE  DESIRED  VALUE  OF  .5. 

RAT  •  TX  RATIO  OF  TX  DUTY  CYCLE  DIFFERENCE  TO  TX  DESIRED  UALUE 

mmtmnmmnt  susroutiws  used  nmmnmammmnuni 


SIGNAL  •  TX  TEST  SICNAL  GENERATOR 

URINOPT  •  CALCULATES  THE  VALUE  OF  VAIN  THAT  PAIRS  UITH  TX  CAL¬ 
CULATED  VALUE  OF  VRAX 

C  ENCODE  I  •  TX  CVSD  ENCODER 


0 - 5U3ROUTIX  START - 

0— —  INITIALIZE  ARRAYS  AND  VARIABLES 

DIMENSION  TSUOOS) 

INTECER  IIKCUTUCSS) 

AIDER!)  •  SLHTUS.  XX ( (DINS  -A. VIS. )  S  .SSI  S  SSS. >  S  SORTIS. ) 
PEAK!  •  ACS. ) 

UMAX  •  IS. 

0—  CALCULATE  TEST  SIGNAL  SAMPLES 


C— 


CALL  S1GNAL(T$,4S9S.FS,SSS..S..PEAK1#S. ) 

—  START  MAX  CALCULATION  LOOP 
CONTIHJC 

—  CALCULATE  ESTIRATED  ENCODER  PARARETERS 
CALL  URINOPT (MAX, UR1N, Ft, TC. RATIO) 

—  PROCESS  TX  TEST  SIUttL 

CALL  CNC0D£l(TS,IIN0UT,4A9S,F$,FCt,FCS,FC3,TC,MAX,MXM.DC ) 

—  FIND  TX  DIFFERENCE  KTVEEN  TX  OUTY  CYCLE  USINQ  TX  ESTIRATED 
MAX  AND  VHIN  A.  J  THE  USIkED  DUTY  CYCLE 

DIF  •  DC  -  .5 
RAT  •  DIF  /  .S 


( 
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ou 


IF  TVC  DUTY  CYCLE  IS  WITHIN  1*  OF  THE  DESIRED  UAUJE  RETURN 
THE  UNAX  AND  WHIN  VALUES  TO  THE  CALL  I  HQ  PHCQRAN 

IF  (ASS (RAT )  .LE.  .*1)  00  TO  9H 

OTHERWISE  REE ST 1  NATE  VNAX  AM)  REPEAT  CALCULATIONS 

WAX  •  UMAX  ♦  .5  »  RAT  *  UNAX 
CO  TO  5 
9N  CONTINUE 
RETURN 
END 
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on  o  o  o o  no  o  o  o  o  t»  rtt» r. 


APPENDIX  J 


V'*IN  Calculatin'?  SiihrouLi  n<"- 


SUBRCUTIX  LNINOPT(VnAX.umN,FS.TC. RATIO) 
'imttuitiiitmttmimitmnnxmumnttmtnutmnumm 


THIS  SUIROUTIX  CALCULATES  THE  UALLt  OP  WIN  THAT  PAIRS  UITH 
THE  VALUE  OF  WAX  T**T  13  U  XfT  SO  THAT  THE  RATIO  CF  TX  fMX- 
INUN  STEP  SIZE  TO  NINIM.il  STEP  SITE  AT  T,  OUTPUT  OF  THE  SVLLABIC 
FILTER  IS  UITrtlN  .61*  Of  THE  VALUE  SPECIFIED. 

uuuittuiuimmuti  variables  (immmmmmmmnmmi 

WAX  •  THE  NAXIflUn  IMVT  VALUE  OF  THE  SVLLABIC  FILTER 

WIN  •  THE  HININUN  VALUE  INPUT  TO  IX  SYLLABIC  FILTER 

FS  •  TX  SAMPLE  RATE 

TC  •  THE  TINE  CONSTANT  OF  TX  SVLLA*IC  FILTER 

BAXSTEP  •  TX  mXINUM  STEP  SIZE  AT  TX  OUTPUT  OF  TX  SYLLABIC  FIL¬ 
TER 

MNSTEP  •  TX  NININUM  STEP  SIZE  AT  TX  OUTPUT  OF  TX  SYLLABIC  FIL¬ 
TER 

BETA  •  TX  DECAY  RATE  OF  TX  SYLLABIC  FILTER 
SUB  •  TX  RUNNING  SUN  OF  STEP  SIZES 

RATIO  -  TX  DESIRED  RATIO  SETVZEN  TX  FAX  I  HUN  STEP  SIZE  AND  TX 

MNinun  step  size  at  the  output  of  the  svuajic  filter  in  db 

C  R  -  TX  VOLTAGE  RATIO  EQUIVALENT  OF  RATIO 
C  KIT  A  -  TX  CURRENT  STEP  SIZE 

cwnmomonnpnntmumn»int»oxiaroomww»«mi 

C -  INITIALIZE  VARIABLES  AND  ARRAYS 

REAL  HAXSTEP.  HINSTEP 
DATA  PI/3. HIS v  :S3«/' 

R  •  It.  TX  (RATIO  /  £6. 1 

€ -  CALCULATE  SYLLABIC  FILTER  DECAY  RATE 

•ETA  •  EXP  (-(3.  »  PI  /  TC  /  FBI) 

C -  ESTINATE  INITIAL  VALUE  OF  WIN 

UNIH  •  UMAX  /  IBB. 

C  START  CALCULATION  LOOP 

S  CONTINUE 

C -  INITIAL  RUNNING  SUM  AMD  STEP  SIZE 


BUN  •  B. 

DELTA  •  B. 


<*-— 


IB 


CALCULATE  AUERAGC  MAXIMUM  ITEP  SIZE 


DO  IB  I  •  t.SBB 
IDT  •  |  3 

IT  •  IDT  (  3 
V  •  WIN 

IF  lit  .EO.  I)  V  •  UMAX 
DELTA  •  IETA  «  DELTA  ♦  (1 
SUN  •  SUN  ♦  DELTA 
CONTINUE 

NAXSTEP  •  SUN  /  SBB. 


BETA)  I  V 


c 


REINITIALIZE  RIPflINO  BUR 
KLTA  •  B. 

BUR  •  B. 


AND  BTEP  SIZE 
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1 


c 


CALCULATE  CURRENT  ESTIMATE  OF  THE  RINIHUN  STEP  SIZE 


IS 

c- 


c - 

c 

909 


do  is  i  •  i.sj* 

DELTA  *  BETA  *  DELTA  ♦  <1.  -  BETA)  X  IA1IH 

Sllfl  •  SUM  ♦  DELTA 

CONTINUE 

ninstep  •  sun  /  so*. 

FIND  THE  DIFFERENCE  BETWEEN  THE  ESTIMATE  AND  T»«  SPECIFIED  RATIO 

THIN  •  NAXSTEP  /  » 

DIF  •  TMN  -  « INSTEP 

RAT  •  ASS  (DIF)  /  THIN  X  1M. 

IF  THE  DIFFERENCE  IS  LESS  THAN  .*1S  RETURN  UNIN 

IF  (RAT  .LE.  .*1)  CO  TO  999 

IF  THE  DIFFERENCE  IS  GREATER,  THEN  REESTIHATE  W1IN  AND  REPEAT 
CALCULATIONS 

UNIN  •  UNIN  ♦  .5  X  DIF 

CO  TO  S 

CONTINUE 

RETURN 

END 


l 


i 

! 


I 
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APPENDIX  K 


CVSP  li'T.t.o-i  c'. i,< •  ’ 


PROGRAM  PULSE  ( INPUT, OUTPUT, TAPE 0* OUTPUT, PLOT ) 

C  THIS  PROCBAM  PLOTS  T«E  OUTPUT  SICNAL  OP  THE  CUSD  TRANSMISSION 
C  SYSTEM  WHEN  TnE  I“«PUT  SIGNAL  IS  ft  003  HERTZ  SINE  UAvE 

C  TVftT  UARIES  1 5  *  -  =  L.  I  TO . :  £  FPfH  ;?"3  TO  0  I  ‘  ;5.  thE  TEST 

C  SIGNAL  CE^EPATCR  AITE^na^ELY  cEMbaTES  SCO  SABLES  AT  -ftp  DON* 

C  AND  0  DC“0  IN  ORDER  TO  IT  -NSTRATe  T-E  SYSTEM'S  ST  TP  rESPONSE 

C  CHARACTERISTICS.  TV-  "  SY-j'Z.N  UN  Oil  TEST  CONSISTS  CP  TV£  INPUT 

C  FILTER, THE  CUSD  ENCCEER  AND  DECOIER.  AND  ThE  OUTPUT  FILTER. 

ctttnottmituttniMi  uariailes  immmiiiiiixmumimium 

C  TSIN  •  AN  ARRAY  CONTAINING  THE  INPUT  TINE  SERIES  SAMPLES. 

C  TSOUT  •  AN  ARRAY  CONTAINING  FIRST  THE  DECODER  OUTPUT  TI"E  SERIES 
C  SAMPLES,  THEN  THE  OUTPUT  TINE  SERIES  SAMPLES  OF  THE  FIR  FILTER. 

C  »  •  AN  ARRAY  CONTAINING  THE  FILTER  COEFFICIENTS. 

C  TINE  •  AN  ARRAY  CONTAINING  THE  TIPS  THAT  THE  FIRST  29#  SANPLES 
C  ARE  TAKEN  SO  THAT  THEY  NAY  BE  PLOTTED. 

C  BINOUT  -  AN  ARRAY  CONTAINING  THE  BINARY  OUTPUT  OF  TIC  CUSD  ENCODER 

C  AMP1  •  THE  AMPLITUDE  OF  THE  TEST  SIGNAL  IN  DBM. 

C  FS  •  THE  SAMPLE  RATE. 

C  FC1,  FC2.  FC3  •  THE  ROLL-OFF  FREQUENCIES  OF  THE  PRIMARY  INTEGRA- 

C  TORS. 

C  TC  •  THE  TIME  CONSTANT  OF  THE  SYLLABIC  FILTERS. 

C  UMAX  t  UNIN  •  THE  MAXIMUM  AND  MINIMUM  INPUTS  TO  THE  SYLLABIC  FIL- 

C  TER. 

C  BETA  -  THE  NORMALIZED  CENTER  OF  THE  TRANSITION  BAND  OF  THE  LOU 
C  PASS  FILTER. 

C  GAMMA  •  THE  NORMALIZED  UIDTH  OF  THE  ROLL-OFF  REGION  OF  THE  OUTPUT 
C  FILTER.  THE  REGION  IS  THE  FREOUENCY  BAND  BETWEEN  THE  B5«  AND 

C  S*  OUTPUT  AMPLITUDES. 

C  PEAK1  •  THE  MAXIMUM  AMPLITUDE  OF  THE  TEST  SIGNAL  IN  UOLTS. 

C  HP  •  THE  NUMBER  OF  FILTER  COEFFICIENTS. 

C  DC  •  THE  DUTY  CYCLE  OF  THE  SLOPE  OVERLOAD  DETECTOR. 

c  ratio  •  the  ratio  between  the  maximum  step  size  and  t»c  minihun 
c  step  size  in  db 

CiiMuuuuMUuii  subroutines  used  immiiiiiuumiamimim 

C  FLTRGEN  •  THE  SUBROUTINE  THAT  GENERATES  THE  OUTPUT  FILTER  COEFFI- 
C  CIEHTS. 

C  FACTOR,  PLOT,  AXIS,  SCALE.  RECT,  LINE,  PLOTE  ■  CALCOMP  PLOTTING 
C  ROUTINES 

C  SIGNAL  •  THE  TEST  SIGNAL  GENERATOR.  PRODUCES  SAMPLES  OF  SINUSOI- 
C  DAL  UAUCS  WITH  AT  MOST  TWO  FREOUENCY  COMPONENTS. 

C  ENCODE 1  •  THE  CUSD  ENCODER  SUBROUTINE  WITH  A  SINGLE  ROLL-OFF 
C  FREOUENCY  IN  THE  PRIMARY  INTEGRATOR. 

C  DECODE1  •  THE  CUSD  DECODING  SUBROUTINE  WITH  A  SINGLE  ROLL-OFF 
C  FREQUENCY  IN  THE  PRIMARY  INTEGRATOR. 

C  FILTER  •  THE  SUBROUTINE  that  FILTERS  THE  INPUT  TIME  SERIES  SANPLES 
C  USING  THE  FILTER  COEFFICIENTS  GENERATED  BY  FLTRGEN. 

C  PLTIHE  •  A  LINEAR  PLOTTING  ROUTINE  TO  PLOT  SIGNAL  AMPLITUDE  US. 

C  TIME. 

C  UMAX OPT  •  GENERATES  UMAX  AMD  VAIN  FOR  THE  CUSD  ENCODER  AND  DECODER 
C  SUBTUTINCS 

cimmtmmmniutuntttttmmmmimnittmttaimimtit 

Ctuuiiiiimm:tiii<imiinnmmuiiuitiiimuiu:atiu>i»iu 
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PROGRAM  START 


C -  INITIALIZE  UPRIABLES  AND  ARRAY9 

DIMENSION  TSINlS0W),TSOUT(S###>,TIX(2SM).B(2M) 

INTEGER  »INCUT(Sce») 

A<DBfl#>  •  SGRTU*.  It <  (DIN#  -4.  I/'t#. )  t  .Ml  <  CM.  )  S  CORTC3.) 
AMPi  •  -48. 

AMP2  •  •. 

PEAK!  •  A  < ANP1 ) 

PEAKS  •  A  (AMP2) 

KN  •  •. 

C -  INPUT  AND  PRINT  TX  WORKING  UARIABLt* 

READ  »,  PS 

READ  F.FCl.  TC.  PATIO 

READ  t.BETA,  GANNA 

PRINT  *,  •  TRANSIENT  RESPONSE  TEST  AT  • ,F9, *  BPS* 

PRINT  *.  ■  WITH  TC  ■  '.TC,'  AND  RATIO  •  *, RATIO 

PRINT  FILTER  PARAMETERS  ARE.  BETA  *  ‘.BETA,*,  GANNA  .  *,OMWA 

C -  GENERATE  FILTER  COEFFICIENTS  AND  CUSD  SYSTEM  PARAMETERS 

CALL  FLTRGEN(BETA.GA«TA,NP,B) 

CALL  UNAXOPT (UMAX, UNIN.FS.FC1.TC, RATIO) 

C -  INITIALZE  PLOTTER 

CALL  FACTOR (.5) 

CALL  PL0T<2.,  8..  -3) 

C -  GENERATE  INPUT  TIME  FUNCTION  SAMPLES 

CALL  SICNAL2(TSIN.S»M.FS.MS..S..PEAK1.PCAK8) 

C -  FILTER  THE  INPUT 

CALL  FILTER (TSlN.SSSS.NP.t) 

C -  PROCESS  THE  INPUT  TIME  SERIES  THROUGH  TX  CVSD  SYSTEM 

CALL  ENCODE1 (TSIN, BINOUT. 5088, FS.FCi.FCP, FC3.TC, UMAX, UMIN, DC  > 

CALL  DECODE1  (3INOUT, TSOUT, 5238,FS,FC1,FC2,FC3,TC, W'flAX.UfllH, DC ) 

C -  FILTER  THE  OUTPUT  OF  THE  DECODER 

CALL  FILTER<TS0UT.46M.NP.») 

C - PLOT  TX  OUTPUT  SIGNAL 


uv  a  1  - 

TIME! I)  •  1  /  FS 
5  CONTINUE 

DO  G  I  •  1.22S 
IK  •  2858  ♦  I 
TSOUT(I)  •  TSOUT(IK) 

C  CONTINUE 

CALL  PLTIFE(TINE,TSOUT ,225.887) 

CALL  PLOTEtN) 

END 

cnttitttmnmmmiiimmitmimiiintmimummnnmnt 

Citmtuittuiiuxtiiuiitiiutiiitouitiiiunummunuaitium 

SUBROUTIX  PLTIMEIX.  V,  N.  NX) 


C . - . -TIME  US.  AMPLITUDE  PLOTTER - 

C  THIS  SUBROUTIX  MAKES  A  LIXAR  PLOT  OF  TIME  US.  AMPLITUDE. 
Cdiiiiiiititttttmutnii  uariailes  lutttttunuiinnnmutmn 

C  X  •  TX  ARRAY  CONTAINING  TX  ORDINATE  UALUE9 

C  V  •  TX  ARRAY  CONTAINING  TX  ABSCISSA  UALUES 

C  M  •  TX  NUMBER  OF  UALUES  IM  TX  X  AND  V  ARRAYS 

C  NX  •  N  ♦  3 

cnmmmmmtmi  subroutixs  used  minimmimnuunitun 

C  SCALE,  AXIS.RECT,  PLOT.  PLOTS,  LIX  ■  CALCOMP  PLOTTING  POUTINES 

cntimtuiiitumiiiiitiitininmaiuittiiiutiniinuunmiui 


11? 


c - Sm ROUTINE  STRUT - 

C - INITIALIZE  UARIABLES  AND  ARRAYS 

DIMENSION  XINX).  V(NX) 

C -  SCALE  THE  X  AND  V  ARRAYS 

CALL  SCALE (X,  It..  N,  1) 

CALL  SCALE ( Y.  6..  N.  1) 

C -  DOX  IN  TX  PLOT 

CALL  RECTI*..  1*..  *..  3) 

C - DRAW  THE  AXES 

CALL  AXIS!*..  *..  1BHTIME  (SEC).  -I*,  t*..  *..  XlfHI),  X(N*8>) 

CALL  AXIS!*.,  13HAMPLITUCE  (U>,  13.  6..  9*..  VIH+l).  V(N+8>) 

C -  PLOT  THE  POINTS 

CALL  IINEIX.  V.  N.  1.  *.  •; 

RETURN 

END 

Cmmtiimmmiimmtnimtmmimmmntmnmnumnt 

Ctnmimmiutiimmummtxmutmuuuiiiiiitmmumiii 


SUBROUTINE  $ICMAL3<OUTPUT,N,FS,FREQ1,FRE02, AMP1 , AHP3 ) 

CUiKKUitmiutmimiuiniuiMiitiuiuiimmiumittiiuQii 

C  THIS  SUBROUTINE  GENERATES  A  SINGLE  FREGUENCY  SINUSOIDAL  SIGNAL 
c  that  alternately  *-<>s  see  samples  at  one  amplitude  then  s»*  samples 

C  AT  A  SECOND  AMPLITUDE. 

cittmimmmiimi  uariables  nuummmmiiitummunmi 
C  OUTPUT  •  THE  ARRAY  CONTAINING  THE  OUTPUT  TINE  FUNCTION  SAMPLES 
C  N  •  THE  NUMBER  Of  SAMPLES  TO  K  PRODUCE# 

C  FS  •  THE  SAMPLE  RATE 

C  FREOl  •  THE  FREGUENCY  OF  THE  TEST  SIGNAL 
C  FRE02  •  UNUSED 

C  ANP1  •  THE  AMPLITUDE  OF  TX  FIRST  SM  SNPLES 

C  AHP2  •  TX  AMPLITUDE  OF  TX  SECOND  54*  SAMPLES 

C  AMP  •  TX  SIGNAL  AMPLITUDE  CURRENTLY  BEING  USED 

ctunnmmitnntintttiimiiitiHiinmmnunmtBmttnini 

c - SUBROUTINE  START - 

C -  INITIALIZE  UARIABLES  AND  ARRAY# 

DIMENSION  OUTPUT(N) 

DATA  PI/3.141S92SS3B/ 

K  •  • 

S  CONTINUE 

C -  SET  CURRENT  SIGNAL  AXLITUDE 


AMP  .  ANPI 
AMP l  •  ANP8 
AMP3  •  AMP 
J  •  • 

1*  CONTINUE 


C -  GENERATE  SB#  TIX  FUNCTION  SAMPLES 

J  •  J  ♦  I 

IF  (J  .OT.  S**>  GO  TO  S 
K  •  E  ♦  I 


C - If  M  SAMPLES  HAUC  BEEN  OCXRATED.  STOP  PROGRAM 


999 


IF  (K  .GT.  N)  GO  TO  999 

OUTPUT(K)  •  AMP  <  SIN  (8. 

CO  TO  IB 

CONTINUE 

RETURN 

END 


B  PI  t  FRCOi  '  FB  <  K) 
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appendix  l 


Idle  Channel  rro.'nn 


PROGRAM  IOLENOI<IIFVr,ajTPUT.TAPE**OUTPU?.PLOT> 

: - IDl£  CHAUCL  NOISE  program - 

S  THIS  PROGRAM  N ENSURES  THE  IDLE  OSFWCL  NOISE  OF  THE  CUSD  TRPNS- 
C  MISSION  SVSTEN  U  IN  IKE  LNCC1  'R  AND  t£CCM »  A“E  CCWCCTED  IACX- 
C  TO-IaOC  AND  THE  INPUT  TO  T>^£  ENCODER  IS  GROUNDED. 

C  THE  SVSTEN  SAIN  IS  ADJUSTED  SO  THAT  AN  PM  KZ  INPUT  SI5NAL  AT 

C  -28  DBN8  PRODUCES  A  -£3  D U*3  SIGNAL  AT  TrZ  OUTPUT  OF  THE  DE- 

0  CODER. 

C 

ctuutuoimtmnin]  uariailes  naiammmnaooaoongn 

C  FRE01  •  THE  REFERENCE  FREQUENCY  USES  TO  SET  TIC  SYSTEM  GAIN. 

C  TSIN  •  AN  ARRAY  CONTAINING  THE  INPUT  TINE  FUNCTION  SAMPLES. 

C  TSOUT  •  AN  ARRAY  CONTAINING  FIRST  TIS  PS CODER  OUTPUT  TIME  FUNCTION 
C  SAMPLES,  THEN  THE  OUTPUT  TIKE  FUNCTION  SAMPLES  OF  THE  FIR  FIL- 

C  TER 

C  *  ■  AN  ARRAY  CONTAINING  TIC  FILTER  COEFFICIENTS. 

C  SINOUT  •  AN  ARRAY  CONTAINING  THE  IINARY  OUTPUT  Of  THE  CUSS  ENCODER 
C  AMP1  •  THE  AMPLITUDE  Of  THE  REFERENCE  SIGNAL  IN  DM18. 

C  FS  •  THE  SAMPLE  RATE. 

C  FC1,  FC2,  FC3  •  THE  ROLL-OFF  FREQUENCIES  OF  THE  PRIMARY  INTECRA- 

C  TORS. 

C  TC  •  THE  TIME  CONSTANT  Of  THE  SYLLABIC  FILTERS. 

C  UMAX  L  WHIN  •  THE  MAXIMUM  AMS  MINIMUM  INPUTS  TO  THE  SYUASIC  FIL- 

«  TER. 

C  SETA  •  THE  NORMALIZED  3  El  FREQUENCY  OF  THE  OUTPUT  FILTER.  TIC 
«  FREQUENCY  IS  NC.VAU2ED  TO  THE  SAMPLE  RATE. 

C  GANNA  •  THE  NORMALIZED  WIDTH  OF  TIT  ROLL-OFF  REGION  OF  TIC  OUTPUT 
C  FILTER.  1>*  f  •  DIC.T  IS  THE  FREQUENCY  (AND  (ETWEEN  THE  S6X  AND 

C  EX  OUTPUT  AMPLITUDES, 

C  PEAKS  •  THE  MAXIMUM  AMPLITUDE  OF  THE  TEST  SIQNAL  IN  W0LT9. 

C  NP  •  THE  NUBER  OF  FILTER  COEFFICIENTS. 

C  DC  •  TIC  DUTY  CYCLE  OF  THE  SLOPE  OVERLOAD  DETECTOR. 

C  PIN  •  TIC  POWER  OF  TIC  IfCUT  SIGNAL  IN  DOW. 

C  POUT  •  THE  POWER  OF  TIC  OUTPUT  SIGNAL  IN  DIMS. 

C  ICN  •  TIC  CALCULATED  IDLE  CHANNEL  NOISE  IN  DSN8. 

?  GAIN  •  TIC  WCLTACC  AMPLIFICATION  OF  TIC  SYSTEM. 


SUBROUTINES  USED 


FLTRCEN  •  THE  SUBROUTINE  THAT  GENERATES  TIC  OUTPUT  FILTER  COEFFI¬ 
CIENTS. 


C  SIGNAL  •  THE  TEST  SIGNAL  GOCRATCA 
C  ENCODE 1  •  THE  CUSD  ENCODER 

C  DECODE t  •  THE  CUSD  DECODER 

C  FILTER  •  THE  SUDROUTIIC  THAT  FILTERS  THE  IICUT  TIME  FIHCTION  SAN- 
C  PLES  USIHO  TIC  FILTER  COEFFICIENTS  GENERATED  1Y  FLTRCEN. 

C  POWER  •  A  ROUTIfC  TO  CALCULATE  THE  POWER  IN  A  SAMPLED  TIME  FUC- 
C  TION  WITH  IMPEDENCE  •  68*  OKfW. 


cmmininnn»nrMtnwmmitttnrt:mnm«umunMU»p 

cmumimmiuimiuaiwmwmNm^m«uninammuj 


lid 


r»e* 


-(TOtW  START- 


N 

C— 


—  INITIALIZE  VARIABLES  AND  ARRAYS 

DIMENSION  TSINCS999).  TS0UT(S9M),  B(BM) 

REAL  ICN 

INTECER  BINOUTISMB) 

A(DW»)  •  SGRTU9.  tJUDSN*  — VI#. )  t  .Ml  t  CM.)  *  9QRT(3. ) 

—  INPUT  AND  PRINT  WORKING  VARIABLES 

READ  t.  FRE01.  ANP1.  FS 
READ  S,FC1,  TC.  r^TIO 
READ  S.BETA,  CA-7A 

PRINT  *,  •  IDLE  CHStfCL  NOISE  TEST  AT  *,FS.'  IPS' 

PRINT  S,  •  WITH  TC  •  •  .TC,  ■  '*  T  ■  aTIO  •  '.PSTIO 
PRINT  *.  •  OUTPUT  FILTER  PARAMETERS  ARE  I  SETA  •  *.»£TA 
PRINT  »,  •  GAiVM  •  •,6A71A 

—  GOCRATE  T>C  FILTER  COEFFICIENTS  AND  CUSD  SYSTEM  PARAMETERS 

CALL  FLTRC5EN(»ETA,C<WA.^,») 

CALL  VnAXCPTIUtAX.WIIN.FS.FCl.TC. RATIO) 

—  GENERATE  INPUT  TINE  FUNCTION  SAMPLES 
PEAK!  •  A  (ATPi) 

CALL  SICNAL<TSIN,S*M.FS.FREQ1.S..P£AK1,C.) 

—  process  nc  urur  tii*  function  tvroucm  the  cvso  system 

CAU  ENCOOE1  (TAIN.  BIXN  JT,  5^*9,  FS.FC1  .FCS.FC3,  TC.VIAX.IJNIN.  DC) 
CALL  DCC0DeU»IN0UT.TS0UT.E6M,FS.FCl.Fca.FC3,TC,UWW,UNIN.DC) 

—  FILTER  THE  OUTPUT  OF  THE  DECODER 
-•ALL  FILTERiTSOUT.SMC.NP.I) 

—  DELAY  THE  INPUT  SIGNAL  START  TO  CORRESPOND  TO  TW 
FILTERRED  OUTPUT. 

DO  39  ID  •  1.4496 

kd  •  aca  ♦  id 
TSIN(ID)  •  TSIN(O) 

CONTINUE 

—  CALCULATE  THE  REFERENCE  SYSTEM  MIN 

CALL  P0L€R<TSIN,40SS.FS.PIN> 

CALL  POl!£R(TSOUT.4;:3,FS,POUT> 

Cain  •  sort  ipin/pout) 

—  GENERATE  A  ZERO  INPUT  SIGNAL  ARRAY 

DO  45  I  •  I, SMS 
TSIN(I)  «■  «. 

CONTINUE 

—  PROCESS  THE  ZERO  SIGNAL  THRCUQM  T*C  SYSTEM 

CALL  ENCOCEIITSIN.BINCUT .59A9.FS.FC1  .FCa.FC3.TC, UNAX, UNIN, DC) 
CALL  OCCOC€l(SINOUT,TSOUT,S<M.FS.FCI,FC<I.FC3,TC,LnAX.WUN.OC) 

—  FILTER  T)C  OUTPUT  SIGNAL 
CAU  FILTERITSOUT.SMS.MM) 

—  ADJUST  THE  OUTPUT  SIGNAL  AMPLITUDE  TO  THE  REFERENCE  VALUE 

DO  46  I  -  1,4996 

TSOUTfl)  •  TSOUT(I)  S  GAIN 

CONTINUE 

—  CALCULATE  THE  IDLE  CHAAtCL  NOISE 

CAU  POUER(TSOUT,4996.FS,PCUT) 

ICN  •  19.  S  ALCG19(PGUT ) 

PRINT  OUT  TIC  RESULTS 

URITEC6.6M)  ICN 

FORNATdX, ’THE  IDLE  OM9CL  NOISE  •  '.Fi.C) 

END 
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00  000  **  00  0000  00  r>  non  000 


APPENDIX  H 


Total  Harnon  ic  Distortion  Fro  "ran 


'*>5**'  MMDISTaVMUT,  OUTPUT,  TAPES.  IVMUT.TAPU»OUTPUT> 

. . total  harmonic  distortion  pro^am - 


THIS  PROGRAM  CALCULATES  TV*  TOTAL  HARMONIC  DT*nCRTtCN  TH  TVM  OlfTPiiT 
-  UMEN  a  SlNCLE  fi-r.  :ncv  TT-t  S!-*AL  IS  r^cclf  'rg rl™  “Hi  IcuJn" 
c  ENCODER  and  CECCC..;»  C'S*  TOTED > --'<--0-:  J,X.  *  D-V’T’-a  fa  ^rtBAT-n 
c  us  Ik;  only  wc -c  ‘  - E c t  - l  cc  pi'-chts  of  tvs  clTi/’ut  tAvt  l*s  1 

t  160  HZ  *V<D  4t‘ v.)  H2.  Tr.C.^  AT  tXACT^V  JM  H2  Cr3  4026  k2  HOT 

C  INCLUDED  IN  If-  CALCULATION.  **  *“  h0T 

cituiiumittimmi  uariailes  uuarammtmnummmmm 


c 

c 

e 

e 


c 

c 

c 

c 

c 

c 


‘“encoder.  1KTE0CT  ****  C0WT*INI*°  TW  »1NA«V  OUTPUT  Of  T>C  CUSD 

OUTPUT  •  A  .C^L,.*p!,*v  CONTAINING  The  TIPS  FUNCTION  OUTPUT  OF  TVS 
OUTPUTIOF:'THi£t£CoIifl.W,D  F,0CE5SIN!3<  THE  THE  FUNCTION 

**  *  C^tM^T^OUTPL'T  SPECTRAL  PC SR  CONPOWNTS 

TWMWORPi  iSSoSn^!  ™  ^SSIW1  w  ™£  £AST  «***« 

IUK,  IRC.  CUE  «  UORKIHO  ARRAYS  USED  »V  THE  FFT  SUWOUTIVC. 

FREOUE  •  A  REAL  ARRAY  CONTAINING  TH*  F~r,L£NCIES  AT  UH1CH  nr  cjrr 
MAS  CALCULATED  THE  SPECTRAL  CCkVCNOTTS.  m 

FREQ  •  TVC  FREQUENCY  OF  THE  TEST  SISNAL  IN  HZ. 

AW  •  THE  AMPLITUDE  OF  THE  TEST  SIGNAL  IN  DIMS. 

FS  •  THE  SAMPLE  RATE  IN  IPS. 

A  •  THE  PEAK  UALUE  OF  THE  TEST  SIGNAL. 

rcv&  si»  ocow? f<*  **  ”inc:pi£  iff7tawToR 

TC  •  THE  COMPANDING  SPEED  OF  THE  SYLLAIIC  INTECRAfTOR  IN  SEC. 

T  •  A  FRE0UENC0  UARIASLE  USED  TO  tETERMIPE  TVS  TEST  StOHW.  gpro- 

cofyorcfirc^^T  ALS0  T0  setewiivc  the  harmonic  shectral 

*•  FREQUENCY  POUO>  W  T*«  OUTPUT  SPECTRAL  COMPONENT  AT  THE  TEST 

REF  •  THE  OUTPUT  SPECTRAL  COMPONENT  POUE*  IN  DIM. 

**•  *  THE  RUHNINO  SIM  OF  TVC  POWER  CF  THE  HARMONIC  CflWCVCMTS. 

THD  •  THE  TOTAL  HARMONIC  DISTORTION  IN  «. 

•ETA  •  NORMALIZED  3  Dt  FREQUENCY  OF  TVC  OUTPUT  FILTER 

•SANA  ■  TVC  NORMALIZED  ROLL-OFF  BANDWIDTH  OF  TVC  OUTPUT  FILTER 

RATIO  •  TVC  RATIO  OF  TVS  MAXIMUM  ST?*  SI7F  TO  THE  MINIUM  RITA 
SIZE  IN  THE  CUSD  ENCODER  AND  CSCticER.  QIUENOM  D». 


•ttttu suircutines  used  nmntmmtnn 

ENCODE  1  •  TVC  CUSD  ENCODER 
DECODE  1  •  TVC  CUSD  DECODER 

W*£0£Jb^«'C*ATES  UMAX  AND  UNIN  USED  IN  THE  CUSD  ENCODER  PND 

FLTERQEN  •  TVC  COCFFICENT  OEVCRATOR  FOR  TVC  OUTPUT  FILTER 
FILTER  •  FILTERS  THE  OUTPUT  SIGNAL. 

fTrU»R(Mr^  F**7  r0W£*  T**Mro*"  SUMOUnvC  FROM  TVC  ZNSL 


Cliiuittttt>mititrrtrti!tTrrtRvnrrnmmnrT?rmmimmin 

CUitititiiiimuimiiuiitiijiiiiiiiimiiimn  mmiiCTOitoinm 
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DO  tS  J  ■  l.KF 

IF  (  F  .re.  FR£3UC(J>)  60  TO  0 
sun  •  sun  ♦  fsx(j) 
ts  continue 

30  CONTINUE 

0 -  CALCULATE  TOTAL  HARMONIC  DISTORTION 

THD  •  SORT  (SITU  /  E*  >  IS*. 

URXTE<e,£«6)  THS 

•M  FORMAT* IX, 'THE  TOTAL  HARMONIC  DISTORTION  It  ‘.Ft.!,**.*) 
END 


I 


\ 

I 


r 


uo  ooo  uu  ooo  u  uu  o  o  u  o  u  o  .. i.  'j  ,j  o  «  o  uu  o  o  uuo  uo 


APPENDIX  N 


Total  Harmonic  "‘intortion  vs,  Tn*^t  co  ml  Po'-.'^r 

PROGRAM  DTHDt INPUT, OUTPUT, TAP£S-1NPUT.TAPE6*0UTPUT,PL0T> 

C - THD  US.  INPUT  POWER - 

C  THIS  PROGRAM  INVESTIGATES  THE  VARIATION  IN  HARMONIC  DISTORTION 
C  AS  THE  SIGNAL  INPUT  PCM*?  IS  VARIED.  THE  INPUT  PCUER  IS  C^tCED 

C  IN  .4  D3  STEFS  FROM  -4e  D3T3  TO  3  DM.  TH;  HARMONIC  DISTORTION 

C  IS  THE  OUTPUT  IS  THEN  TEAS  WED  AND  PLOTTED. 

C  THE  PROGRAM  IS  REPEATED  THREE  TirES.  STEPPING  THE  STEP  SIZE 

C  RATIO  FROM  32  DO  TO  3  3  08.  THE  TWEE  SETS  OF  DATA  ARE  THEN 

C  PLOTTED  ON  THE  SA,  i  CRAPH. 

ottmmtmmmm  variables  mmtRunitmuuinmmmnn 


INPUT  -  AN  INTEGER  ARRAV  CONTAINING  THE  BINARY  OUTPUT  OF  THE  CVSD 
ENCODER. 

OUTPUT  •  A  REAL  ARRAV  CONTAINING  THE  TIME  FUNCTION  OUTPUT  OF  THE 
TEST  SIGNAL  GENERATOR  AND  AFTER  PROCESSING,  THE  TIRE  FUNCTION 
OUTPUT  OF  THE  DECODER. 

POWER  •  A  REAL  ARRAY  CONTAINING  THE  POWER  THAT  EACH  SAMPLE  IS 
TAKEN. 

PSX  ■  A  REAL  ARRAY  CONTAINING  THE  OUTPUT  SPECTRAL  POWER  COMPONENTS 
OF  THE  DECODER  OUTPUT  AFTER  PROCESSING  BY  THE  FAST  FOURIER 
TRANSFORM  SUBROUTINE. 

I UK,  UK,  CMC  •  WORKING  ARRAYS  USED  BY  THE  FFT  SUBROUTINE. 

FREQUC  -  A  REAL  ARRAV  CONTAINING  TV?  FREQUENCIES  AT  WHICH  THE  FFT 
HAS  CALCULATED  THE  SPECTRAL  COMPONENTS. 

N  •  THE  NUMBER  OF  TIME  SAMPLES  TO  BE  TAKEN. 

FREQ  •  THE  FREQUENCY  OF  THE  TEST  SIGNAL  IN  HZ. 

AMP  ■  THE  AMPLITUDE  OF  THE  TEST  SIGNAL  IN  DBMQ. 

FS  •  THE  SAMPLE  RATE  IN  BPS. 

A  •  THE  PEAK  VALUE  OF  THE  TEST  SIGNAL. 

UA.  UK,  CUK  •  WORKING  ARRAYS  USED  BY  THE  FFT  SUBROUTINE. 

FREQUC  •  A  REAL  ARRAY  CONTAINING  THE  FREQUENCIES  AT  UHICH  THE  FFT 
MAS  CALCULATED  THE  SPECTRAL  COMPONENTS. 

II  •  THE  NUMBER  OF  TIME  SAMPLES  TO  BE  TAKEN. 

FREQ  •  THE  FREQUENCY  OF  THE  TEST  SIGNAL  IN  HZ. 

AMP  .  lMt  AMPLITUDE  OF  THE  TEST  SIGNAL  IN  DBMO. 

FS  •  THE  SAMPLE  RATE  IN  BPS. 

A  •  TX  PEAK  VALUE  OF  THE  TEST  SIGNAL. 

FC1.  FC2,  FC3  •  ROLL-OFF  FREQUENCIES  FOR  THE  PRINCIPLE  INTEGRATOR 
IN  THE  CVSD  ENCODER  AND  DECODER. 

TC  •  THE  COMPANDING  SPEED  OF  THE  SYLLABIC  INTEGRATOR  IN  SEC. 

I,  J,  K  •  COUNTING  INDICES  FOR  THE  VARIOUS  *DO*  LOOPS. 

F  •  A  FREQUENCY  VARIABLE  USED  TO  DETERMINE  THE  TEST  SIGNAL  SPEC¬ 
TRAL  COMPONENT  AND  ALSO  TO  DETERMINE  TX  HARMONIC  SPECTRAL 
COMPONENTS. 

£•  •  THE  RMS  POWER  OF  THE  OUTPUT  SPECTRAL  COMPOfCNT  AT  THE  TEST 
FREOUEMCY. 

C  REF  •  THE  OUTPUT  SPECTRAL  COMPONENT  POWER  IN  DBM. 


119 


oo  o  o  oo 


C  sun  •  THE  RUNNING  SUM  OF  THE  POWER  OF  THE  HARMONIC  COMPONENTS. 

TKD  «  AN  ARRAY  CONTAINING  THE  VALUE  Of  HMfWHIC  DISTORTION  AT  EACH 
LEVEL  Of  INPUT  POWER. 

I  *  AN  ARRAY  CONTAINING  THE  OUTPUT  FILTER  COEFFICIENTS. 

HP  •  THE  NUMBER  OF  OUTPUT  FILTER  COEFFICIENTS. 

BETA  •  THE  NORMALIZED  CENTER  OF  THE  TRANSITION  BAND  FOR  THE  OUTPUT 
FILTER. 

C  GANNA  -  THE  NORMALIZED  UDTH  OF  THE  OUTPUT  FILTER  TRANSITION  BAND. 


c- 


■PROCRAM  START' 


C -  INITIALIZE  VARIABLES  AND  ARRAYS 


DIMENSION  INPUT( ) , OUTPUTfSCA* ) , POWER (202 > . PSX(  ISO ) 
l«IUK(20),l-K<r>a>,Ffv£CLE<15*>,  THD(2v>2>.  I(22a) 


COMPLEX  CUX(BJa) 
mv DBMO )  •  SCAT (10. 


»  ( ( DBMO  -  4.)  /  ID.)  *  .Ml  *  600. )  *  SORT 


1 1 

1CHAR  •  -1 


INPUT  WORKING  VARIABLES 

PRINT**, ^DYNAMIC  RANGE  TEST  AT  ‘.PS,’  BPS  AND  *,FREOi,*  HZ* 
READ  *.  FC1 .  TC 

READ  *,  XLEN,  YUEN,  XHIH.  XMAX,  YMIN,  YMAX 

XSTEP  •  (XT'.  AX  -  XniN)  /  XLF.N 

VSTEP  •  (VflAX  -  YNIN )  /  VLEN 

PRINT  *,*  TC  •  *,TC,\  FC1  •  *,FC1 

READ  *.BETA,  GAW1A 

PRINT  *,*  EETA  •  '.SETA.',  GANNA  •  ‘.CAMMA 
CALL  FLTRGENtBETA, GANNA, NP.B ) 


START  LOOP 

DO  1000  NR  -  2,6,2 
RATIO  *30.  ♦  NR 

CALL  VMAXOPTf UMAX, VHIN,FS,FC1,TC, RATIO) 

DO  SOO  IS  •  1,100 
POUER(IS)  •  -AO.  ♦  .A  *  IS 

DETERMINE  PEAK  VALUE  OF  TEST  SIGNAL 

AflPl  •  A(POUER(IS)) 

GENERATE  INPUT  TINE  FUNCTION 

CALL  SICNAL(OUTPUT,5MO»FS.FRE01.O.,AMPI,#. ) 

PROCESS  THE  TINE  FUNCTION  THROUGH  THE  CVSD  SYSTEM 

CALL  ENCODE 1 (OUTPUT, IwPUT,SC80,FS,FCl,FC2,FC3,TC, VNAX,VMIN,DC ) 
CALL  DCC0CC1 ( INPUT, OUTPUT, 5500, FS.FCl , FC2, FC3, TC, UMAX, VNJ N,DC ) 


C -  FILTER  TW  OUTPUT 

CALL  FILTER(OUTPUT,SOOO,NP,B) 

C -  DETERMINE  THE  SPECTRAL  COMPONENTS  OF  THE  OUTPUT 

CALL  F7FPS(0UTPUT,DUM,A«9S,2SS,0,PSX,DUM,DUM,IUK,UK,«X,IER) 

c -  DETERMINE  THE  COMPONENT  AT  THE  TEST  FREQUENCY  AND 

C  CALCULATE  REFERENCE  VALUES. 
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DO  3  K  •  3,123,2 
KF  •  Kf  ♦  1 
f  •  »<K-1  )/2S6. )  t  FS 
FREOUE(KF)  •  f 

IF  <PSX(K>  .  LE.  S.E-19)  PSXOC)  -  S.E-10 
PSXIICF)  ■  PSXOO 
IF  (F  .NE.  FREQ1 1  GO  TO  8 
EO  •  SQRT<  PSXOC)  ) 

•  CONTINUE 

C -  CALCULATE  POWER  AT  HARflONIC  FREQUENCIES 

SUH  -  0. 

DO  20  I  •  2,10 
F  •  FREOl  /  I 

IF  <<F  .LE.  180.)  .OR.  <F  .GT.  4000.))  CO  TO  20 
DO  IS  J  •  l.KF 

IF  (F  .NE.  FREOUE(J))  GO  TO  15 

sun  •  sun  ♦  psxu) 

15  CONTINUE 

20  CONTINUE 

DO  30  1  >  2,30 
F  •  FREOl  *  I 

IF  ((F  .LE.  1C0.)  .OR.  (F  .GC.  4000.))  GO  TO  30 
DO  25  J  •  1.ICF 

IF  (  F  .NE.  FREQUE(J))  GO  TO  2S 
SUN  •  SUN  ♦  PSX(J) 

25  CONTINUE 

30  CONTINUE 

C -  CALCULATE  TOTAL  HARflONIC  DISTORTION 

THD(IS)  •  SORT  (SUfl)  /  E3  *  100. 

IF  <THD(1S)  .GT.  100.)  THD(IS)  •  100. 

500  CONTINUE 

C -  PLOT  RESULTS 

IC8A&  •  ICMAP  ♦  1 

CALL  PLTRANG  ( POWER,  THD.100, 1  CHAR,  XHIN.XLEN,  XSTEP.  YNIN,  VLEN,  YSTEP  ) 
1000  CONTINUE 

CALL  PLOTE(N) 

END 


SUBROUTINE  PLTRANOtX, V,N,  ICHAR.XH IN, XLEN, XSTEP, VWN. VLEN, YSTEP ) 
C - DVNANIC  RANGE  PLOT  SUSROUTI* - 


C  THIS  SUBROUTINE  CREATES  A  PLOT  OF  THE  V  ARRAV  VERSUS  THE  X  ARRAY. 

ctutttitittttttttttttttx  variables 

C  .<  •  AH  ARRAV  CONTAINING  THE  ORDINATE  UALUES 

C  V  •  AN  ARRAV  CONTAINING  THE  UALUES  TO  BE  PLOTTED. 

C  H  •  THE  NUMBER  OF  VALUES  IN  THE  ARRAYS 

cmtmimmminumimmmmiimtmmmitmuinmtnn 

ctmmiiminmmmnnmnnmjnimjmmmnminmmu 

C - SUBROUTINE  START - 

C -  INITIALIZE  ARRAYS 

DIMENSION  X( 4096 ),  V(4096) 

X(N+1 )  •  XNIN 

X(N+2)  •  XSTEP 

V<NM)  •  VniN 

V(N42)  •  VSTEP 

IF  (ICHAR  .QT.  0)  GO  TO  500 

C -  ESTABLISH  NEW  PAGE  ORIGIN 

CALL  FACTOR (.5) 

CALL  PL0T<2..  2.,  -3) 
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C -  BOX  IN  THE  GRAPH 

CALL  RECT<0. ,  0..  YLEN,  XLEN,  0.,  3) 

C -  DRAW  THE  AXES 

CALL  AXIS<0. 0,0.0, 1SHTNPUT  POVER  (CDf10),-18,XLEH,e  0.XNIN.XSTEP) 
CALL  AXIS<e.0,0.0,14HDISTCRTION  (*>,H,YL£N, W.0.YNIN.  YSTEP) 

C -  PLOT  VALUES 

600  CONTINUE 

CALL  LItiE(X,Y,N,l,l$,lCHAR) 

RETURN 

END 
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"  i  r: r;  ti'li<?ri  Vo!  :i]  l|:r”on ;  ■'  Di  r. ‘orl:  i on  vr  .  In10-:1-  P-i  'n~il  Po v;n r 


PROCRAfl  MMTHD( INPUT, OUTPUT, TAPES-INPUT,TAPE6«0UTPUT, PLOT) 

C - MISMATCHED  THO  US.  INPUT  POWER - 

C  THIS  PROGRAM  INVESTIGATES  THE  VARIATION  IN  HARMONIC  DISTORTION 
C  AS  THE  SICNAL  I>rUT  PCUER  IS  UARIED.  THE  I'PUT  POWER  IS  CHANGED 

C  IN  .4  DB  STEPS  FROM  -41  D0“3  TO  8  C:T3.  THE  HARMONIC  DISTORTION 

C  IS  THE  OUTPUT  IS  THEN  MEASURED  AND  PLOTTED. 

C  THE  PROGRAM  IS  REPEATED  THREE  TIMES.  WHILE  THE  ENCODER  PARAMETERS 

C  ARE  HELD  CONSTANT.  THE  DECODES  STEP  SICE  RATIO  IS  ALLOWED  TO  CATV 

C  FROM  32  D3  TO  33  D3.  THE  THREE  SETS  OF  DATA  A St  THEN  PLOTTED  ON 
C  THE  SANE  GRAPH. 

cummmitmumt  uariables  mmimmtmtimmmmtmm 

C  INPUT  -  AN  INTEGER  ARRAY  CONTAINING  THE  BINARY  OUTPUT  OF  THE  CUSD 
C  ENCODER. 

C  OUTPUT  •  A  REAL  ARRAY  CONTAINING  THE  TIME  FUNCTION  OUTPUT  CF  THE 
C  TEST  SIGNAL  GEf-SRATCR  AND  AFTER  PROCESSING,  THE  TIFE  FUNCTION 

C  OUTPUT  OF  THE  DECODER. 

C  POUER  •  A  REAL  ARRAY  CONTAINING  THE  POUER  THAT  EACH  SAMPLE  IS 
C  TAKEN. 

C  PSX  *  A  REAL  ARRAY  CONTAINING  THE  OUTPUT  SPECTRAL  PCUER  COMPONENTS 
C  OF  THE  DECODES  OUTPUT  AFTER  PROCESSING  BY  THE  FAST  FOURIER 

C  TRANSFORM  SUBROUTINE. 

C  I  UK,  UK.  CUK  •  WORKING  ARRAYS  USED  BY  THE  FFT  SUBROUTINE. 

C  FREQUE  •  A  REAL  ARRAY  CONTAINING  THE  F^EGUEMCIES  AT  WHICH  THE  FFT 

C  HAS  CALCULATED  THE  SPECTRAL  COMPONENTS. 

C  N  •  THE  NUMBER  OF  TIKE  SAMPLES  TO  BE  TAKEN. 

C  FREO  •  THE  FREQUENCY  OF  THE  TEST  SIGNAL  IN  HZ. 

C  AMP  •  THE  AMPLITUDE  OF  THE  TEST  SIGNAL  IN  DBN*. 

C  FS  -  THE  SAMPLE  RATE  IN  BPS. 

C  A  •  THE  PEAK  UALUE  OF  TF€  TEST  SIGNAL. 

C  FC1,  FC2,  FC3  •  ROLL-Cf F  FRCOUENCIES  FOR  THE  PRINCIPLE  INTEGRATOR 
C  IN  THE  CUSD  ENCODER  AND  DECODER. 

C  TC  •  THE  COMPANDING  SPEED  OF  THE  SYLLABIC  INTEGRATOR  IN  SEC. 

C  I,  J,  K  •  COUNTING  INDICES  FOR  THE  UARIOUS  'DO*  LOOPS. 

C  F  •  A  FREQUENCY  UARIABLE  USED  TO  DETERMINE  THE  TEST  SIGNAL  SPEC- 

C  TRAL  COMPONENT  AND  ALSO  TO  DETERMINE  THE  HARMONIC  SPECTRAL 

C  COMPONENTS . 

C  £•  •  THE  RMS  POUER  OF  THE  OUTPUT  SPECTRAL  COMPONENT  AT  THE  TEST 
C  FREQUENCY. 

C  REF  •  THE  OUTPUT  SPECTRAL  COMPONENT  POUER  IN  DBM. 

C  SUM  •  THE  RUNNING  SUM  OF  THE  POUER  OF  THE  HARMONIC  COMPONENTS. 

C  THD  •  AN  ARRAY  CONTAINING  THE  UALUE  OF  HARMONIC  DISTORTION  AT  EACH 
C  LEVEL  OF  INPUT  POWER . 


cmmmtnmntttmmnmmmmttmmntmnnnmmmn 

cniuutttmitmumiumitimmtuiimimimmmitmiuiu 
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■PROGRAM  START 


C 

C -  INITIALIZE  VARIABLES  A NO  ARRAYS 

DIMENSION  INPUT (EC’D ) .OUTPUT (5309  )  ,  POL'ER(E02),PSX(  150  ) 
1,U;!C<M),UK(1S0>.EREGUE(1S0>,  THDS£F(£92),  TWD(202>.  BOM) 
CONPLEX  CUKC3C9) 

A(DBN0)  •  SCfiTtie.  XX  ( (DBH0  -  4.)  /  IB.  )  X  .Ml  X  EM.  )  X  SORT (2. 
1) 

I  CHAR  •  -1 


C -  INPUT  WORKING  VARIABLES 

READ  X.  EREOl ,  ES 

PRINT  X.*  DYNAMIC  RANGE  TEST  AT  '.ES.1  BPS  AND  ’.FREOi, *  HZ* 
READ  X,  EC1.  TC 

READ  X.  XL£N,  VLEN,  XNIN,  X’AX,  YHIN,  YMAX 

XSTEP  •  (X  *»X  -  XNIN)  /  X'.EN 

ySTEP  •  (YMAX  -  VP IN )  /  VLEN 

PPIHT  t, *  TC  •  ‘.TC.*,  EC1  •  *,EC1 

PEhD  L.BETA.  GANNA 

PRINT  t,*  BETA  •  ‘..ETA,*,  GANNA  •  *, GANNA 
ChIL  FLTRGEN(BETA,CA«V1A.NP,B> 


C - SThRT  loop 

DO  10M  NR  •  2,6.2 
RATIO  •  33.  +  NA 

CALL  VNAXOPTI VNAX.VNIN.ES, ECl.TC, RATIO) 
IE  ( I  CHAR  .GE.  8)  GO  TO  2 
EVHX  •  UMAX 
EVHN  •  VNIN 
2  CONTINUE 

DC  500  IS  •  1,100 

POUER(IS)  •  -40.  ♦  .4  X  IS 

C -  DETE RHINE  PEAK  VALVE  OF  TEST  SIGNAL 


ANPi  «  ALPOUERUSn 

C -  GENERATE  INPUT  TINE  FUNCTION 

CALL  SIGNAL<OUTPUT,5Oe0,FS, EREQ1.0. .ANPI,®. ) 

C -  PROCESS  THE  TINE  EUNCTION  THROUGH  THE  CVSD  SYSTEM 


CALL  ENCODE! (OUTPUT, INPUT, 5O30.ES, EC1,FC2,EC3,TC, EVNX.EVNN, DC ) 
CALL  DECODEK  INPUT, OUTPUT,  5000, F5,FC1.FC2,FC3,TC,  UMAX,  VMIN.DC) 


C -  FILTER  THE  OUTPUT 

CALL  FILTER  ( OUTPUT.  SCM.NP.B) 

C -  DETERMINE  THE  SPECTRAL  COMPONENTS  OF  THE  OUTPUT 

CALL  FTFPS (OUTPUT, DUN, 4096, 256, 0.PSX, DUN, DUM, IUK,UK,CUK,IER) 

C -  ELIMINATE  THE  COWENTS  AT  ODD  MULTIPLES  OF  THE  SAMPLE  RATE. 

C  DETERMINE  THE  COMPONENT  AT  THE  TEST  FREQUENCY  AND 

C  CALCULATE  REFERENCE  VALVES. 

KE  •  • 

DO  8  K  ■  3,129.2 

KF  •  KE  ♦  1 

F  •  ( (K-l  )/256. )  X  ES 

FREGUE(KF)  •  E  „ 

IF  (PSX(K)  .LE.  6.E-10)  PSX(K>  •  6.E-1® 

PSX(KE)  •  PSX(K) 

IF  (F  .NE.  EFE01 )  GO  TO  8 
E®  •  SORT (  PSX(K)  > 

•  CONTINUE 

C -  CALCULATE  POUCH  AT  HARMONIC  FREQUENCIES 
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IS 

29 


25 

» 


sun  •  «. 

00  20  I  ■  2,1$ 

F  •  FFJE01  /  I 

00  1SF  t‘bEi  ir0*’  •°*'  lF  *0T*  oo  TO  29 

If Jf  ;NE.V&3UE(J»  CO  TO  IS 
85n\$"  *  PSX,J) 

CONTINUE 
00  33  I  •  2.30 
F  -  FRE01  X  I 

00  ’  '°fi'  if  ,Qe*  4^®e*»  CO  TO  39 

IF l  F  ,NS.'fk30UE(J))  CO  TO  2S 
SU1  •  SUN  ♦  PSX< J ) 

CONTINUE 

CONTINUE 


C— —  CALCULATE  TOTAL  HARMONIC  DISTORTION 


tHD(IS)  ■  SORT  (SUN)  /  E9  *  109. 

_  IF  (THD(IS)  .CT.  I09.)  THD(IS)  •  199. 
590  CONTINUE 


C - PLOT  RESULTS 


ICHAR  •  ICHAR  ♦  i 

CALL  PLTRAnGIPOUER.THO,  109. ICHAR.xniN.XLEN.XSTEP.YniN.YLEN.YSTEP) 

Ifw  wflTXnub 

CALL  PLOTE(N) 

END 


CXtXXXXXtXXXXlXXXIXtXXItXXXXTXXXXXTXXTXTXtXTXXXXXXXXXXXXIXXltXSXTXXX  XX* 

6txxxxxxxxxxxxxxxxxxxxxxxxzxzzxxxxxxxxxxxxxxzxxzxzs.xxizxxzxxxtxizixx7  ux 

SUBROUTINE  PLTRANQtX, Y,N,  ICHAR,X?1IN,XL£N,XSTEP,  VMIN,  YLEN,  YSTEP ) 

c - - - dynamic  range  plot  subroutine - 

C  THIS  SUBROUTINE  CREATES  A  SEMI-LOG  PLOT  OF  THE  X  AND  Y  ARRAYS. 

C txtxtxttxtxxxxxxxxxxtxxx  UARIA8LES  xxxxtxxxt xxxxxxxzxxxxzxxxxxxxxxxxxn 
C  X  •  AN  ARRAY  CONTAINING  THE  ORDINATE  UALUES 

C  Y  •  AN  ARRAY  CONTAINING  THE  UALUES  TO  BE  PLOTTED. 

C  N  •  THE  NUMBER  OF  UAUES  IN  THE  ARRAYS 


cnxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxtxxxtnxxxxxxxxxxtxx 

cminutmmu>x»mimi»mnminmmuimummxxmmx 

C - ■ - SUBROUTINE  START - - 

C -  INITIALIZE  ARRAYS 

DIMENSION  X(4096  >,  Y<499€> 

X(H+1)  •  XMIN 

XCN+2)  •  XSTEP 

Y<NM)  •  VN!N 

Y(N*2 )  «  YSTEP 

IF  (ICHAR  .GT.  9)  GO  TO  S99 

C -  ESTABLISH  NEU  PAGE  ORIGIN 

CALI  FACTOR (.5) 

CALI  PL0T(2.,  2.,  -3) 

C -  BOX  IN  THE  GRAPH 

CALL  RECT(0.,  9.,  YLEN,  XLEN,  9.,  3) 

C -  DRAU  THE  AXES 

CALL  AXIS(9.9,9.9, 18HINPUT  PCUER  <DBN9).-18,XLEN,9.9,XNIN(XSTEP) 
CALL  AXIS(0.9,0.0.  HHOISTCRTICN  <*>,  14,  YLEN,  90.0,  VftIN,  YSTEP  > 
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PLOT  UAUJES 


CALL 

RETURN 

END 


APPENDIX  P 


In  to  modulation  rinlort  ion  I’m;." 


MOCMfl  l NTERNDt  INPUT,  OUTPUT,  TAPES*  If^uT.TAPC**  OUTPUT  > 

C . --INTERNODULATION  DISTORTION  PROOR/* - 

C  THIS  PRO  GRAN  CALCULATES  TVS  iNTFC-OCULATtON  DISTO»Tir«t  OP  A 
c  cuss  system  Li  z  txi  £•«.'•:«  P-.3  z  :cc:*.*  F-1  CCrvECTtD  BACK-TO- 
C  SACK.  DISTORTION  IS  r~t>:  rj  SV  I'PUTTtrO  A  TT.GT  SI*  ,*1  CO  "POSED 

c  op  tuo  cgual  (■  -litl  *  .airs  at  i : .  .1  *z  c-j  7*i  ki.  t>  - 

C  AMPLITUDE  CP  T-"  F."  --JCT  IS  Ty  N  t  "ASL  -;g  *,ND  COMPARED 

C  TO  TX£  INPUT  SIGNAL  TO  UTZSMIni  THE  PERCENT  D1STC/.T1CM. 

cmuiimunilttltu  VARIABLES  mimnmmnmttmtmiminn 


C  IWVT  •  AN  INTERS  ARRAY  CONTAINING  T»€  UNARY  OUTPUT  OF  THE  CU6D 
C  ENCODER . 

C  OUTPUT  ■  A  F.TAL  A7«AV  CONTAINING  THE  TIPS  FUNCTION  OUTPUT  OP  THE 
C  TEST  SIG'UiL  CLM7ATCG  AND  **TER  PROCESSING.  THE  TINE  FUNCTION 

C  OUTPUT  OF  THE  CECCCZR. 

C  PSX  •  A  REAL  ATAAY  CONTAINING  Tw£  CLNTUT  SPECTRAL  POLCR  COMPCTCNTS 
C  OF  THE  IZCCIER  OUT-I/T  ( FTTR  PROCESSING  BY  THE  FAST  FOURIER 

C  TRANSFORM  < FFT  J  SLl^CUTINE. 

C  I  UK,  UK,  CUK  «  UORKINO  ARRAYS  USED  »V  THE  FFT  SUBROUTINE. 

C  FREQUE  -  A  REAL  ARRAY  CONTAINING  Tvc  FCE'-JENCIES  AT  UHICH  THE  FFT 

C  HAS  CALCULATED  THE  SPECTRAL  CCLPCNENTS. 

C  N  •  TFC  NUNBER  OF  TINE  SAMPLES  TO  BE  TAKEN. 

C  AflPt  «  THE  AMPLITUDES  Of  TPC  TEST  SIGNALS  IN  DIM. 

C  rs  •  THE  SAMPLE  RATE  IN  BPS. 

C  PEAK  •  T>C  PEAK  UALUES  Of  THE  TEST  SIGNAL  COMPOFCNTS. 

C  A  *D*  OR  'E*  P7EFIX  Oh  THE  NEXT  TK’EE  SETS  OF  UARIA5LE  INDICATES 
C  TVC  VARIABLE  IS  USED  BY  EITHER  THE  DECODER  CH  ENCODER,  RESPEO- 

C  TN*LY 

C  FC1.  FC8.  FC3  •  P*U-C5T  FPfQLSNCIES  FOP  THE  PRINCIPLE  INTEGRATOR 
C  IN  THE  CVSD  ENCODER  AND  DECODER. 

C  TC  •  THE  COMPAND  INO  SPEED  OF  T»C  SVUABIC  INTEGRATOR  IN  SEC. 

C  RATIO  •  THE  TAX  I  NUN  STEP  SIZE  TO  RININUN  STEP  SIZE  RATIO  IN  Ct 

C  r  •  A  FREQUENCY  VARIABLE  USED  TO  TETTERIM  THE  TEST  SIGNAL  SPEC- 

C  TEAL  CO  -CfCNT  AND  ALSO  TO  DETLaAINE  THE  FWRNONIC  SPECTRAL 

C  COMPONENTS. 

C  REF  •  THE  OUTPUT  SPECTRAL  COMPONENT  POUER  IN  DBA. 


'  SUM  •  T>C  SUN  OF  T>«  POWER  IN  THE  TEST  SIGNAL  COMPONENT*. 

•*  NP  •  THE  NUMBER  OF  FILTER  COEFFICIENTS 

0  B  •  AN  ARRAY  CONTAINING  T>e  INPUT  AND  OUTPUT  FILTER  COEFFICIENTS 

iiuiKiittimmuut  subroutines  USED  luatniannnuuuuunn 

C  ENCODE  1  •  THE  CUSD  ENCODER 

C  DECODE 1  •  THE  CUSD  DECODER 

C  VR1AXOPT  •  DETERNtfCS  THE  VALUES  Of  UMAX  AND  UNIN  USED  IN  T»C 
C  ENCODER  AND  DECODER 

C  fLTRGEN  •  GENERATOR  Of  TW  OUTPUT  FILTER  COEFFICIENTS 

C  FILTER  •  FILTERS  THE  OUTPUT  USING  THE  COEFFICIENTS  PRODUCED  *Y 
C  FLTRGEN 

C  SIGNAL  •  THE  TEST  SIGNAL  GENERATOR 

C  FTFPS  •  THE  FAST  FOURIER  TRANSFORM  FROM  THE  IMSL  LIBRARY 

cmwimiiiiumtaiinmmunmmctamHnmmmmmo 
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c - -  'Program  start - - - 

C -  INITIALIZE  VARIABLES  AND  ARRAYS 

DIMENSION  INPUT(SP3«).OLTPUT(Se«4>.PSX(lS#> 

1, IUK(28  >,UK(lS6>.F;iiOUE(294>.  1(244) 

REAL  IMiCD 
COMPLEX  CUK ( 304 ) 

AlCSMO)  •  SORT  ( 14.  (t  ( (DBA4  •  4.V14.)  t  .441  I  644.)  *  SQRT(2. > 

c - imrr  and  print  ucwxinq  variables 

READ  *,  A.MPI,  FS 
READ  *,  EFC1 ,  ETC.  E^ATIO 
READ  t.  f?Cl,  DTC.  CCAT10 
READ  *.  SETA,  r  .'"A 

PRINT  *,•  INTEND  T^ST  FOR  INPUT  SIGNALS  OP  7F4  AND  1444  HZ  AT  \ 
1AMP1,*  D6A0,  AM)  S.UPLE  KATE  •  *,FS.*  0/3* 

PRINT  *,*  ETC  •  ‘.ETC,  *,  E.-Cl  •  •.E.-Cl.*,  EPATIO  •  '.ERATIO 
PRINT  *,*  DTC  •  • ,  DTC ,  • .  DEC!  •  *,  DFCi,  * ,  C-iATIO  •  ‘.JRATIO 
PRINT  *,•  SETA  ■  ".SETA,*,  CArtflA  •  * .GAMMA 

©—  DETERMINE  PEAR  VALUE  OP  TEST  SIGNAL 

PEAK  .  A  (ANPl) 

C -  GENERATE  INPUT  TINE  FUNCTION 

CALL  SIGNAHOl/TPUT ,  5444,  FS,  ?S4..  1444..  PEAK,  PEAK) 

C -  GENERATE  FILTER  COEFFICIENTS  AND  CVSD  SYSTEM  PARAMETERS 

CALL  UNAXOPT(Fm*X,rvW,FS,F»ei,ETC.ET><»TIO) 

CALL  VNAXOPTf  CVTX,Cv:'N,FS.DPCl, DTC, DRAT 10 ) 
vriLL  FLTRCENi 6tTA,GAnMA,NA,B ) 

' - PROCESS  THE  TINE  FUNCTION  TVROUGH  THE  CUSD  SYSTEM 

Call  ENC0DE1  (OUTPUT,  INPUT, F304.F9.EFC1.FC3.FC3, ETC. EVfK.EVMN, DC) 
CALL  DCCODEK  INPUT. OUTPUT, St44,FS,DfCl,FCa,FC3.DTC,DVR<.DVW1.DC) 

C - FILTER  THE  OUTPUT 

CALL  FILTER (OUTPUT, S444.NP.I) 

c -  subtract  the  average  value  from  the  sample  string 

sum  •  4. 

S>0  6  I  -  1,4006 
SUM  «  SUM  ♦  OUTPUT  (I ) 

6  CONTINUE 

AUER  •  SUM  /  4996 

DO  7  1  •  l,4ecs 

OUTPUT!  I )  •  OUTPUT!  I  )  -  AVER 

7  CONTINUE 

C -  CALCULATE  THE  SPECTRAL  COMPONENTS  OF  THE  OUTPUT 

CALL  FTFPS  ( OUTPUT,  DUM,  4496, 256, 4,  PSX.DUfl.DUM,  XUT.UC.CUC,  IER) 

C  DETERMINE  THE  COMPONENT  AT  THE  TEST  FREQUENCY  <WD  CALCULATE 
C  THE  (STS  VOLTAGE. 

SUM  •  4. 

DO  S  K  •  2,129 
F  •  ( (K-l  )/cS6.  )  S  FS 
FREOUE(K-l)  •  F 
PSXOC-l)  •  P5X(K ) 

IF  < (F  .NE.  714.)  .and.  (F  .N£,  1444.))  00  TO  S 
SUM  •  SUM  ♦  SCRT  (  PSX(K)  ) 

REF  •  14.  S  AL0014(PSX(K )  /  644.  /.til) 

•  CONTINUE 

C -  CALCULATE  INTERMODULATION  DISTORTION 

DIF  •  2S4. 

DO  16  I  •  I.12S 

IF  (F REOCJE(I)  OIF)  GO  TO  16 

DIS  •  PSX(I) 

»  CONTINJE 

I  NOD  •  SORT  (DIS)  /  SUFI  t  144. 

IF  (IflOO  .OT.  104.)  IHOD  •  144. 

WRITE (6,646 )  IPOD 

444  FORMAT! IX, ‘THE  INTERMODULATION  DISTORTION  IS  '.FV.I,**.*) 

M 


v 
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r 


I 


k 


M’p;:ririi>;  q 

In  I  nr  O'-.:' ;  1  :i 1  i  c.r  r~  ■  -io:''  ion  v  • .  Tr.r  -i'..  :•! 


program  oinoot  i^uT,cvmrr,TMt«-<xrr»vT,puyn 
c- - interhod  v*.  :*vr  power - 

C  THIS  PROCRAM  MEASURES  I NTt*“C DELATION  DISTORTION  N«  A  FUNCTION  or 

C  Input  strrAi  pc  ••!  in  a  c.  o  % '  .~ZA  :*-  i  t,<?  r.:jcc;:.'»  rw9  trccrai 

C  PARAMETERS  A*E  F-LFECTLV  PlTO'D.  C '  LU,  LA  Tic  hi  *..£  ^..EaTED  Th»SX 
C  TINES  NS  THE  STEP  SIZE  RATIO  13  WftSIEJ  F*.JN  3J  09  TO  _j  Pi. 

c  rx  n«ti  sets  or  dpt*  wt  then  plotted  on  the  sane  graph. 

cnmiimunumntt  umi *ius  omnmmmmitrmnonmm 

c  X  •  AN  ARRAY  COKTAINIwQ  the  TEST  SIGNAL  INPUT  POWER  AT  WHICH  NEA- 
C  SURENENTS  ARE  PACE. 

C  V  ■  AN  ARRAY  CONTAIN  I  NO  THE  INTERNODULATION  OISTOtmON  HEASURt- 
C  NEWTS. 

C  a  •  AN  ARRAY  CONTAINING  TX  OUTPUT  FILTER  COEFFICIENTS. 

C  I  NOD  •  INTERNODULATION  DISTORTION  AT  THE  PRESENT  TEST  SIGNAL  POWER 

C  FS  •  THE  SAHPLE  RATE 

CSS  THE  FOLLOWING  UARIAILES  USE  AN  ■£■  AND  A  'D*  ("“"Ft*  TO  INDICATE 
C  USE  SY  EITHER  THE  ENCODER  CR  DECODES.  6ESPECTIULLY. 

C  FC1  •  ROLL-OFF  FREQUENCY  OF  THE  PRINARY  INTEGRATOR 

C  TC  •  THE  TIRE  CONSTANT  OF  THE  SVLLASIC  FILTER 

C  UMX  •  THE  NANINUH  INPUT  TO  TX  SVLLASIC  FILTER 

€  UNN  •  THE  NININUN  INPUT  TO  THE  SVLLASIC  FILTER 

C  SETA  -  THE  NORMALIZES  CENTER  FREQUENCY  OF  THE  OUTPUT  FILTER  TKANS- 
C  IT10N  IANS 

C  GANNA  •  THE  NORMALIZED  WIDTH  IF  THE  OUTPUT  FILTER  TRANSITION  SAND 

cmuimwwnmta  susroutines  used  umumnmmmmmm 

C  FLTRQCN  •  T»C  OUTPUT  FILTER  COEFFICIENT  GENERATOR 

0  FILTER  •  THE  sursouri*  THAT  FILTERS  THE  DECODER  OUTPUT  SIGNAL 
?  USING  THE  COEFFICIENTS  CALCULATED  |Y  FLTWGEN. 

:  IRIA-KOPT  .  CALCULATES  OTAR  AND  UNIN  FOR  THE  CWSD  ENCODE* 

C  AND  DECODER  SUBROUTINES. 

C  INTCRND  •  THE  SUBROUTINE  THAT  CALCULATES  THf  INTERNODULATION  DIS- 

C  T  OAT  I  ON  IN  T>£  CV3D  SYSTEN  OUTPUT  f  ZONAL. 

C  FACTOR.  PLOT.  SCALE.  AXIS,  RECT,  LINE  •  CALCCNP  PLOTTER  ROUTINES 

cssssmimssssssssssssssstsssssssssssssssssssssssssssssssssssssssosss 

C - - PROGRAM  START - - 

C -  INITIALIZE  VARAISLES  AND  ARRAYS 

DIMENSION  X(1SA).  VCIMI,  HIM) 

REAL  INOD 
I  CHAR  .  -I 

C - READ  AMS  PRINT  TIC  WORKING  UARIA1LES 

READ  S.  FS 
READ  S.  EFCl.  ETC 
READ  S.  BFC1.  DTC 
READ  >,  SETA,  CaNNA 
PRINT  *.*  F9  •  \FS 

PRINT  I,'  EFCl  •  '.EFCl.*.  ETC  •  *.CTC 
PRINT  «.*  Df Cl  •  \D-C1»*.  DTC  ■  • ,BTC 

PRINT  I,'  ICTA  •  *,irrA.*.  6ANNA  •  ' .GANNA 

C~—  GENERATE  LOU-PASS  FILTER  COEFFICIENTS 
CALL  FLTRGEMSETA.QAMNA.W.S) 

C -  START  CALCULATION  LOOP 

so  som  j  •  t.s.a 


IPO 


r.r 


calculate  0*9  system 


ICHM  •  ICHAA  *  I 
EMTIO  •  5*.  ♦  J 

coll  vmxorricvNx.cvm.FS.EFCl.fTC.ERRTloi 
1M  CCMT1HLZ 

DRATIO  *  ~+.  *  J 

Cou.  vnaxopt(Dl*«.9u»w.fs.dfci.d,’C.9Ratio> 

6 - calculate  intianodulation  distortion  vs.  nwr  pouch 


DO 


I  I  • 

-4*. 


UN 

♦  .41 


1 


X(  I  >  ♦  «W1 

call  intermits,  anpi.efci. etc, EVMX.iuMN.DPCi.Drre.oww.DuNH.NP.i.iA 
10D> 

vi  i  >  •  moo 

CONTINUE 


C - PLOT  T>€  MSULTP 

IP  ( I  CHOP  ,QT.  01  00  TO  DM 
CALL  FfcCTCRI .5) 

CALL  PLOT<3.,a.,-3> 
call  scALEix,i«..ic»,n 

CALL  SCALEIY.6..1- J.l> 

CALL  AXlsIO.. *..!••  U  VTCYTION  (*» .  14. S. .00. .VC101 OVUM) I 

call  rxisi«..s.,i..-.  vutue  a.M>,-i«.iA..*.,xUui>,x<iW)) 

CALL  A?CT(S.,#.,g.,H..4.,2) 

MO  CONTI  ML’S 

CALL  LI'  EtX.V.lM.l.lS.lCNAR) 

MM  CONTIKLi 

CALL  PLOTE(N) 

END  _ 

SUBROUTINE  INTEANDIFS.ANPl.LFCt.CTC.EUNX.EVNH.DF-Cl.DTC.lXMW.DlMW.N 
tP.D.lnODI 


C— — - INTEEMODULATION  01STCRTI0N  SUSROUTII* - 

C  TNI*  PROGRAM  CALCULATE*  TV.  INTTLYenULATION  DISTCTTICY  OF  A 
C  0*0  SVSTSN  »••••:.£  tri  f-fC  — 9  S  ■!  I  X.  I*  fit  CL  .  TCTLO  JAOK-TO- 
c  BACK.  DX0TC4TJ ;x  is  r  .•  x .  .o  »v  i* -urr:.-.i  A  r.  JT  .ll  cc-.tc50 
C  OP  TWO  e.'.'.tl  e  ,'LXtl-.-  j  at  i;„v  ri  *  t  t-  rx.  5;  : 

C  AMPLITUU  Of  Tv?  DI?v  :*?'*:£  F  C  >'JCT  10  T-'l  T.  > "X-  3  CC.PAAES 

C  TO  THE  INPUT  SIC.NAL  TO  L*T«0:«Ih£  THE  F6XUHT  DI-TC-irtON. 


IMUT  •  AM  INTEGER  ARRAY  CONTAINING  THE  BINARY  OUTPUT  OP  TtC  0*9 
ENCODER. 

OUTPUT  •  A  PEAL  A*"»AV  CONTMHJNO  TNE  TIP?  FUNCTION  OUTPUT  OP  TVC 
TEST  SIC-.AL  Cl:  .  ATCA  AND  AFTER  PROCCSflNG.  THE  TIf*  FUNCTION 
OUTPUT  OF  TNE  UCCStA. 

pyx  •  A  REAL  #-°A*  CCNTAINIPQ  TV*  OUTPUT  SPECTRAL  POL'S  CCMPOMPFT* 
OP  TV?  t  :cc  A  OU-  -  jT  ALT’-*  AtfOCESilNO  IY  THE  FAST  P0UH1CR 
TRANSFORM  <FFT  1  SL..VOUTIN*. 

I  UK.  UK.  CUK  •  U0RK1NG  ARRAY*  USES  IY  T>*  PPT  SUSR0UT1NE. 

FREttUC  •  A  REAL  ASSAY  CONTAINING  Twf  FBTOUENCIE*  AT  UNICN  TNE  PPT 
MM  CALCULATED  THE  SPECTXAL  C0APCNCN1S. 

C  ARPl  ♦  T>*  AMPLITUDE*  OP  THE  TEST  SIGNALS  IN  DBMB. 

<  P*  •  TTC  SAMPLE  RATS  IN  IP*. 

C  PEAK  •  THE  PEAK  VALUES  OP  TNI  TEST  SIGNAL  COMPOFCNT*. 

C  •  •**  OP  **•  PREFIX  OM  TNE  NTXT  TV»TE  SET*  CP  VARIABLE  INDICATES 
C  Ttc  VARIABLE  IS  U4XD  »Y  EITHER  THE  DECODER  OR  CNC0CCR.FX6P£C> 

«  TIVCLV 

C  PCI.  PCI,  FC3  •  POLl-CPT  FR'OUFNCIES  FOR  TW  PRINCIPU  INTEGRATOR 
C  IN  THE  0*9  C/ICCDER  AMS  KCCCER. 

C  PC  •  TX  COMPASSING  SPEED  OP  T>€  SYLLA9IC  INTEGRATOR  IN  SEC. 

«  RATIO  •  TNE  MAXIMUM  STEP  SIZE  TO  MIN!**  STEP  SIZE  RATIO  IN  Cl 

C  P  •  A  PREOUCNCY  VARIABLE  USED  TO  rETTWINE  T>*  TEST  *ICNAt  SPEC- 
0  TRAL  CC  VCMNT  AMO  ALSO  TO  DETERMINE  TNE  HARMONIC  SPECTRAL 

C  COMPONENTS. 

C  ACT  •  TV*  OUTPUT  SPECTRAL  COMPOKNT  POUCR  IN  99M. 

C  SUN  ■  TNE  OW  OP  TVC  POUER  IN  TNE  TEST  SIGNAL  COMCMWTS. 
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c  •  the  Nina  or  filter  coefficients 

C  •  •  Ml  ARRAY  CONTAINING  THE  INRUT  MO  OUTPUT  FILTER  COCFTICIOfTS 
ammnmnmmn  mmjriKi  u«a  ttmmmmmromana] 
e  D«C  DC  1  >  THE  CVSD  ENCODER 

e  sc  cocci  ■  nc  cvsd  decoder 

c  rvnwacn  «  generator  or  nc  output  filter  coefficients 

filter  •  filters  nc  output  using  nc  cccfficiewts  produced  *y 
nrnocn 

c  iimm.  •  nc  test  sional  generator 

C  FTFPS  •  nc  FAST  FCURICR  TRANSFORM  FROM  TV*  IBSL  LIMARY 


I 


c . . . SUSROUTINE  START - - - 

v - INITIALIZE  UARIMLES  AND  ARRAYS 

DIMENSION  I NPUt ! S9#A >, OUTPUT! 5**A),PSX! IBS) 

1 .  lUKl  cl  > ,  W(  ( 151  >,F,.5QU£  (cAA  >,  Ki/A) 

REAL  l NOD 
COMPLEX  CVKOtt) 

RiosN#)  •  sut  as.  n  news  -  4.»/ts.)  s  .ssi  t  sos.i  s  semi.) 

C -  DETERMINE  PEAK  VALUE  OF  TEST  SIONAL 


PEAK  ■  A  (ANPil 

C -  COCRATE  INPUT  TINE  FUNCTION 

CALL  SIGNAL (OUTPUT.  MAS.  FS.  TOO..  IASS..  PEAK,  PEAK) 
C -  PROCESS  THE  TINE  FUNCTION  TMWUGH  THE  CUSS  SYSTEM 


CALL  CNCOCEt (OUTPUT. IlWT.!P-'»A.FS.tFCl.FC3.FC3.ETC,EUre«.EUnn. DC) 
CAU.  DECODE!  ( InfVT.  OUTPUT.  6WR,  Ft,  OF  Cl .  Fit.  FC3.  PTC,  DWU.  DWSM.DC  > 

C -  FILTER  THE  OUTPUT 

CALL  FILTER (OUTPUT, SARA, NP, S’ 

c —  surntACT  the  aucrage  ualue  from  the  sanpu  string 

SUN  •  S. 

DO  C  I  ■  1.4F3S 
SUN  •  SUN  ♦  OUTPUT! I ) 

•  CONTINUE 

aucr  •  sum  /  4<n « 

PO  7  I  •  1,4  :;4 
OUTPUT! 1 )  •  OUTPUT! I)  -  AUER 
T  CONTINUE 

0 -  CALCULATE  THE  SPECTRAL  CONPONENTS  OF  THE  OUTPUT 

CALL  FTFPS  (OUTPUT,  DUN,  4AOS.KS.S.PSX,  SUN.  SUN,  tUK.tAC.CUC.lCR) 

«  OETERMIl*  THE  COMPONENT  AT  TW  TEST  FREQUENCY  AND  CALCULATE 

«  THE  RNS  VOLTAGE. 


c— - 


If 


ft*  •  0. 

SO  ■  K  *  3.139 
F  •  ((K-DA.5S.)  ■  FS 
FRfgiC(K-l)  •  F 
PSXlK-l)  •  P$X(K)  • 

IF  t  «F  .NC.  751.1  .AND.  <F  .HE. 
Sun  •  Sun  ♦  s^t  (  p;x«>  i 
REF  •  1A.  I  AL0G1S(PSX!K)  /  SAA 
CONTINUE 


IGOO.Il  00  TO  • 
/.SOI) 


calculate  internosulation  distortion 


SIF  •  Kl.  _ 

DO  IS  I  •  1,133 

IF  (FRZOLC(t)  .ME.  DIF)  00  TO  If 
SIS  •  PSX(I) 

CONTINUE 

I NOD  •  SORT  (SIS)  /  SIN  S  ISO. 

!F  (INOD  .07.  ISO.)  I NOD  ■  ISO. 

RETURN 

IMS 
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APPENDIX  R 


Sip,nal-to->!oi i.  .1io  Pro1  rii’. 

program  smm  an  <  input,  output,  tapes«output> 
i. - - — — filtepred  son  program - - - - -  ■ 

i  THIS  PROGRAN  MEASURES  T>€  SIONAL-TO-NOISE  PBfOWW  or  A  CUSD 
C  ENCODER  AND  DECODER  CONNECTED  >'i:CX-T0*?®C>C.  A  3T-;U  F'T'-EM CY 
C  SINE  WAGE  IS  I.-UT  TO  T,T  SYSTEM  AND  THE  DIFFERENCE  BETWEEN  TVC 
C  OUTPUT  «SD  INPUT  COWTED. 

C 

C  A  NAXZfWLLV  PLAT  LINEAR  PHASE  FIR  FILTER  IS  PLACED  Cft  TVC  OUTPUT 
C  Of  THE  DECODER  TO  REMOUE  SICNAL  CCMPOVCNTS  ABOUE  360*  HZ. 

c 

Clttniltlliummum  UARIABLES  »nrnnmHwmnnmnmn»ii 
C  fREOl  •  THE  TEST  SIGNAL  FREQUENCY 
C  SNR  •  T>€  MEASURED  SM  UAUJE 

C  TSIN  •  AN  ARRAY  CONTAINING  THE  INPUT  TINE  SERIES  SAfftES. 

TSOUT  •  AN  ARRAY  CONTAINING  FIFST  Tvr  CECCTSR  OUTPUT  TINE  SERIES 
SAMPUS.  THEN  THE  OUTPUT  Tlf.S  SERIES  SAMPLES  CP  TVC  PIR  FILTER. 

ERR  •  AN  ARRAY  CONTAINING  THE  DIFFERENCE  BETWEEN  TVC  OUTPUT  SAM¬ 
PLES  AND  THE  INPUT  SA.FLES. 

•  •  AN  ARRAY  CONTAINING  THE  FILTER  CCEPFICIENTS. 

•  INCUT  •  AN  ARRAY  CONTAINING  Tr€  UNARY  OUTPUT  Of  THE  CUSD  ENCODER 
ANP1  .  THE  ANPLITUDC  Of  THE  TEST  SIGNAL  IN  DM*. 
fS  •  THE  SAMPU  RATE. 

SCI.  fca.  FC3  •  THE  ROLL-OfP  FREQUENCIES  Of  THE  PRHlflflY  INTEGRA¬ 
TORS. 

TC  •  THE  TINE  CONSTANT  Of  THE  SYLLABIC  FILTERS. 

UMAX  i  UNIN  .  THE  NAXINUN  AND  MINIMUM  IWVTS  TO  THE  SYLLABIC  FIL¬ 
TER. 

BETA  •  TVC  NORMALIZED  3  DS  FREOUENCV  OF  THE  OUTPUT  FILTER.  T>C 
FREOl  ENCY  IS  KCNMALIZED  TO  THE  SANPU  RATE. 

OANNA  •  TVC  NOffTALITEO  WIDTH  C*  TVC  POLL-OFF  REGION  OF  TVC  OUTPUT 
FILTER.  TVC  F  I?N  IS  THE  FREQUENCY  BAND  BETWEEN  TVC  BSP  AT© 

R  OUTPUT  AMPLITUDES. 

PEAKt  •  TVC  NAXINUN  AMPLITUDE  OF  THE  TEST  SIONAL  IN  UOLTf. 

I#  •  TVC  NUNBCR  OF  FILTER  COEFFICIENTS . 

DC  •  The  duty  CYCU  of  THE  SLOPE  OVERLOAD  DETECTOR. 


itHtmititutnmmmtmmnmminmmnnnmnnimuu 

FLTRGEN  •  TVC  SUBROUTINE  THAT  GENERATES  TVC  OUTPUT  FILTER  COEFFI¬ 
CIENTS. 

SICNAL  •  TVC  TEST  SICNAL  GENERATOR.  PRODUCES  SAMPLES  OF  BINJSOI- 
DAL  WAULS  WITH  AT  POST  TWO  FflZOUENCY  COMPONENTS. 

ENCODE t  •  TVC  CUSD  ENCODER  SlT»0UTINE  WITH  A  SINOU  ROLL-OFF 
FREQUENCY  IN  THE  PHIPPRY  INTEMNTCR. 

DECODE!  •  THE  CUSD  CECODJNO  SLSrOUTIHE  WITH  A  SINOU  ROLL-OFF 
FREQUENCY  IN  THE  PRIIWRY  INTEGRATOR. 

FILTER  •  THE  SUBROUTINE  THAT  FILTERS  TK  INPUT  TI»S  SERIES  SAFFUS 
USINQ  THE  FILTER  CCEPFICIENTS  GENERATED  BY  FLTRGEN. 

POWER  •  A  ROUTINE  TO  CALCULATE  THE  POWER  IN  A  SAMPLED  Tim  SERIES 
WITH  INFERENCE  •  SOS  CM'S. 


cnmmwuoimmugmipuumuuunntuuimmamoni 


132 


c - procran  start - - - 

C -  INITIALIZE  UARIAILES  AKO  ARRAY* 

DIMENSION  TSIn<5«*#),  TSOCTIMM),  DIR ( MM ) 

1,11200) 

JNTEGtR  »I«0LrT(5*M) 

Aiosnai  •  sortu*.  nunnt  -a.  via.)  t  .aai  >  aaa.i  <  south.) 
c - iwvr  vokino  uariajles 

KAO  t,  FS 

kao  i.fFci.  rrc,  e^atio 

kao  a.t-ci.crc.::^T:o 

CALL  WV'XC^T tf.  •x,r.'-N,FS.tFCl.rrC,FRAT10) 

CALL  X.wL>«.FS.WCl,DTC,UlATIO) 

REAO  I.I£TA.  C  .  ,<A 
PEAX1  •  A(-<  i.  ) 

PRINT  *,  •  S‘R  TEST  AT  ESS  MZ  AND  SAPPD!  RATE  .  •  ,FS 

PRINT  *.  •  £FC1  •  •,^-Cl.*.  VZ  •  ,,'r0,',  C'VTIO  •  ‘.'RATIO 

print  *,  •  ffci  •  -C!.‘,  r*c  .  :  .atio  •  *,:  atz o 

PRINT  *.  *  FILTER  P>*AfCTE*S  ARE,  »£TA  •  *,KTA,  *,  QANNA  •  ‘.QAftl 


6- -  ODCRATE  OUTPUT  FILTER  COEFFICIENT* 

CALL  FLTRC£N<»tTA,OAflNA,)#,»> 

C - GENERATE  INPUT  TINE  SERIES  SATPLES 


CALL  5ICNAL<TS!N.Sa#*,F*,S#a. .PZAEl.t. > 

C— -  PROCESS  T)C  INPUT  TINE  SERIES  TVCOUOH  THE  CWSD  SYSTEA 

CALL  CNCOO€l(TSIN,tI‘OUT.f'M.FS.EFCl,F^,Fr!.E*C.Ejr^.D^T<.DC) 
CALL  BECOC£l<»lN0LrT.TSC^T.i,v^«,Fi.£FCl,FCa,F(.3.0;c,.W*<.0WT<.DC) 


C -  FILTER  TVC  OUTPUT  OF  TVC  CECCDCR 

CALL  FILTER(TSOUT.S«**.NP,l) 

C— —  DELAY  THE  INPUT  SIGNAL  START  TO  CORRESPOND  TO  TW 
C  FILTER*!!)  OUTPUT. 


M 


DO  M  ID  •  l.WI 

ro  -  £'#  ♦  is 

TS.N(ID)  •  TSIN(ICO) 
CONTINUE 


C -  ADJUST  OUTPUT  AMPLITUDE  SO  INPUT  POWER  •  OUTPUT  POWER 


CALL  P0UER(TSIN. ACTS. FS. PIN) 
CALL  POUERITSCU'T.a.  S.FS.POUT) 
OAIN  -  SwJT  (PIM,2CUT) 

DO  «  I  *  l.A.'-n 
TSOUT(I)  •  TSCUT <  1 )  S  CAIN 
4*  CONTINUE 

0 -  CALCULATE  TVC  NOISE  POWER 


f* 


DO  M  I  ■  1.A9M 

CRR( I )  ■  TSOUT(I)  -  TSIN(I) 

CONTIffJE 

CALL  P0UER(ERA,4M*,Ff,ERAP> 


CALCULATE  TVC  */* 


SNR  •  1*.  *  AL001S  (PIN  /  ERRP) 


PRINT  TVC  RESULTS 

PRINT  S,‘  TVC  S10D1AL  TO  NOIK  RATIO  • 
CONTIMJE 

oc 
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no  on  on  on 
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r.  ;  1  —  L?i  — *To  i  l‘~.1  i  i  •  ,  .  r!  *• 1  •  i  •ry~il  r,ri  ■ ;  ■  'icy  :'po  ~  r-1  v 


C 

c 

e 

c 


program  s»*rm«input, output, tapes«output, plots 
- FIirtMtJ  SNR  WOMM*— 


THIS  PROC#AH  NFASURES  TK  SIGNPL-VO-HOITE  PERFCRPANCt  OF  A  CVSD 

encoder  and  t£co: :a  con  *.cte5  :  '-cx-to-:--c<. _ «  *.>:••  die  f-'Scuehcv 

SIFC  WAVE  IS  Iv-UT  ’C  T  :  DySTLN  AMS  THE  DIFFERENCE  SETWEEN  THE 
output  and  input  cotv>ute5. 

A  fVWIWtLLV  FLAT  LINFA#  PHASE  FIR  FtLTfR  IS  PLACED  OH  TVC  OUTPUT 
Of  TVC  DECODE#  TO  RU/TOCE  SIGNAL  CCiVON-UNTS  A»CUC  36»R  HZ. 

THE  PROGRAM  IS  REPEATED  TVREE  TITES  AS  TV*  VALUE  OF  TVC  STEP  SIZE 
RATIO  IS  CVWNGED  FxUfl  32  TO  JS  DJ  IH  STEPS  Of  £  DJ. 


cimuitmtmtitmiu  uaiiaiui  ttxtntxxtmztuxuxnmxxxxtxtm 

fSTOl  •  Mi  ARRAY  CONTAINING  THE  FrtEOUEPCIES  THAT  THE  SNR  HAS  SEEN 
HCASUR1D  AT.  Th£  RANGE  IS  J.S  hS  TO  3UOO  HZ. 

SW  •  AH  ARRAY  CONTAINING  TVC  REASURED  SNR  UAUXS. 

TSIH  •  AH  ARRAY  CONTAINING  THE  INPUT  TIFE  FUNCTION  SAfPLES. 

TSOUT  •  AH  MTV  AY  CONTAINING  FI*ST  THE  DECCTER  OUTPUT  TIRE  FUNCTION 
SAHPLES.  TVC  THE  FIL72R  OUTPUT  TIT*  S'-^LES. 

ERR  •  AN  AHRAY  CONTAINING  TyE  DIFFERENCE  IE  TWEEN  TVC  OUTPUT  SAM¬ 
PLES  AND  THE  INPUT  SA'N»LES. 

■  .  AN  ARRAY  CONTAINING  TVC  FILTER  COEFFICIENTS. 

•INOUT  •  Mi  MIRAY  CONTAINING  THE  UNARY  OUTPUT  OF  THE  CVSD  ENCODER 
ANP1  .  THE  AfPLITUDE  OF  TVC  TEST  SIGNAL  IN  DINS. 

FS  .  TVC  SRPU  RATE. 

FC1.  FC2.  FC3  •  THE  ROLL-OFF  FREGUDCIES  OF  TVC  PRIIWRY  INTEGRA- 
TOW« 

C  TC  •  TVC  TINE  CONSTANT  OF  TVC  SYLLAlIC  FILTERS. 

C  UNAX  t  UNIN  •  TVC  HAXIHUH  AND  HINIHUH  INPUTS  TO  TVC  SYLLAlIC  FIL- 

C  TER. 

C  KTA  *  the  NOrr»LI7E5  3  Dl  FREQUENCY  Of  TVE  OUTPUT  FILTER.  THE 
C  FREQUENCY  IS  HJ^NALIZEO  TO  TVS  SAMPLE  RATE. 

C  MM  .  TV*  VT-'HPUZtD  WIDTH  Of  Tv*  ROLL-OFF  REGION  of  TVC  OUTPUT 
C  FILTER.  Tv*  >  I  .UN  IS  TVC  FREQUENCY  IAND  ItTVEEN  TVC  OS*  AND 

C  It  OUTPUT  A/*UTUDCS. 

e  PEARS  •  THE  HAXIMJH  A**LITUDE  Of  TVC  TEST  SIGNAL  IN  VOLTS. 

C  HP  •  TVC  NUHSCR  OF  FILTER  COEFFICIENTS. 

C  KM  •  THE  NUN1ER  OF  TEST  FREQUENCIES. 

C  1C  •  TVC  DUTY  CYCLE  OF  TVC  SLOPE  OVERLOAD  DETECTOR. 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


iiiiuittuntittiiitiuiiitiiuutmmtmiimsiitstmnitfimnn 
CutUHiiiitiiMiiii  suinout ires  vviD  iiimimuuiimuaiiiumu 

0  fLTRCEN  •  TVC  SUIROUTIVC  THAT  GENERATES  THE  OUTPUT  FILTER  COEFFI- 
•*  CIENTS. 

C  PLOT.  SCALE.  AXIS.  RCCT,  LINE.  PLOTE,  •  CALCO**  PLOTTING  ROUTI>CS. 

SIGNAL  •  TVC  TEST  SIGNAL  GENERATOR .  PRODUCES  ‘APPLES  OF  SINUSOI¬ 
DAL  WAVES  WITH  AT  POST  TWO  FREQUENCY  CCNAONENTS. 

ENCODCl  •  TVC  CVSD  ENCODER  SUDRCUTIHE  WITH  A  SINGLE  ROLL-OFF 
FREOUENCY  IN  THE  PHITARY  INTEGRATOR. 

It CO DC 1  •  TVC  CVSD  CCCC3INQ  SUT ROUTINE  WITH  A  (INGLE  ROLL-OFF 
FREOUENCY  IN  THE  PRIMARY  INTEGRATOR. 

FILTER  •  TVC  SUIROUTINE  THA*  FIL*rRS  TVC  IN*UT  TINE  fUYTION  SAN- 
PUS  USING  THE  FILTER  COEFFICIENTS  CiN*  RATED  IY  FLTRCEN. 

e  POWER  •  A  ROUTINE  TO  CALCULATE  TrC  POWER  IN  A  SAMPLED  TINE  FUNC- 
C  TION  WITH  INVEDENCE  •  M  OhN*. 


1  'Ll 


c - ppccpaa  start - 

C -  INITIALIZE  UARIASLES  AND  ARRAYS 

DINENSION  FVFQlUeS),  TSJN(BMS),  TSOUT(SMS),  ODi(UM) 

i,S(2?9>.  GAIrUc?*) 

INTEGER  »IfcXIT<ECeS) 

A!DSN*>  •  SCa!T<l».  M((ESH*  -4.  )/lf.  )  I  .Ml  t  CM.)  t  SORT <3. ) 

I CHAR  •  -1 

C -  IMVT  AND  PRINT  WORKING  UARIASLES 

READ  *.  ANP1,  FS 
READ  S.FC1 ,  TO 
READ  I. SETA,  CPMAA 
PEAK1  •  A(A/P1) 

PRINT  *.  *  s>:«  TEST  AT  \AflPl«*  DIM  AND  •,F»,‘  IPS* 

PRINT  I.  •  WITH  TC  ■  '.TO 

PRINT  *.  *  SETA  •  ‘.SETA.*.  SANAA  »  ’ , OAHPA 

C -  COCRATE  FILTER  COEFFICIENTS 

CALL  FLTRC€N(»CTA.OAPHA,NP.»> 

C -  IN1T1AL2E  PLOTTER 

CALL  FACTOR!  .$) 

CALL  PLOTO..  a..  -3) 

C -  START  OF  SICNAL-TO-NOISE  LOOP 

DO  1MO  NT1NCS  •  a.S.i 

KN  •  • 

10,4)  •  ICHAR  «  1 

RATIO  •  30.  ♦  Nrirrs 

.•mLL  UNAXOPT <OAAX. UNIN. FS.FC1.TC. RATIO) 

KWH  •  3M.36M.1M 
KN  «  AN  ♦  1 

C -  GENERATE  TEST  S1QNAL  FREQUENCY 

FREQUKN)  •  K 

C- -  GENERATE  INPUT  TIRE  FUNCTION  SARPLES 

CALL  SICNALCTSIN.SOM.FS.FREQl  (KN),S,  .PEAK1.8. ) 

C -  PROCESS  THE  input  rim  FUNCTION  throuch  the  cusd  SYSTER 

CALL  ENC0DEl(TSIN.lIN0UT.S^»,FS,FCl.FC2,FC3.TC,LroX.WHIN.DC) 

CALL  D£C00£l(SlNCUT,TS0UT.S0o#,FS,FCl,FCa.FC3.TC.UfiAX.WNIN,CC) 

C -  FILTER  THE  OUTPUT  OF  THE  DECODER 

CALL  FILTER<TSOUT,6M*.NP,») 

C -  DELAY  THE  INPUT  SIGNAL  START  TO  CORRESPOND  TO  THE 

C  FILTERRED  OUTPUT. 

DO  M  ID  •  1,4«M 
KD  •  2M  ♦  ID 
TSIN(ID)  •  TSINdCD ) 

N  CONTINUE 

C -  ADJUST  OUTPUT  AMPLITUDE  SO  INPUT  POWER  ■  OUTPUT  POUCR 

CALL  P0UER<T$:n,4MS.FS.P!N> 

CALL  POLER(TSOUT,4. .s.FS.POUT) 

OAIN(KN)  •  SO  IT  (PJh/VOUT) 

DO  4«  I  •  l,4tfJS 

TSOUT(I)  •  TSOUT(I)  S  GAIN(KN) 

40  CONTINUE 

©—  CALCULATE  THE  NOISE  POWER 

DO  M  I  •  1,4AM 

CRR(I)  •  TSOUT(I)  -  TSIN(I) 

M  CONTINUE 

CALL  POWER! ERR, 40M, PS. ERRP) 

C -  CALCULATE  THE  SYH 

•Ml (04)  ■  1«.  «  AL0C1S  (PIN  /  ERRP) 


1% 


f 


nor  the  «ajm 


M 

INI 


CONTINUE 

IF  (ICK.-R  .OT.  •)  00  TO  CM 
CPU  PlOT(13..B.,-3> 

CPU  tC:Cr.L(F-'3l,l*..«« 

CPU  CCFLtlt'  V*.6.,»LN.l  ) 

cpu  u,wis(».,t.,i«*+«:aucNCv  <«>, 

IN*2  •  l 

cpu  Axiii*.,*..sMs m  (ci).a.c.,M., 

CpU  P£CT<O.,0.. 6. .10. 

CONTINUE 


-14>l*..*..FHCai(DHI).FRE01(K 

*m(KN»n,«NH(OfH)) 


CpU  IGUN£(FRE01,S»«.ICN.I*,ICHP«,-1> 

CONTINUE 

CPU  PLOTEIN) 

CNO 


* 
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H  grnritchcrt  E'i  ■  n.-i I-io-i'Ioi ^  Palin  ts,  In-  u!  i 

PROGRAM  fWS^(U^. OUTPUT, TAPES-OUTPUT. PLOT) 


C- . .  — . MISNATCHED  6NR  PROGRAM - - - 

0  THIS  PROGRAM  MEASURES  TVS  SIGNAL-TO-NOIS5  PERFORMANCE  OF  A  CUSD 

o  encoder  and  cican  :s  coasted  b;csc-to-!a«.  a  sivsu  (-tzquency 

C  SINE  WAVE  IS  l  .-PUT  TO  7T  :  SYSTEM  AND  T«  DIFFERENCE  LETWLEN  THE 
0  OUTPUT  AM®  INPUT  COMPUTED. 

C  THE  CALCULATIONS  APE  PERFORMED  AS  THE  UAUJE  OF  THE  ENCODER  STEP 
C  SIZE  RATIO  IS  KEPT  CC:-  ;T.*NT  AND  THE  CECCl.'R  STEP  SIZE  RATIO  IS 

C  CHANCED.  THE  ENCODED  RATIO  IS  02  Cl.  WHILE  THE  DEC Cu£R  IS  STEPPED 

C  FROM  32  TO  32  DB  IN  STEPS  CF  a  Cl. 

c 

C  A  MAXIMALLY  FLAT  LINEAR  PHASE  FIR  FILTER  IS  PLACED  ON  THE  OUTPUT 
C  OF  THE  DECODER  TO  KETOUE  SIGNAL  COMPONENTS  A*  CUE  2SC*  KZ. 

C 

cwnniimimanan  uariajles  immmmimtnummwmin 

C  FRE01  •  AM  P.PRAV  CONTAINING  THE  nXOUEHCIES  THAT  THE  SNR  HAS  SEEK 
C  MEASURED  AT.  THE  RANGE  IS  SCO  KZ  TO  2.>3»  KZ. 

C  SMt  •  AM  ARRAY  CONTAIN  I MC  THE  MEASURED  CMS  UAUXS. 

C  TSIM  •  AM  ARRAY  CONTAINING  THE  INPUT  Tlfl  SERIES  SAMPLES. 

C  TSOUT  ■  AM  ARRAY  CONTAINING  FIRST  THE  DECODER  OUTPUT  TirE  SERIES 
C  SAMPLES.  THEN  THE  OUTPUT  TlhE  SERIES  SAMPLES  CF  THE  FIR  FILTER. 

C  ERR  *  AN  ARRAY  CONTAINING  TOE  DIFFERENCE  BETWEEN  THE  OUTPUT  BAM- 
C  res  AND  TIC  INPUT  SAMPLES. 

C  I  -  AM  ARRAY  CONTAINING  THE  FILTER  COEFFICIENTS. 

C  TIME  •  AM  ARRAY  CONTAINING  THE  TINE  THAT  THE  FIRST  20*  SAMPLES 
C  ARE  TAKEN  SO  THAT  THEY  MAY  BE  PLOTTED. 

C  BINOUT  «  AN  ARRAY  CONTAINING  THE  BINARY  OUTPUT  OF  THE  CUSD  ENCODER 

«  AfPl  •  THE  AMPLITUDE  OF  THE  TEST  SIGNAL  IN  DBM*. 

€  Ft  •  THE  SAMPLE  RATE. 

C  FCi#  FC2.  FC3  •  THE  ROLL-OFF  FREQUENCIES  OF  T«  PRIMARY  IMTEGRA- 
C  TORS. 

C  TC  •  THE  TIME  CONSTANT  OF  THE  SYLLABIC  FILTERS. 

C  UMAX  B  UMIN  •  THE  MAXIMUM  AND  MINIMUM  INPUTS  TO  THE  SYLLABIC  FIL- 

€  TER. 

C  BETA  •  THE  NORMALIZED  3  D*  FREQUENCY  OF  THE  OUTPUT  FILTER.  THE 
C  FREQUENCY  IS  NORMALIZED  TO  THE  SAMPLE  RATE. 

C  QAMMA  •  THE  NORMALIZED  WIDTH  OF  T>S  ROLL-CFF  REGION  OF  THE  OUTPUT 
C  FILTER.  Tt-I  >  ■  , A  IB  THE  FREQUENCY  BAND  BETWEEN  THE  SS*  AMD 

C  (>  OUTPUT  AMPLITUDES. 


C  HP  •  the  NUMBER  OF  FILTER  COEFFICIENTS. 

KM  •  THE  NUMBER  OF  TEST  FREQUENCIES. 

C  DC  •  nc  DUTY  CYCLE  OF  THE  SLOPE  OVERLOAD  DETECTOR. 


{iumii<miuutimmntmtmmwn<utitMunmnitn»;iitt 
Ctmtmmimmu  subroutines  USED  uu»itn»muim:muii:ui 

C  FLTRCCN  •  THE  SUBROUTINE  THAT  GENERATES  THE  OUTPUT  FILTER  COEFFI- 
C  CIENTS. 

C  PLOT.  SCALE.  AXIS,  RCCT,  LINE,  PLOTE.  •  CALCOMP  PLOTTING  ROUTINES. 

C  SIGNAL  •  THE  TEST  SIGNAL  GENERATOR.  PRODUCES  SAMPLES  OF  SIHUSCI- 

C  DAL  WAUES  WITH  AT  MOST  TWO  FREQUENCY  COMPONENTS. 

C  CHCOXI  •  THE  CUSD  ENCODER  SUBPCUTIPC  WITH  A  SINGLE  ROLL-OFF 
C  FREQUENCY  IN  THE  PH  If ARY  INTEGRATOR. 


C  DECODE 1  •  THE  CUSD  CECOOINQ  SUBROUTINE  WITH  A  SINGLE  ROLL-OFF 
C  FREQUENCY  IN  T>€  PRIMARY  INTEGRATOR. 
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filter  ■  tv?  sir^wrrvr  -vat  filttps  rvr  zypur  nr?  series  s**ples 

USING  THE  FILTER  CCCFFIG1ENTS  LOERATED  8Y  FLTSIEN. 

POWER  •  A  ROUTIfr:  TO  CCLCL'LATE  TV£  POWER  IN  A  SAMPLED  TINE  SERIES 
WITH  IKPEEENCE  •  604  CHNS. 

PL TINE  •  A  LINEAR  PLOTTING  ROUTINE  TO  PLOT  SIGNAL  F^PUTUDE  OS. 
TINE. 

cmttntuinnnmmimmuuiuiuuiiistmniinnnuiiuuMi 

C - PROCRAH  START - - 

C -  INITIALIZE  UARIA8LES  AND  ARRAYS 

DIFENSICM  FPE01MMI,  rWR(lT»>,  TSINIMOS).  TSOUT (6*0*).  ERRtSDM) 
1.8(203).  TI?  :<£  ')>.  CAIH(EG3) 
lNTECER  8Ii;C'JT(fO*) 

AIDSW)  •  SCRT(1*.  M((C8N*  -4.  VI*.  >  8  .Ml  t  6**.  >  *  SGRTC8. ) 

I CHAR  •  -1 


INPUT  WORKING  WWIAJLES 


READ  t.  AflPl,  FS 
READ  t.FCl,  TO 
READ  *.E£Tn.  C.VTVA 
PEAK1  •  A(F.'Pl) 

PRINT  t,  •  EM  TEST  AT  •,APP1.'  CTN*  AND  *,FS,* 
PRINT  *.  •  WITH  TO  ■  '.TC.1  C.'  D  .'<3TI0  *  •.RATIO 
PRINT  *,  •  I £TA  •  ‘.SETA,*,  CASNA  •  •.GA.WA 


IPS* 


C -  GENERATE  OUTPUT  FILTER  COEFFICIENT* 


CALL  FLTR0EN<*ETA.CAflNA.W.8) 
C -  INITIALZE  PLOTTER 


CALI  FACTOR  I. S) 

call  plot < 2. ,  a..  -31 

C— —  START  of  sicnal-to-noisc  LOOP 

DO  1K»  NTIHES  •  3.6,9 
KN  •  C 

ICHAR  •  ICHAR  ♦  1 

RATIO  •  30.  ♦  NTirES 

CALL  WMXCPT(UP.AX,  WIN.  FS.FCl.TC, RATIO) 

IF  (ICKAS  .GT.  •!  CO  TO  3 
tu»W  •  TO* 

EUNN  *  WIN 
3  CONTINUE 

DO  394  K  •  369,3044, It* 

KN  •  KN  ♦  1 

C -  GENERATE  TEST  SIGNAL  FREQUENCY 

FREOi(CN)  •  K 

C -  GENERATE  If*>UT  TINE  SERIES  SAHPLES 

CALL  SICNAl(TSIN,SM*,FS,FRE01  (KH  ),*. ,PFA<1,*.  ) 

C -  PROCESS  THE  IffUT  TINE  SERIES  THROUGH  THE  CUSD  SVSTEN 

CALL  ENCCDC1 (TSIN, 8INCUT,5?A3,PS, FC1 ,FCT, FC3, TC,£UHX,EVNN,DC  > 
CALL  DECODEKIINOUT.TSOUT ,  SCGt.FS.FCI  ,FC3,FC3,TC,  UNAX.WIN, DC ) 

C -  FILTER  THE  OUTPUT  OF  THE  DECODER 

CALL  FILTER < TSOUT, S*M, > 

C -  DELAY  THE  INPUT  SIGNAL  START  TO  COfiRESPW®  TO  T»C 

C  FILTERRED  OUTPUT. 

DO  3*  ID  «  1,4*96 

ICD  •  28*  ♦  ID 

TSIHUD)  •  TSIHOCD) 

3*  CONTINUE 

C -  ADJUST  OUTPUT  AMPLITUDE  $0  INPUT  POWER  •  OUTPUT  POWER 

CALL  POWER  ( TS  IN,  4604,  FS,  PIN) 

CALL  P0LER(Tr'VT,<?)6,FS.PCVr) 

GAIN(ITN)  •  S'  (T  (PINxACUT) 

DO  4*  I  •  I,-  '■* 

TSOUT(I)  •  TSOUT(I)  <  CAIN(KN) 

4*  COrCTIMJE 


CALCULATE  THE  NOISE  POUCH 


M 


DO  5#  I  •  1.4979 

EKRd)  •  Tsocrru)  -  tsinid 

CONTIKT 

CAU.  FX*JZR(EJW.4®««,F*,EJa»P> 


0. - CALCULATE  THE  */-N 


5NH<KNl  •  19.  t  ALO019  (PIN  /  ESHP) 
0 -  POINT  AND  PLOT  THE  RESULTS 


««e 

309 


999 

1999 


u0ITE(6,€9«)  FRECKtH),  S.NO(KH>,  OAINIOI) 
F0PHAT(?0X,FI9.1,!>X.F18.«,i;x,F19.9> 

CONTIS 

IF  < I CHAP  .CT.  8)  CO  TO  S04 
CALL  LCSCAL(FnC1.10,.C() 

CALL  SCAL£<Sf«,6.,CN,  1) 

.£?kV.lCA*Ii<0,'*"1^SEOteWCV  (KZ>,-14,19.,9.,FR€ai(Xmt  ),FU£01(IC 

InTC J  I 

<**•}■  JSI*55-*J"8HS'W  99.. 6HMICN*l».SNR(m**>» 

CALL  R'CT(8.,8.,6..tA.,9.,3) 

CONTIh-.E 

CALL  LCLIMt  (FRE9I ,  STS,  ICN,  19,  XCHAS.-l  ) 

CONTINLC 
CALL  FLOTE(N) 

END 
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APPENDIX  ’J 


Sir  no  1 -to-iToi '■■■■  jp  y~  T'V,P:  ''P  'i  i1  Pov-sr  Dro."r. i- 

PROCPAN  SNRC<  I  NPVT,OUPUT,TA«l«  OUTPUT.  PLOT) 

. - . SNR  US.  INPUT  POWER  PROGRAM - 

C  THIS  PROGAM1  MEASURES  THE  SICNAL-TO-NOISE  performance  op  p  cuso 

•'  ENCODER  AND  CCCC'iR  CONf  ECTtD  i'ACK-TO-S'-CK.  A  *U  FRECUENCV 

C  SINE  UAUC  IS  P--UT  TO  T-  SYSTEM  PND  THE  DIFFERENCE  EETWEEN  THE 
C  OUTPUT  AND  INPUT  COMPUTED. 

C  A  MAXIMALLY  FUAT  LINEAR  PHASE  FIR  FILTER  IS  PLACED  CM  THE  OUTPUT 
C  OF  THE  DECODER  TO  REMOVE  SIGNAL  COMPONENTS  AS  OLE  DGW  KZ. 

C  TW  INPUT  POWER  IS  U*fl!ED  FRCP  -A*  DDM#  TO  •  DT«9.  TVE  PROGRAM 

C  IS  REPEATED  AS  THE  STEP  SIDE  RATIO  IS  VARIED  FRCP  32  Cl  TO  3o  D». 

C  T>C  RESULTS  ARE  THEN  PLOTTED  ON  A  SINGLE  GRAPH. 


cunimnmtimimn  wwiaius  mumnmtmmmmammn 

C  FREQ1  •  AN  ARRAY  CCNTAIN1MQ  THE  FREOUEHCIES  THAT  TVS  SNR  HAS  SEEM 
C  MEASURED  AT.  THE  fo\nGE  IS  3DR  hZ  TO  HZ. 

C  SNR  •  Mi  ARRAY  CONTAINING  THE  MEASURED  SNR  UAUJES. 

C  TSIN  •  Ml  ARRAY  CONTR1HINO  THE  INPUT  TIME  SERIES  SAMPLES. 

C  TSOUT  •  Ml  ARRAY  CONTAINING  FIRST  THE  DECODER  OUTPUT  TIME  SERIES 
C  SAflPLES,  THEN  THE  OUTPUT  TITLE  SERIES  STAPLES  OF  THE  FIR  FILTER. 

C  ERR  •  Ml  ARRAY  CONTAINING  THE  DIFFERENCE  irruEEN  THE  OUTPUT  SAM- 
C  PUS  AND  THE  INPUT  SA/FLES. 

C  »  •  AN  ARRAY  CONTAINING  THE  FILTER  COEFFICIENTS. 

C  TIME  •  AN  ARRAY  CONTAINING  THE  TIME  THAT  THE  FIRST  ZSG  SAMPLES 
C  ARE  TAKEN  SO  THAT  THEY  NAY  JE  FLOTTED. 

C  SINOUT  ■  AN  ARRAY  CONTAINING  THE  UNARY  OUTPUT  OF  THE  CUSD  ENCODER 

C  ArtPl  .  THE  AMPLITUDE  OF  THE  TEST  SIGNAL  IN  DSMG. 

C  FS  •  THE  SAMPLE  RATE. 

C  FCt.  FCZ,  FC3  •  THE  ROLL-OFF  FREQUENCIES  OF  THE  PRIMARY  IMTEGRA- 
C  TORS. 

C  TC  •  THE  TIME  CONSTANT  OF  THE  SYUASIC  FILTERS. 

C  UMAX  l  UNIN  •  THE  MAXIMUM  AND  MINIMUM  INPUTS  TO  THE  SYUASIC  FIl- 

C  TER. 

C  SETA  •  T»C  MIDPOINT  OF  THE  OUTPUT  FILTER  TRANSITION  SAND. 

e  OAHMA  •  THE  NORMALIZED  WIDTH  OF  Tvr  POLL-OFF  REGION  OF  THE  OUTPUT 

C  FILTER.  TV-  REGION  IS  THE  FREOUENCY  SAND  IETWEEN  THE  Six  AND 

C  SX  OUTPUT  AMPLITUDES. 

C  PEAK!  •  TV*  MAXIMUM  AMPLITUDE  OF  THE  TEST  SIGNAL  IN  UOLTS. 

C  MP  •  THE  NUMSCR  OF  FILTER  COEFFICIENTS. 

C  OC  •  T>«  DUTY  CVCU  OF  TVC  SLOPE  OLCRLOAD  DETECTOR. 

Ciiumttutumtitmtitimtmnmmunimmtimnmmnm 
ouiiitiiiimtmiii  susRGWTiiseS  used  uiiua»»unuju3uumua 

0  FITRGEN  •  THE  SUSROUriNE  THAT  GENERATE*  THE  OUTPUT  FILTER  CCEFFI- 
C  CIEMTS. 

v  PLOT.  SCALE,  AXIS.  RECT,  LINE.  PLOTt.  •  CALCOMP  PLOTTING  ROUTINES. 

C  SIGNAL  •  TK  TEST  SIGNAL  GENERATOR.  PRODUCES  SAMPLES  OF  SIMUSOI- 

C  DAL  UAWES  WITH  AT  MC9T  TWO  FREOUENCY  COMPONENTS. 

C  ENCODE!  •  THE  CVSD  ENCODER  SUSROUTIHE 

C  DECO  DEI  •  THE  CVSD  DECODER  SUJROUT 

C  FILTER  •  THE  SUSROUTIHE  THAT  FI  IT  ,/S  THE  INPUT  TIDE  SERIES  SAMPLES 
C  USING  THE  FILTER  CCEFFICIEJfTS  GENERATED  »Y  FlTr.aiJi. 

C  POWER  •  A  ROUTIK  TO  CALCULATE  THE  POWER  IN  A  8AFFUD  TIME  SERIES 
C  WITH  IMPECENCE  •  SR*  ChmS. 

CtittimituittiutaunmtmmiuntniniuisumnnmuiMut 
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C - PROGRAM  START - 

C -  INITIALIZE  UAflIAiLES  AND  ARRAYS 


DIMENSION  AMPl (SOS ),  SNR(2M>,  TSIM(E«M>,  T60UT(SMS>.  EERCSMS) 
l-SCc^l 

IN1ECER  SINOUTISCM) 

A(DIT4)  •  SGRTUS.  XX  (DSN*  -A.  }/!•.)  t  .Ml  S  CM.)  S  SORT<2. ) 

I  CHAR  •  -1 

C -  INPUT  AND  PRINT  THE  WORK I NO  VARIABLES 

READ  *,  FREOl,  FS 
READ  S.FCl,  TC 
READ  t.BETA.  CAN** 

PRINT  *.  *  SNR  TEST  AT  ■ , FREOl , *  DEM*  AND  ‘.FS.‘  IPS* 

PRINT  *,  *  WITH  TC  •  •  ,TC.  *  3  RATIO  •  ‘.RATIO 

PRINT  *.*  FILTER  PARAMETERS  ARE,  JET A  •  META.'.GAflMA  •  '.GANNA 

0 -  OEHERATE  OUTPUT  FILTER  CCEFFIC1ENTS 

CALL  FLTRCENIBETA,  CAWA.NP.S) 

C -  INITIAL2E  PLOTTER 

CALL  FACTOR <.S> 

CALL  P1.0T ( 2 • .  2..  -3) 

C -  START  LOOP 

DO  lfrM  rn  •  8. A. 2 
RATIO  •  S3.  ♦  r.T 

CALL  VMAXOPTIUFIAX.  WHIN.  FS.FCl.TC. RATIO) 

C- -  START  OF  SIGNAL-TO-NOISE  LOOP 

DO  300  K  •  MM 
ANPHKI  •  -4*.  ♦  .4  *  K 
PEAK!  •  A(ANPKK)) 


0— —  GENERATE  1HVT  TINE  SERIES  SAMPLES 

CALL  SIQNAL(T$IH,5fM,FS,FR£Ql,6.  ,PEAK1,S. ) 

C -  PROCESS  THE  INPUT  TIME  SERIES  THROUGH  THE  CVSD  SYSTEM 

CALL  ENCODE 1 < TS IN, B I N0UT,SNMS,FS.FCi,FC3.FC3,TC. UMAX, UNIN, DC) 

CALL  DeC0DEl<BIN0UT,TS0UT,5CvJ,F3,FCl,FC2,FC3.TC.U.'!AX.UmN.DC) 

C -  FILTER  THE  OUTPUT  OF  THE  DECODER 

CALL  FILTER(TSOUT,SMS.NP.I) 

C -  DELAY  ns  input  SIGNAL  START  TO  CORRESPOND  TO  THE 

C  FILTERREO  OUTPUT. 

DO  36  ID  •  1.4*96 
KD  •  200  ♦  ID 
TSIN(ID)  •  TSINUCS) 

36  CONTINUE 

C -  ADJUST  THE  OUTPUT  SIOMAL  AMPLITUDE  SO,  OUTPUT  POWER  •  INPUT  POWER 

CALL  POUCRtTSIN,  4006, FS, PIN) 

CALL  POL^»<TSOUT,4.;4,FS,POUT) 

GAIN  •  SORT  (PIN/POUT) 

DO  4t  I  •  l,4J.i 
TSOUTCI )  •  TSOUT(I)  S  GAIN 
46  CONTINUE 

C -  CALCULATE  THE  NOISE  POWER 

DO  66  I  •  1,4*88 

ERR(I)  •  TSOUT(I)  -  T8IHCI ) 

66  CONTINUE 

CALL  POWER! ERR, 4M6,FS,ERRP) 

C -  CALCULATE  THE  %A\ 

SNR(K)  •  16.  8  ALOQIS  (PIN  /  ERRP) 

IF  <  ICHAR  ,LT.  t )  GO  TO  c«* 

IF  <SNR(K>  .LT.  SHNltlll  SH»(K>  ■  GMKltt) 

MO  CONTINUE 

3M  CONTINUE 
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PSINT  AMD  nor  Ttt  RESULTS 


ICHIW  •  !W«  ♦  l 
IF  (ICHW  .CT.  t)  CO  TO  990 
CALL  £CALC(f:'J»t.lS..l«#.l> 

call  scau<£.'*,s..::*.i> 

^ call  AXi8(s.,#.,iSKArtjLm®c  (Dsnt>.-is.is..s..AflPiusi>,MPi(iM> 

CAU  AXISct.,*.,8MSNR  ( D»  ) ,  t,«.  ,8f.  ,SM(  1*1 ),  SNRUtfl) ) 

CALL  RECTiS.. IS. .A. ,3) 

MS  CONTINUE 

:aLL  LINE ( AAPl  ,SHR,  IM,  1»  lt>ICHAR ) 

IMS  CONTINUE 

call  PLOTE(N) 

END 
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APPKND1X  V 

Misnatci'.c:  r,.i--n~)  1  -  tn-;lo  ;  :-il;o  v. ,  In:' 
Si  tk  1  ' 'o r  i  i-  >  r\- . 


PROGRAM  ICYSMDC  INPUT, OUTPUT, TAPES»OUTPUT .PLOT ) 

C - MISNATCXD  SNA  VS.  INPUT  POWER  PROGRAM - 

«'  THIS  PROCRAM  MEASURES  THE  SICNAL-TO-NOl  ft  PERFOTAHCE  CP  A  CVSD 

c  ehcodir  and  d--.cc:.  :a  covs  :cted  fack-to-i;  -:k.  a  c:  ;ls  ff-guehcv 

0  SINE  WAVE  IS  If/UT  TO  T-:  SYSTEM  PND  THE  DIFFERENCE  IETUEEN  TX 
C  OUTPUT  AND  INPUT  COMPUTED. 

tf 

C  P  fttXlNAUY  FLAT  LINTAR  PVAEE  FIR  FILTER  IS  PLACED  CN  THE  OUTPUT 
C  OF  TX  DECODER  TO  FENOuC  SIGNAL  CCVCXNTS  PiOUE  3C3C  HZ. 

C  TX  INPUT  POVER  IS  <JP  R I  ED  FROM  -4R  DIPS  TO  *  DO'S.  TX  PROWPH 

C  IS  REPEATED  PS  TX  STEP  SIZE  RATIO  IS  UPRIED  F  .M  22  DG  TO  ES  D» 

C  IN  TX  DECODER.  WHILE  THE  ENCODER  STEP  SIZE  RATIO  IS  VELD  CQN- 
C  STANT. 

C  THE  RESULTS  ARE  THEM  PLOTTED  OH  A  SINGLE  GRAPH. 

C 

cttmtsitmtmmtmt  variables  nuinnnnnnmttttnutnnm 

C  FREQl  •  AN  ARRAY  CONTAINING  TVS  FT!ECUZPCIE9  THAT  THE  SW  HAS  SEEN 
C  flEASURED  AT.  THE  RANGE  IS  SCO  HZ  TO  S.C»  HZ. 

C  SNA  ■  AN  ARRAY  CONTAINING  THE  MEASURED  SNR  UAUiES. 

C  TSIN  «  AN  ARRAY  CONTAINING  THE  INPUT  TINE  SERIES  SAMPLES. 

C  TSOUT  •  AN  ARRAY  CONTAINING  FIRST  TX  DECODER  OUTPUT  TirE  SFRIES 
C  SAMPLES,  THEN  TVS  OUTPUT  TIPS  SERIES  S-VPLSS  CF  THE  FIR  FILTER. 

C  ERR  •  AN  ARRAY  CONTAINING  THE  DIFFERENCE  IETUEEN  THE  OUTPUT  Sttt- 
C  PLES  AND  THE  INPUT  SATtPLES. 

C  >  •  AN  ARRAY  CONTAINING  THE  FILTER  COEFFICIENTS. 

C  TIME  •  PH  ARRAY  CONTAINING  THE  TIME  THAT  THE  FIRST  EC*  SAMPLES 
C  ARE  TAKEN  SO  THAT  THEY  MAY  BE  PLOTTED. 

C  BINOUT  •  AM  ARRAY  CONTAINING  TX  UNPRY  OUTPUT  OF  THE  CUSD  ENCODER 

C  AHP1  .  TX  AMPLITUDE  OF  THE  TEST  SIGNAL  IN  DIN*. 

C  FS  •  THE  SAMPLE  RATE. 

Ctt  NOTE*  A  *D*  OR  ’E1  PXTIX  ON  TX  FOLLOWING  VPRIASLES  INDICATES 
«  THAT  THE  W-RIR51E  IS  U'ZD  EY  TVS  TICCICR  CR  EHCC'ER.STIPSCTIUE- 

C  LY.  UMAX  AND  WIN  MAY  BE  CONTRACTED  TO  UHX  AND  VH1,  SESPECTIVE- 

«  LV. 

C  FCi,  FC3.  FC3  •  THE  ROLL-OFF  FREQUENCIES  OF  THE  PRIMARY  INTEGRA- 
C  TORS. 

C  TC  •  THE  TIME  CONSTANT  OF  THE  SYLLABIC  FILTERS. 

C  UMAX  I  UNIN  «  THE  MAXIMUM  AND  MINIMUM  INPUTS  TO  THE  SYLLABIC  FIL- 

C  TER. 

C  SETA  •  TVS  NORMALIZED  MIDPOINT  OF  THE  TRANSITION  SAND  OF  THE  OUT- 
C  PUT  LOU  PASS  FILTER. 

C  GAMMA  •  THE  NOPMALIZED  WIDTH  OF  THE  POLL -OFT  REGION  OF  TX  OUTPUT 
•*  FILTER.  THE  AEGICN  IS  THE  FREQUENCY  SAND  BETWEEN  TX  8S*  AMD 

C  S\  OUTPUT  AMPLITUDES. 

0  XXI  •  TX  MAXIMUM  AMPLITUDE  OF  THE  TEST  SIGNAL  IN  UCLTS. 

C  HP  •  TX  NUMSER  OF  FILTER  COEFFICIENTS. 

C  KM  •  TX  NUMSER  OF  TEST  FREQUENCIES. 

C  DC  •  TX  DUTY  CYCLE  OF  TX  SLOPE  OVERLOAD  DETECTOR. 

cmuuiumtimitnmuimtuttniitimtnmnitmsmimnn 
Ciituiuoitiitniu  sucaouTixs  used  »mwm«xmuma«amu» 

C  FITRCEN  •  TX  SUWOUTINE  THAT  QEXRATES  TX  OUTPUT  FILTER  CCEFFI- 
C  CIENTS. 

C  PtOT.  SCALE,  AXIS,  RE CT,  LINE,  PLOTE.  •  CALCOTP  PLOTTING  ROUTINES. 

e  SIGNAL  •  TX  TEST  SIGNAL  GCXRATOR.  PRODUCES  SAMPLES  OF  SINUSOI- 

C  DAL  WAVES  WITH  AT  MOST  TWO  FREQUENCY  COHPOXNTS. 
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C  ENCODE!  •  TX  CVSO  ENCODER  SUIROUTIX 

C  DECODE!  •  THE  CVSO  DECODER  SUIROUTIX 

FILTER  •  THE  SUBROUTINE  THAT  F  HTTPS  THE  INPUT  TT*E  SERIES  SAMPLES 
USING  TX  FILTER  COEFFICIENTS  CONERATED  IV  FLTRSEN. 

POWER  •  A  ROUTINE  TO  CALCULATE  TX  POWER  IN  A  SAMPLED  TIME  SERIES 
WITH  IMPEECfaCE  •  Cii  CHtS. 

onnomtuuBmmimnnMttmmmnimmtHnmtmimq 

conHPMMKwaunitauunmnmninnBaamanmitma 

C - PROGRAM  START - . - 

C -  INITIALIZE  VARIASLES  AND  ARRAYS 

DIMENSION  AMPKZSS),  SNS(aM),  TSIN(SMS),  TSOUNSSM),  CRRCSMS) 
i.I  ) 

INTECER  IINOUTCSIM) 

AIOLMA)  •  &-iXT<  It.  XSUDINS  -A.)/1S.)  t  .Ml  t  644. )  t  SORTIE.  > 

I  CHAR  •  -1 

C -  IXUT  AMO  PRINT  WORKING  WaRIAILES 

READ  t,  FRE01.  PS 
READ  *.FC1,  TC 
READ  t.LETA.  CV*A 

print  *.  •  r-i  T.E3T  at  \frmi,*  dim  AMD  *,F8» *  IPS* 

PRINT  *.*  WITH  TC  •  * . *C 

PRINT  *.•  FILTER  PARAMETERS  ARE.  BETA  .  • . SETA.  •  .QM*M  •  ‘.OAfTM 

C -  GENERATE  OUTPUT  FILTER  COEFFICIENTS 

CALL  FlTRQENdCTA.QAfVM.NPa> > 

C - -  INITI ALZE  PLOTTER 

•’ALU  FACTORi  .Si 
CALL  PL0T12..  2..  -31 

C - START  LOOP 

oo  teas  nr  •  a.s.a 

RATIO  •  3».  ♦  IP 

CALL  VNAXOPTl UMAX, VAIN. FS.FCl.TC, RATIO) 

IF  «  ICHX  .C£.  •>  00  TO  S 
CUKX  •  UMAX 
Cum  •  VAIN 
•  CONTINUE 

c -  START  OF  SIQNAL-TO-NOISE  LOOP 

H3MK  •  1.1M 
AMPMK)  •  -44.  ♦.AIK 
PEAK1  •  A(AVIU) ) 


C -  CEICRATE  INPUT  TIME  SERIES  SAMPLES 

CPLL  SICHALITSIN, SMS. FS.FREQl.Q., PEAK1, >. ) 

C -  PROCESS  TX  INPUT  TIME  SERIES  THROUGH  TX  CVSO  SYSTEM 

CALL  ENCODFIITSIN, UNCUT, SF0A,FS,FC1.FC3.FC3,TC,EWX.EUNN. DC) 
CALL  DCCOD€ldINOUT,TSOUT.Se*»,FS,FCt.FCa,FC3,TC,V«AX,VHIN.DC) 

C -  FILTER  TX  OUTPUT  OF  TX  DECODER 

CALL  FILTERITSOUT.WM.NP.I) 

C -  DELAY  TX  INPUT  SIGNAL  START  TO  CORRESPOND  TO  TX 

C  FILTERRED  OUTPUT. 

00  30  ID  •  1.4S9S 
KP  •  EM  ♦  ID 
TSINI ID )  •  TSINIK9) 

30  CONTINUE 

C -  ADJUST  OUTPUT  SIQNAL  AfV>LITUSC  SO.  OUTPUT  POWER  •  IXUT  POWER 

CALL  POWERITSIH, 4S9S.FS.PIN  > 

CALL  POUUMTSOUT,<;;S,F8.PCUT) 

QAIH  •  S«T  (PIN/PCUT) 

DO  AS  I  •  !.4;.S 
Tsourin  •  rsouTit)  ■  gain 
AS  CONTI** 
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c - calculate  ne  noi*c  pome* 

W  H  1  •  J.40I8 
EPRfll  •  TSOUTU)  -  TSIK(I) 
H  CONTINUE 

CALL  P<XCX<ER*.4M«.Ff,EXftA> 

C -  CALCULATE  TW  ts* 


SfMiKi  •  l«.  «  ALOO10  (PTN  /  CRAP) 

IF  » I CHAR  ,LT.  0)  CO  TO  f9 
M  •tT*  SNR< 101)1  S.NROC)  •  *Mt<101> 

••  cOHTIf<u€ 

3M  CONTINUE 

C - PRINT  ANO  PLOT  THE  RESULT* 

I  CHAR  •  I  CHAR  ♦  t 

IF  i  ICHAR  .CT.  ()  GO  TO  CM 

CALL  SCALE(/;0»t.l9..IC»,I) 

CALL  SC  ALE  ( £.  1*‘3,1) 

^CALL  AXISIO.  /♦. , i&vV^LITUDE  CCSH* IS.  1*. ,AAP1  ( 101  ),ANP1 C  1*2) 

cAUKECT!^:^:^!S^,3^•''9•••^a•1,'^,1•^,, 

000  CONTI  ALE 

CAU  UTE(A.TPlaSNRa100alal*aICW«) 

1000  cottimje 

CALL  PLOTE(N) 
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APPKND1X  W 


K];it  i  c’.  :  >r-c;  liG  <  ■/  or,  1  ro-j-in 


PROGRAM  DIFGAIN<  INPUT, OUTPUT, TAPE6*0UTPUT, PLOT ) 

«. . . —DIFFERNTIAL  CAIN  PROGRAM - 

v  THIS  PROGRAM  MEASURES  THE  SVSTEM  CAIN  US.  FREQUENCY  RESPONSE  FOR 
•-  H  C.'SO  DIGITAL/AMALOG  SVSTEM  CONNECTED  BACK-TO-BACK. 

C  H  NmXIMAUV  flat  LINEAR  PHASE  fir  filter  is  placed  on  the  output 
C  OF  THE  DECODER  TO  REMOVE  SIGNAL  COMPONENTS  A30UE  3SB0  HZ. 

cuuMuminmimut  variables  tmmmmmununmuumu 

C  FREOl  •  AH  ARRAY  CCNTAININQ  THE  FREQUENCIES  THAT  THE  SHR  HAS  BEEN 
C  MEASURED  AT.  THE  RANGE  IS  300  HZ  TO  3300  K2. 

C  TSIN  •  AN  ARRAY  CONTAINING  THE  INPUT  TIME  FUNCTION  SAMPLES. 

C  TSOUT  •  AN  ARRAY  CONTAINING  FIRST  THE  DECODER  OU’H’UT  TIME  FUNCTION 
C  SAMPLES,  THEN  THE  OUTPUT  TIME  FUNCTION  SAMPLES  OF  THE  FIR  FIL- 

C  TER. 

C  B  -  AN  ARRAY  CONTAINING  THE  FILTER  COEFFICIENTS. 

C  BINOUT  •  AN  ARRAY  CONTAINING  THE  BINARY  OUTPUT  OP  THE  CUSD  ENCODER 

C  AMP1  •  THE  AMPLITUDE  OF  THE  TEST  SIGNAL  IN  DBM0. 


C  FS  •  THE  SAMPLE  RATE. 


C  FC1,  FC2,  FC3  «  THE  ROLL-OFF  FREQUENCIES  OF  THE  PRIMARY  INTECRA- 
C  TORS. 

C  TC  •  THE  TIME  CONSTANT  OP  THE  SYLLABIC  FILTERS. 

C  UMAX  l  VMIN  •  THE  MAXIMUM  AND  MINIMUM  INPUTS  TO  THE  SYLLABIC  FIL- 

C  TER. 

C  BETA  •  THE  NORMALIZED  CENTER  PREOUENCY  OP  THE  OUTPUT  FILTER  TRANS- 
C  ITION  BAND. 

C  GAMMA  -  THE  NORMALIZED  WIDTH  OF  THE  ROLL-OFF  REGION  OF  THE  OUTPUT 
C  FILTER.  TrE  KEGICN  IS  THE  FREQUENCY  BAND  BETWEEN  THE  9SX  AND 

C  5*  OUTPUT  AMPLITUDES. 

C  PEAK!  •  THE  MAXIMUM  AMPLITUDE  OF  THE  TEST  SIGNAL  IN  UOLTS. 

C  HP  •  THE  NUMBER  OF  FILTER  COEFFICIENTS. 

C  KM  •  THE  NUMBER  OF  TEST  FREQUENCIES. 


C  DC  •  THE  DUTY  CYCLE  OF  THE  SLOPE  OVERLOAD  DETECTOR. 


cutinu»unnnunminnn»ntmmm:>m!nm»mmmni 
Citmtmmutnni  subroutines  USED  xxxxxxxxxxxtxxxxxxxxxxxxxxxxxxxxx 


C 

A 

l 

c 


FLTRCEN  •  THE  SUBROUTINE  THAT  GENERATES  THE  OUTPUT  FILTER  COEFFI¬ 
CIENTS. 


»l"T,  SCALE,  AXIS,  RECT,  LINE,  PLOTE,  ■  CALCCMP  PLOTTING  ROUTINES. 

MGNAL  •  THE  TEST  SIGNAL  GENERATOR.  PRODUCES  SAMPLES  OF  SINUSOI¬ 
DAL  WAVES  WITH  AT  MOST  TWO  FREQUENCY  COMPONENTS. 


ENCODE!  •  THE  CUSD  ENCODER  SUBROUTINE 


V 

C 


DECODE 1  •  THE  CUSD  DECODER  SUBROUTINE 

FILTER  •  THE  SUBROUTINE  THAT  FILTERS  THE  INPUT  TIME  FUNCTION  SAM¬ 
PLES  USING  THE  FILTER  COEFFICIENTS  GENERATED  BY  FLTRGEN. 


C  POWER  ■  A  ROUTINE  TO  CALCULATE  ThE  POWER  IN  A  SAMPLED  TIME  FUNC- 
C  TIOH  WITH  IMPEDENCE  ■  600  OHMS. 


cttxttxxxxxxxitxttxsxixxxxtttttxttxixttttttxtxttxxttxxattxxxxtxtttttttxt 

ZtttttXXXXXXXXXiXXXXXXXXXXXXXXtlXttXXXXXXXXXXXttiXXXXXiXXXXXXXilXXXXtXXt 
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C - PROCRAN  START - 

C -  INITIALIZE  VARIABLES  AND  ARRAYS 

DIMENSION  FREO1C100),  TSIMSMO),  TSOUT(SOM) 
i,B<20O),  CA1  HCc‘:d) 

INTEGER  BINOUTC5C30 ) 

A<DBN0>  •  SGRTU0.  XXUDBNO  -4.)/10.)  X  .Ml  X  6M. )  X  SQRT<2.) 
I CHAR  •  -1 

C -  INPUT  AND  PRINT  THE  WORKING  VARIABLES 

READ  X,  AMP1 .  FS 
READ  X,FC1,  TC 
READ  X.BETA,  GAMMA 

PRINT  X.  ■  SNR  TEST  AT  '.ANPl,*  DBAS  AND  ‘.FS,*  BPS' 

PRINT  X,  •  WITH  TC  ■  \TC 

PRINT  X.*  OUTPUT  FILTER  PARAMETERS  ARE.  BETA  •  ‘.BETA.*  GANNA  • 
INNA 

C -  GENERATE  OUTPUT  FILTER  COEFFICIENTS 

CALL  FLTRGEN< BETA. GANNA. NP.BJ 

C -  INITIALZE  PLOTTER 

CALL  FACTOR (.5 1 
CALL  PL0T(2. .  2..  -3) 

C -  START  LOOP 


DO  2OO0  NTINFS  •  2.6,2 
ICHAR  ■  ICHA3  +  1 
RATIO  •  SO.  ♦  NTirES 
hLL  1 41A.NOPT  ( UMAX,  uniN.FS.FCl.TC,  RATIO) 

»:n  •  ■) 

* -  start  Of  SIGNAL-TO-NOISE  loop 

00  300  K  ■  300,3600.100 

KN  •  KN  ♦  1 

C -  GENERATE  TEST  SIGNAL 

FREOIOCN)  •  K 
PEAK1  •  A<A."?1) 

CALL  SlCNAL(TSIN,50e0,FS,FREQl (KN).O. .PEAK1.0. ) 

C— —  PROCESS  THE  INPUT  TINE  FUNCTION  THROUGH  TV€  CUSD  SYSTEN 

CALL  ENC0DE1 <TSIN,BINOUT,S330,FS,FC1 ,FC2,FC3,TC,VNAX,VNIN  DC) 
CALL  DECODE 1 ( BINCUT, TSCUT, 5C30,f S.FCl .FC2.FC3.TC, vnAX.UNIN, DC  > 

C -  FILTER  THE  OUTPUT  OF  THE  DECODER 

CALL  FILTER<TS0UT.5Cee,NP.B> 

c -  DELAY  THE  INPUT  SIGNAL  START  TO  CORRESPOND  TO  THE 

C  FILTERRED  OUTPUT. 

DO  30  ID  •  1.4096 

W>  •  200  ♦  ID 

TSIN(ID)  •  TSIN(KD) 

30  CONTINUE 

C -  FIND  SYSTEN  GAIN 

CALL  P0WER(TSIN,4OS6,FS,PIN) 

CALL  POUER(TSOUT,40SG,FS.PCUT> 

GAIN(KN)  •  10.  X  ALCGlOtPOUT/PIN) 

3M  CONTINUE 


.CA 
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■6  C -  ADJUST  CAIN  VALUES  70  88*  HZ  REFERERENCE 

DO  6  X  •  l.KN 

IF  (FREQKI )  .EC.  808.)  REFCAIN  •  GAIN(I) 

6  CONTINUE 

DO  8  I  •  1,KN 

CAIN(I)  •  GAINU)  -  REFCAIN 
S  CONTINUE 

C -  PLOT  THE  RESULTS 

IF  < I  CHAR  .GT.  0)  CO  TO  S88 
CALL  SCALE(CAIN.6.,ICN,t) 

CALL  LCCC«LCfAEJ1.13..'CN) 

CALL  LCAXISC0..O..1W^CUENCV  (HZ  >,-H,  18.  ,0.  ,FR£Q1  (ICN+1  ),FRECl(tC 
1N>2 ) ) 

CALL  AXIS(0./8..£2HDIFFERENTIAL  CAIN  <D»  >,22,6.,88.  .GAINOCN+l ), 
1CAINOCN+2 ) ) 

CALL  RECTO. ,0.. 6. ,18. .8. .3) 

988  CONTINUE 

CALL  LGLIN£(FR£Q1,CAIN,ICN,  10, 1  CHAR,  “1 ) 

2888  CONTINUE 

CALL  PLOTE(N) 

END 
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an  no  o  o  no  r.  o-,  •»  a  non  on  no  o  no  o  o  o  o  non  o  on  nnononn  o 


APrKMDiy  x 


’  ^  i  t  ' '  ^  1  .  i  t  ':V'  '  |'pr*  ■*' :  ]  ■  Mr- ,r  ■  r,  ‘'r>  '  *^o 

PROGRAM  MMDCA IN ( INPUT, OUTPUT ,TAPE6* OUTPUT, PLOT ) 

— . . MISMATCHED  SYSTEM  MIN  RESPONSE - 

THIS  PROCRAM  MEASURES  THE  SYSTEM  GAIN  US.  FREQUENCY  RESPONSE  FOR 

A  CVSD  DIGITAL/ ANALOG  SYSTEM  CONNECTED  BAOC-TO-SACK. 

A  MAXIMALLY  FLAT  LINEAR  PHASE  FIR  FILTER  IS  PLACED  ON  THE  OUTPUT 

OF  THE  DECODER  TO  REHCUE  SIGNAL  COMPONENTS  ABOVE  3S88  H2. 

ttitittmtttmmiiiit  variables  tmtmtmmimnntimmnm 

FRE01  •  AN  ARRAY  CONTAINING  THE  FREQUENCIES  THAT  THE  SNR  HAS  BEEN 
MEASURED  AT.  THE  RANGE  Is  3M  H2  TO  3=38  HZ. 

TSIN  •  AN  ARRAY  CONTAINING  THE  INPUT  TIME  FUNCTION  SAMPLES. 

TSOUT  •  AN  ARRAY  CONTAINING  FIRST  THE  DECODER  OUTPUT  TIME  FUNCTION 
SANPLES,  ThCN  THE  OUTPUT  Tlf.E  FUNCTION  SAMPLES  OF  THE  FIR  FIL¬ 
TER. 

•  •  AN  ARRAY  CONTAINING  THE  FILTER  COEFFICIENTS. 

•INOUT  •  AN  ARRAY  CONTAINING  THE  BINARY  OUTPUT  OF  THE  CUSD  ENCODER 

AflPl  -  THE  AMPLITUDE  OF  THE  TEST  SIGNAL  IN  DBMS. 

FS  •  THE  SAMPLE  RATE. 

FC1,  FC2.  FC3  •  THE  ROLL-OFF  FREQUENCIES  OF  THE  PRIMARY  INTEGRA¬ 
TORS. 

TC  -  THE  TIME  CONSTANT  OF  THE  SYLLABIC  FILTERS. 

UMAX  l  UNIN  •  THE  MAXIMUM  AND  MINIMUM  INPUTS  TO  THE  SYLLABIC  FIL¬ 
TER. 

BETA  -  THE  NORMALIZED  CENTER  FREQUENCY  Of  THE  OUTPUT  FILTER  TRANS¬ 
ITION  BAND. 

GAMMA  •  THE  NORMALIZED  UIDTM  OF  THE  ROLL-OFF  REGION  OF  THE  OUTPUT 
FILTER.  THE  REGION  IS  THE  FREQUENCY  BAND  BETWEEN  THE  95X  AND 
EX  OUTPUT  AMPLITUDES. 

PEAK1  •  THE  MAXIMUM  AMPLITUDE  OF  THE  TEST  SIGNAL  IN  UOLTS. 

HP  •  THE  NUMBER  OF  FILTER  COEFFICIENTS. 

f.il  •  THE  NUMBER  OF  TEST  FREQUENCIES. 

PC  •  THE  DUTY  CYCLE  OF  THE  SLOPE  OVERLOAD  DETECTOR. 

iiuuitmtuttniimsimmmtmummnmtnmtmiimmn 
Ktimitmumm  subroutines  used  KCJisiuiutsummmumi 

FLTRCEN  •  THE  SUBROUTINE  THAT  GENERATES  THE  OUTPUT  FILTER  COEFFI¬ 
CIENTS. 

PLOT,  SCALE,  AXIS,  RECT,  LINE,  PLOTE,  •  CALCOMP  PLOTTING  ROUTINES. 

SIGNAL  •  THE  TEST  SIGNAL  GENERATOR.  PRODUCES  SAMPLES  OF  SINUSOI¬ 
DAL  WAVES  UITH  AT  MOST  TWO  FREG.  ENCY  CCrPONENTS. 

ENCODE1  •  THE  CVSD  ENCODER  SUBROUTINE 

DECODE 1  •  THE  CUSD  DECODER  SUBROUTINE 

FILTER  •  THE  SUBROUTINE  THAT  FILTERS  THE  INPUT  TIME  FUNCTION  SAM¬ 
PLES  USING  THE  FILTER  COEFFICIENTS  GENERATED  BY  FLTRGEN. 

POWER  •  A  ROUTINE  TO  CALCULATE  THE  POWER  IN  A  SAMPLED  TIME  FUNC¬ 
TION  UITH  IMPEDEHCE  •  680  OHMS. 

Cxxxxtxxxxxtxxxxxxtxxxxxxxxxxxxxxxxxxxtxxxtxxxxxxxxxxxaxxxxxxxxxxxxxxtxx 

CimtmuKttttimimiiimiHitmmiiuumnmmmmunm 


- - PROGRAM  START - - - 

C -  INITIALIZE  VARIABLES  AND  ARRAYS 

.Di^!I0N  TsiNcseeo),  Tscurcseee) 

INTEGER  BlNOUTCSJdd) 

ICMAR3.5  :xSMJU<d‘  «“»»"•  t  .991  8  600.)  *  S0RT(2.  ) 

; -  INPUT  AND  PRINT  THE  WORKING  VARIABLES 

READ  *,  AMP1,  FS 
READ  t,F Cl.  TC 
READ  X.8ETA,  GF.-fSA 

mint  i: :  fnW?cT.‘'T-.,Tcm,l‘‘ Dsm  ™  tPS' 

IMa”1  ***  0UTPUT  FILTER  PARAMETERS  ARE.  BETA  •  '.BETA.1  GAMMA  *,CA 

: -  GENERATE  OUTPUT  FILTER  COEFFICIENTS 

CALL  FLTRGEN< BETA, GAMMA. NP.B) 

-  INITIALZE  PLOTTER 

*LL  FACTOR(  .S) 

•  hLl  puma.,  a..  -3) 

-  SThRT  loop 


do  2000  NTiriES  •  2.6.2 

ICHAR  •  ICHAR  ♦  1 
RATIO  •  30.  ♦  NTirSS 

CALL  UM AXOPT ( UMAX, VH IN, FS.FCl.TC, RATIO) 

kn  •  e 

IF  (ICHAR  .GT.  0)  CO  TO  199 
EUMX  •  UMAX 
EUMN  •  UMIN 
tea  CONTINUE 

C -  START  OF  SIGNAL-TO-NOISE  LOOP 

DO  399  1C  •  390,3606, tee 

ICN  •  KN  ♦  1 

C -  GENERATE  TEST  SIGNAL 

FREQl(ICN)  -  K 
PEAX1  -  A(AMPl) 

CALL  SIGNAL<TSIN.5000.FS,FREQI(ICN),e..PEAICl.e.) 

C -  PROCESS  THE  INPUT  TIME  FUNCTION  THROUGH  THE  CVSD  SYSTEM 

CALL  ENCODEKTSIN.BINOUT, 5OD0.FS ,FC1,FC3,FC3,TC,EVNX,EVMN,DC ) 
CALL  DEC0DE1  ( BIISCUT,  TSOUT,  5200,  FS,FC1»FC2,FC3,TC,  UMAX,  UMIN,  DC) 

C -  FILTER  THE  OUTPUT  OF  THE  DECODER 

CALL  FILTER(TSOUT,S200,NP,I) 

C -  DELAY  THE  INPUT  SIGNAL  START  TO  CORRESPOND  TO  THE 

C  FILTERRED  OUTPUT. 

DO  30  ID  •  1,4096 

KD  •  220  ♦  ID 

TSIN(ID)  •  TSIN(KD) 

30  CONTINUE 

C -  FIND  SYSTEM  GAIN 

CALL  POWER* TSIN, 4096, FS, PIN) 

CALL  POWER (TSOUT , *0SS,FS,POUT ) 

GAIN(KN)  •  10.  X  ALCG10(POUT/PIN) 

300  CONTINUE 

C -  ADJUST  CAIN  VALUES  TO  800  KZ  REFERERENCE 
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If  ^fREai(pN-£G»  80®.  )  REFGA1H  •  W1H(I) 

6  continue 

CA1N(J)  •1CAIN(I)  ~  REEGAIN 
8  CONTINUE 


c -  plot  THE  RESULTS 


w* 

■UAA 

CvW 


Arr'f'I  TV 


Tent.  Si  'nst  C'~  ;  ‘ r  ■  lor 


5UW0UTINE  SIGNAL!OUTPUT,N.rS,FREQt,FR£G«.A»l.A«»*) 
Citumuimuuiimmnmtmnimnnmmninmmmmnn 
C  THIS  5USP0UTINE  GENERATES  *  TEST  SICNAL  COMPOSED  OF  UP 

c  to  two  sine  uaues  of  differs  frequencies  and  amplitudes. 

cmiuuuiiutitimni  maihus  mtuimunimummiuinnn 
C  OUTPUT  •  AN  ARRAY  CONTAINING  THE  OUTPUT  TINE  FUNCTION  SMFUS 
C  N  •  TX£  HPUER  Of  SAMFOES  Of  THE  TINE  FUNCTION  DESIRED 

C  FS  •  THE  SWli  RATE  IN  K»/S 

C  FRE01  •  THE  FREOUENCY  Of  THC  FIRST  SIGNAL  COMPONENT 

C  FREQ2  •  THE  FREOUENCY  OF  TX  SECOND  SICF4U  C OPPONENT 

C  «SP1  •  THE  PEAK  AMPLITUDE  OF  THE  FIRST  SIGNAL  CCfPOFCfT 

C  MVZ  •  THE  PEAK  AMPLITUDE  OF  THE  SECOND  SIGNAL  COMPONENT 

CtumMtuutmuuiuusuuuuimnsxjxtmiwn'  usunisum 

C -  INITIALIZE  UAflIAJLES  AND  ARRAYS 

DIMENSION  OUTPUT! N) 

DATA  PlF3.1SlSSaSS3SF 

C -  GENERATE  OUTPUT  SAMPLES  Of  TEST  SIGNAL 

DO  S«  I  -  t.N 

OUTPUT! I)  -  ANP1  t  SIN!  2.  *  PI  *  FREQt  /  FS  I  1) 

!♦  FhPS  t  SIN  (2.  *  PI  *  FHEG2  /  FS  S  I) 

6*  CONTINUE 

RETURN 
END 
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OO 


Arn-rm;:  /. 


’ml 


1  il  •: !  i  or.  ’v. .ro  i  \  ■ 


SUSR0UT1«  POVER(X.N.FS.P) 

. --AUCRACE  POUER  SU**OUTI*. 


THIS  5U3PCVTINC  CALCULATES  THE  ATTRACT  PCt*R  or  THE  INPUT  TI?* 

fur*TiOH  s-wus  ac.- css  a  *:•>  c~n  ii .*  ♦:?«■;«£.  tv*  itmrr  tire  FLEC¬ 
TION  SAMPLES  A^E  1I*HJT  TO  TvE  iUIRCUTlAt  AS  A  1  X  N  A.-RAY. 


ciimmmimumtu  uttmtumiummnmmmnti 

C  X  •  AN  ARRAY  CONTAINING  THE  INPUT  TIKE  FUNCTION  SAMPLES 
C  N  ■  THE  HUNSER  Of  SA«»US  TO  IE  PROCESSED 
C  FS  •  THE  SAMPLE  RATE 
C  P  •  T>€  CALCULATED  SIGNAL  POWER 

cmmnumionnroioniomminntoaaoaaointamnon 

C - SUBROUTINE  START - 

C -  INITIALIZE  U ARIAS LES  AND  ARRAYS 

DIMENSION  X(N ) 

SUN  • 


C -  SUN  TVC  SQUARES  Of  T«  INPUT  TINE  FUNCTION  SAMPLES 

DO  IS  I  •  l.N 

sun  •  sun  ♦  xtn  «  a. 

IS  CONTINUE 

c -  CALCULATE  THE  AWE RACE  POWER  ACROSS  A  CSS  OHN  IWPEDENCE 

P  •  SUN  /  N  /  GSS. 

RETURN 

END 
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Item  PO  continued. 
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j n  the  icaft  standard,  The  model  is  then  exorcised  by  varyint  t!.e 
parameters  to  the  limits  imposed  by  the  standard  ana  the  resui'inp  per¬ 
formance  compared  to  the  previously  determined  idea]  system  performance 
The  results  show  that  the  performance  characteristics  measured  are  nosh 
sen:  : t i ve  to  the  primary  in terra tor  response  and  output  filter  response 
when  the  system  siarameters  are  restricted  to  the  ran, re  allowed  by  the 
‘draft  t'ATO  standard. 
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